Electrical measurement of cross-sectional position of particles flowing through a microchannel

The first high-throughput system for the electrical detection of cross-sectional position and velocity of individual particles flowing through a rectangular microchannel is presented. Lateral position (along channel width) and vertical position (along channel height) are measured using two different sets of coplanar electrodes. In particular, the ratio of travel times measured with electrodes generating a current flow transverse or oblique with respect to particle trajectory yields lateral position. The relative prominence and transit time of a bipolar double-Gaussian signal obtained with a suited electrode configuration, respectively, supply vertical position and velocity. The operating principle is presented by means of finite element numerical simulations. The method is experimentally validated by comparing the electrical estimates of position and velocity of polystyrene beads with optical estimates obtained by processing high-speed images. The system is used to observe bead focusing at different particle Reynolds numbers. This system, providing a fully electrical characterization of single-particle motion, represents a powerful tool, e.g. to understand fluid motion at the microscale, in particle separation studies, or to assess the performance of particle focusing devices. Moreover, it can be simultaneously used to perform single-cell impedance spectroscopy, thus achieving an unprecedented multiparamteric characterization.


Introduction
Microfluidic applications are continuing to expand in the fields of nanotechnology, manufacturing and life sciences, and therefore, there is a pressing need for methods to characterize fluid behaviour at the microscale and suspended particle motion (Williams et al. 2010;Winer et al. 2014;Tasadduq et al. 2015). As an example, detecting cross-sectional position of flowing particles is essential in the investigation of self-ordering mechanisms (Gao et al. 2017), inertial particle focusing (Di Carlo 2009) or active particle focusing [e.g. dielectrophoretic (Shaker et al. 2014) and acustophoretic (Muller et al. 2013)], which have applications in single-cell analysis and sorting (Wyatt Shields IV et al. 2015). In fact, under specific force fields, position can directly indicate the property of target cells (Wang et al. 2017).
Numerous optical methods have been developed to obtain flow field and particle position in microfluidic channels. Conventional video microscopy depicts particle motion projected onto a plane defined by the focus of the objective. For many applications (e.g. micromixers or particle focusing devices), visualization of motion in the third dimension is required to fully understand the trajectory of particles (Tasadduq et al. 2015). To that aim, multi-camera (e.g. stereoscopic imaging, tomographic imaging) and single camera (e.g. confocal scanning microscopy, anamorphic/astigmatic imaging, digital holographic microscopy, deconvolution and defocusing microscopy) approaches have been proposed [for a description and comparative analysis see e.g. the reviews Lee and Kim (2009), Williams et al. (2010) and Cierpka and Kähler (2012)]. However, those methods require expensive and intricate optical setups involving numerous lenses, aperture designs and light sources, or demand for refined 41 Page 2 of 13 and intensive processing of massive image data (Winer et al. 2014;Lin and Su 2016).
An attractive alternative to optical systems is represented by electrical impedance-based approaches (Wang et al. 2017;Brazey et al. 2018). Impedance-based microfluidic devices have been extensively studied in the last decade due to their simplicity and potential for low-cost and portable implementation. In particular, microfluidic impedance cytometers for label-free analysis of single particles and cells at high throughput (Cheung et al. 2010;Sun and Morgan 2010;Petchakup et al. 2017) have been used in different biological assays, including particle sizing and counting, cell phenotyping and disease diagnostics [e.g. Haandbaek et al. (2016), McGrath et al. (2017), Rohani et al. (2017), Esfandyarpour et al. (2017) and Rollo et al. (2017)]. Those systems use electrodes integrated in a microchannel to measure variation in the electric field caused by the transit of particles. The accuracy of the technique is challenged by the positional dependence issue, i.e. identical particles flowing in the microchannel along different trajectories provide different signals (Spencer and Morgan 2011). This is due to fringing electric field and may introduce blurring in the estimated particle properties. Compensation strategies aimed at solving the positional dependence issue have been proposed by our group, thus achieving high accuracy particle sizing. To this end, new metrics have been devised, correlating with either lateral (i.e. along channel width) (Caselli et al. , 2018 or vertical (i.e. along channel height) (Spencer et al. 2016;De Ninno et al. 2017) particle position.
Those metrics are here exploited for quantitative particle localization. The first impedance-based device for high-throughput determination of cross-sectional (i.e. both lateral and vertical) position and velocity of single particles flowing through a rectangular microchannel is presented. The system comprises two consecutive regions, each equipped with a different set of coplanar electrodes (Fig. 1). As detailed in Sect. 2, (i) an electrical estimate, X, of lateral particle coordinate (x-coordinate) is obtained from the ratio R of the travel times measured using two sets of electrodes placed in lateral channels (first region, Fig. 1b Fig. 1 a Microscope image of the electrical sensing area. The region highlighted with a yellow box, comprising five pairs of lateral electrodes, is used for the determination of particle lateral position. The region highlighted with a blue box, formed by five electrodes spanning the channel width, is used for the determination of particle vertical position and velocity. b, c 3D rendered models of the two meas-uring regions. Lateral position is determined by comparing signals relevant to transverse and oblique current flows (yellow and green field lines in panel (b), respectively). Vertical position and velocity are determined exploiting an electrode configuration that generates a peculiar current distribution (blue field lines), characterized by alternating high-and low-field regions (colour figure online) (iii) particle velocity is extracted from the latter signal, dividing the spatial length of the sensing zone by the transit time.
A thorough validation of the proposed approach is reported. Impedance signals for 6-and 10-μ m polystyrene beads have been collected. The electrical estimates of particle cross-sectional position and velocity are compared with the estimates derived by image processing of optical frames recorded using a high-speed camera. Moreover, the method is used to monitor the bead flow under different flow rates, and partial inertial focusing is observed.
The paper is organized as follows: the operating principle enabling the electrical estimate of particle cross-sectional position and velocity is presented in Sect. 2, the experimental methods are described in Sect. 3, and validation and results are presented in Sect. 4. Finally, conclusions are drawn.

Electrical estimate of particle lateral position
The chip region used for the electrical estimation of particle lateral position (i.e. x-coordinate) exploits the concept of liquid electrodes (Demierre et al. 2007). They are vertical equipotential surfaces generated on apertures of the main channel walls by recessed coplanar electrodes patterned at the bottom of dead-end lateral chambers (Fig. 1b). Liquid electrodes provide a nearly homogeneous electric field over the total channel height while keeping a simple process flow with a single planar metal layer (Mernier et al. 2012).
In the present work, five pairs of coplanar electrodes are considered and the following wiring scheme is proposed ( Fig. 2a): an AC voltage is applied to two stimulation electrodes on one side of the main channel (the 2nd and the 4th electrodes along z, at x > 0 ), and two differential current signals are measured using four sensing electrodes located on the opposite side (i.e. at x < 0 ). The first signal, denoted I TSV , is the difference in electric current flowing through the sensing electrodes with the same z-coordinate as the stimulation electrodes (the 2nd and the 4th electrodes along z) which features a current path transverse to particle flow. The second signal, denoted I OBQ , is the difference in electric current flowing through two sensing electrodes with a z-coordinate different from the stimulating electrodes (the 1st and the 5th electrodes along z) featuring a current path oblique with respect to particle flow. An analogous wiring scheme was introduced in Spencer et al.  and an oblique measurement ( I OBQ ). b, c FEM simulations of the differential current signals (real part) generated by a flowing particle (6-μm-diameter insulating bead). b When a particle travels near the stimulation electrodes (i.e. in the positive x region, trajectory 1) the peak-to-peak distances Δz TSV and Δz OBQ are similar and their ratio is close to one. c When a particle travels near the measuring electrodes (i.e. in the negative x region, trajectory 2) Δz OBQ is greater than Δz TSV and their ratio is greater than one. d Simulated peak-to-peak time ratio R (equivalent to Δz OBQ ∕Δz TSV ) versus particle x-coordinate. The one-to-one relationship between R and x shows that the ratio R can be used to estimate particle lateral position Figure 2b [resp. c] shows the trace (real part) obtained in simulation when a dielectric bead travels along trajectory 1 [resp. 2] depicted in Fig. 2a. The differential signals I TSV and I OBQ exhibit a bipolar profile, and the relevant peak-to-peak distances are denoted by Δz TSV and Δz OBQ , respectively. While Δz TSV is the same for both trajectories, Δz OBQ depends on particle lateral position: it is higher for particles travelling near the sensing electrodes (trajectory 2, Fig. 2c) than for particles travelling near the stimulation electrodes (trajectory 1, Fig. 2b). From the experimental point of view, the peakto-peak times Δt TSV and Δt OBQ are available, instead of the peak-to-peak distances Δz TSV and Δz OBQ . Although Δt TSV and Δt OBQ depend on particle velocity, their ratio coinciding with Δz OBQ ∕Δz TSV , is a dimensionless quantity, independent of velocity, that correlates with particle lateral position. Consequently, R can be used to obtain an electrical estimate X of particle x-coordinate.
A simple geometric argument, assuming straight current paths ( Fig. 2a), yields the following linear estimate: where W denotes the channel width. Simulation results ( Fig. 2d) confirm that such a simple model is valid, except for a slight saturation at extremal lateral positions. A refined model, involving quadratic and cubic corrections, is introduced to cover the whole position range: The calibration procedure of the parameters b 1 and b 2 is described in Sect. 4.3.

Electrical estimate of particle vertical position and particle velocity
In order to estimate particle vertical position (i.e. y-coordinate), the sensing zone comprising five electrodes spanning the main channel width (Fig. 1c) is used. The relevant wiring scheme, introduced in De Ninno et al. (2017), is depicted in Fig. 3a: an AC voltage signal is applied to the central electrode, and the difference in electric current flowing through the lateral electrodes is measured. Intermediate electrodes are left floating. This wiring scheme results in a non-homogeneous electric field distribution along the channel axis (z-direction), characterized by weak-field regions in front of the floating electrodes, which in turn are reflected as local minima in the collected current. Therefore, the resulting signal trace exhibits a bipolar double-Gaussian (BDG) profile and is denoted as I BDG .
(1) R = Δt OBQ ∕Δt TSV , Electrical metric encoding particle vertical position. a Diagram showing the electrode design and the wiring scheme. Compared to the standard three-electrode measurement configuration [e.g. Gawad et al. (2001)], two additional floating electrodes are inserted between the stimulation and the measuring electrodes. These electrodes create a non-homogeneous electric field distribution, shunting particles travelling close to them. b, c FEM simulations of the differential current signal (real part) generated by a flowing particle (6-μm-diameter insulating bead). The signal exhibits bipolar double-Gaussian (BDG) shape which is characterized by the relative prominence P of the peaks compared to the saddle. b When a particle travels far from the electrodes (i.e. in the positive y region, trajectory 1) the relative prominence P is low. c When a particle travels near the electrodes (i.e. in the negative y region, trajectory 2) the relative prominence P is high. d Simulated relative prominence P versus particle y-coordinate. The one-to-one relationship between P and y shows that the relative prominence P can be used to estimate particle vertical position Figure 3b [resp. c] shows the trace (real part) obtained in simulation when a dielectric bead travels along trajectory 1 [resp. 2] depicted in Fig. 3a. The prominence of the two peaks with respect to the saddle in between is higher for particles travelling close to the electrodes (trajectory 2, Fig. 3c) than for particles travelling away from the electrodes (trajectory 1, Fig. 3b). Because signal amplitude also depends on particle size, the following normalized metric, referred to as relative prominence, was introduced in De Ninno et al. (2017): where m and M are the signal amplitude at the saddle and peaks, respectively.
Finite element simulations (Fig. 3d) show that the relative prominence P correlates with the height of particle trajectory (y-coordinate): the higher the former, the lower the latter. This claim has been experimentally supported in De Ninno et al. (2017) by means of a quantitative defocusing approach (Wu et al. 2005). As a consequence, the relative prominence P can be used to obtain an electrical estimate Y of particle y-coordinate. In particular, a quadratic model is used in this work: where H denotes channel height. As shown in , the parameters c i depend on the experimental setup (e.g. buffer conductivity, frequency of the AC stimulation, electrode double-layer capacitance). The relevant calibration procedure is described in Sect. 4.3.
An electrical estimate V of particle velocity is obtained from the signal I BDG as follows: where Δz BDG is the distance between the centres of the floating electrodes ( Δz BDG = 80 μ m for the present geometry) and Δt BDG is the relevant transit time ).

Microfluidic chip
The impedance chip consists of a glass microscope slide with integrated gold microelectrodes, bonded to a PDMSembedded microchannel (Fig. 1). The device was fabricated following standard microfabrication techniques, as reported elsewhere ). The main channel was 40 μ m wide and 21.5 μ m high. The lateral channels in the lateral position determination region were 30 μ m wide and 30 μ m apart from each other, with electrodes recessed by 20 μ m with respect to the main channel. The five electrodes in the vertical position determination region, spanning the main channel width, were 30 μ m wide with a 10-μ m gap between them. For fluidic access, Teflon tubing (OD: 1/16″) was inserted into chip inlet and connected to a 500-μ l glass syringe (Hamilton) loaded on a syringe pump (Elite 11, Harvard Apparatus). The chip was connected electrically to a custom chip holder using pogo pins.

Sample preparation
The system was tested using polystyrene beads with diameters of 6 or 10 μ m (Sigma-Aldrich). The beads were resuspended to a concentration of approximately 10 6 beads/ ml in PBS containing 0.1% Tween 20 and enough sucrose to match particle density (1050 kg/m 3 ). The conductivity of the final medium was 1.1 S/m. The samples were sonicated prior to the experiments, in order to prevent bead aggregation. Samples were pumped through the device at a flow rate of 3, 10 or 30 μl/min.

Impedance data acquisition and processing
Electrical signals were measured using two transimpedance amplifiers (HF2TA, Zurich Instruments) and an impedance spectroscope (HF2IS, Zurich Instruments). Details of the electrical connections are reported in Supplementary material, Figure S1. Excitation signals of 8 V, 1 MHz and 6 V, 884.8 kHz were used for the determination of lateral and vertical positions, respectively. Differential currents were sampled at 115.2 kSa/s, with 6, 20 or 60 kHz filter bandwidth, respectively, for 3, 10 or 30 μl/min flow rate.
Event detection in the data streams was performed with a previously reported algorithm , and a simple MATLAB script was used for event feature extraction by template fitting (Spencer et al. 2016;De Ninno et al. 2017). In particular, for each detected event the ratio R (Eq. 1) was computed from the peak-to-peak times of the oblique and transverse signals I OBQ and I TSV , whereas the relative prominence P (Eq. 4) and the transit time Δt BDG were extracted from the bipolar double-Gaussian signal I BDG . Parameters R , P and Δt BDG were then used to compute electrical position X, electrical height Y and electrical velocity V by means of Eqs. (3), (5) and (6), respectively.

Optical data acquisition and processing
For validation purposes, independent estimates of particle cross-sectional position and velocity were obtained by means of an optical approach. The optical recording region was located downstream the electrical sensing region. Images of a 1024 × 112 pixels area around the channel were recorded using a high-speed camera (Photron FASTCAM Mini UX100, 4000 fps, 4 μ s shutter time) mounted on an inverted microscope (Zeiss Axio Observer, 20× objective). Images were first saved on the camera local memory (10.6 GB) and then moved to a desktop computer. These settings allowed for a continuous optical recording of 25 s.
In order to ease the image segmentation step (i.e. the identification of beads in the recorded frames), the electrical velocity V of each event detected in the impedance data streams was used to predict particle entrance time in the optical region, thus reducing the computational cost of image processing. For each particle, the optical lateral position X opt was computed as the distance between the particle centre and the channel z-axis (Fig. 4a), whereas the optical velocity V opt was computed as the product of the distance travelled by the particle across two consecutive frames and the camera frame rate (Fig. 4b). The optical vertical position Y opt (Fig. 4c) was obtained by using a quantitative defocusing approach (Wu et al. 2005;Winer et al. 2014) based on calibration data relating the y-position of a particle to its observed central intensity ).

Particle inertial focusing
The onset of particle focusing can be observed in the laminar regime, when inertial effects start to become relevant. It can be predicted by particle Reynolds number, defined for a rectangular channel as (Di Carlo et al. 2007): where Re is channel Reynolds number, a is particle diameter, U is fluid maximum velocity, D h is cannel hydraulic diameter (defined as 2 WH∕(W + H) ) and μ and are fluid dynamic viscosity and density, respectively.
At Re p ≥ 1 , randomly dispersed particles in a rectangular microchannel tend to order in specific equilibrium positions. They first migrate from the channel bulk towards equilibrium positions near the walls (fast migration), and then, they migrate parallel to channel walls into wall-centred equilibrium positions (slow migration) (Zhou and Papautsky 2013). In a square channel, there are four stable position in correspondence with the middle of each wall, whereas in a rectangular channel the positions close to the short walls become unstable and only two equilibrium position are left (Zhou and Papautsky 2013;Amini et al. 2014).

Information provided by the electrical metrics
R and P Figure 5 shows the density plot for 6-μm-diameter beads pumped at 3 μl/min, with the relative prominence P plotted against the peak-to-peak time ratio R. The particle Reynolds number Re p , equal to 0.11, implies that little inertial focusing is established, so that particles adopt a random distribution in the rectangular channel cross section. Since a one-to-one monotonic decreasing relationship prevails between R and x [resp. P and y], topological information on particle position in the cross section is provided by the (R, P) density plot.
In particular, particles are expected to visit a rectangular region, as experimentally confirmed. Examples of experimental single-particle impedance signals are shown in Fig. 5, along with arrows indicating the position of the particle on the density plot. Particles (i), (ii), (iv) and (v) have trajectories near the corners of the cross section. As an example, particle (i) flows near the lowerleft channel corner, because its signal exhibits high relative prominence P (i.e. low y-coordinate) and high peak-to-peak  Fig. 1a. a Lateral position X opt is determined by measuring the distance between particle centre and channel z-axis. b Velocity V opt is determined by measuring the distance travelled by a particle along the z-axis between two consecutive frames and multiplying it by the camera frame rate (1∕Δt) . c Vertical position Y opt is determined by measuring particle central intensity and referring to a quantitative defocusing calibration set time ratio R (i.e. low x-coordinate). Particle (iii), exhibiting intermediate values of R and P and shortest Δt BDG (i.e. highest velocity), flows along the channel axis.

Inertial focusing detection
The topological information provided by the electrical metrics R and P, along with the electrical velocity V, was used to investigate inertial focusing effects. To this aim, different degrees of focusing were obtained experimentally by varying flow rate and beads diameter, thus varying Re p . In particular, Re p values from 0.11 (6-μm-diameter beads, 3 μl/min) to 3.14 (10-μm-diameter beads, 30 μl/min) were considered (Fig. 6). As shown in the (R, P) plane, at Re p = 0.11 particles were randomly distributed in the cross section (Fig. 6a). By increasing Re p , particles migrated away from the channel centre towards lines parallel to the four channel walls (Fig. 6, first column). The higher the particle Reynolds number Re p , the more pronounced the hydrodynamic focusing. This effect was analysed on the (R, V) plane as well (Fig. 6, second column): at low Re p , the complete parabolic velocity profile could be observed, whereas when Re p was increased, the high-velocities region (relevant to particles travelling close to the channel centre) became progressively less populated and visited by fewer particles. At Re p = 3.14 , the parabolic velocity profile collapsed to the zones corresponding to velocities of the preferential positions parallel to channel walls. The same behaviour could be observed on the (P, V) plane (Fig. 6, third column).
At all values of Re p , the ranges of R and P values visited by the 10-μm-diameter beads were slightly narrower than those of the 6-μm-diameter beads, since a narrower cross section is available to the centres of the former particles, due to their larger size. Moreover, compared to 6-μ m beads, 10-μ m beads showed more pronounced focusing when operated under the same hydrodynamic conditions (i.e. same Re p ). This might be explained considering that for 10-μ m beads the ratio of bead diameter to channel hydraulic diameter, a∕D h , is close to 0.4 and additional steric effects become relevant (Amini et al. 2014).

Calibration procedure used to map the (R, P)-plane onto the (X, Y)-plane
In order to transform the electrical metrics (R and P) into electrical estimates (X and Y) of the Cartesian coordinates of particle centre (x and y), model Eqs. (3) and (5) were used. The relevant model parameters were calibrated using the experiment at Re p = 0.11 (6-μm-diameter beads at 3 μ l/min flow rate), where particles adopt a random distribution in the channel cross section (Fig. 6a-c). In particular, Density plot for 6-μm-diameter beads (flow rate 3 μl/min), with the relative prominence P plotted against the peak-to-peak time ratio R. i-v Are experimental single-particle signals for the data points in the density plot ( I TSV : orange line; I OBQ : green line; I BDG : blue line).
Fitting templates are also shown (grey lines). The histogram of the root-mean-square error of the fit, normalized by peak amplitude, is reported in Supplementary material, Figure S2 6 Density plots of relative prominence P versus peak-to-peak time ratio R (first column); electrical velocity V versus peak-to-peak time ratio R (second column); electrical velocity V versus relative prominence P (third column). Six-μm-diameter beads (odd rows) or 10-μm-diameter beads (even rows) were pumped through the device at flow rates of a-f 3 μl/min, g-l 10 μl/min, and m-r 30 μl/min. The higher the particle Reynolds number Re p , the more pronounced the hydrodynamic focusing. At least 6500 events are reported in each density plot and colour-bar values are normalized with respect to the highest bin count (colour figure online) accounting for finite particle diameter a = 6 μ m, and allowing for a small gap, g ≈ 2 μ m, between particle boundary and channel walls, particle centre coordinates satisfy the following inequalities: where s = (a∕2) + g . Accordingly, parameters b 1 and b 2 appearing in Eq.
(3) were determined by imposing the conditions: where R L and R H are the 1 and 99% percentile values of R, respectively (Fig. 7a, red lines). Their values are reported in Table 1, along with the resulting calibration coefficients b 1 and b 2 . An analogous procedure was used to calibrate parameters c 0 , c 1 and c 2 appearing in Eq. (5). However, an additional condition was needed in this case. Figure 7b shows the density plot of the electrical velocity V versus the relative prominence P. The curve V max (P) that gives the maximum velocity V measured at each relative prominence value P is also plotted (green curve). It was computed using cubic 95% quantile regression instead of maxima search, in order to obtain a robust algorithm against outliers. Because, for any y, the maximum particle velocity in the channel cross section is attained at x = 0 , and a one-to-one correspondence between y and P prevails, the curve V max (P) represents the velocity profile on the line x = 0 , as a function of P. That curve attains its maximum value at a relative prominence value P M which corresponds to y = 0 , since the maximum velocity is attained at channel centre.
Therefore, the parameters c 0 , c 1 , c 2 entering Equation (5) were determined by imposing the conditions: where P L and P H are the 1 and 99% percentile values of P, respectively (Fig. 7a, green lines). The values of P L , P M and P H are reported in Table 2, along with the resulting calibration coefficients c 0 , c 1 and c 2 .

Optical validation of electrical position (X, Y)
and electrical velocity V The electrical estimates of particle lateral position X (Eq. 3), vertical position Y (Eq. 5) and velocity V (Eq. 6) were compared with the optical estimates X opt , Y opt and V opt in Fig. 8 (6-μ m beads, 10 μl/min flow rate). The density plot of the electrical lateral position X against the optical lateral position X opt shows a very good agreement between the two estimates along the whole position range (Fig. 8a, 0.97 correlation coefficient, 2.2 μ m rootmean-square difference, that is 5.6% of channel width). It is pointed out that even the very simple estimate X lin proposed in Eq. (2), based on a geometric argument and not requiring any calibration, provides a reasonably accurate estimate of particle lateral position as shown by the density plot of X lin versus X opt (Supplementary material, Figure S4).
A reasonable agreement is found between the electrical vertical position Y and the optical one Y opt (Fig. 8b, 0.94 correlation coefficient, 1.7 μ m root-mean-square difference, that is 8.1% of channel height). As a matter of fact, the procedure used to estimate the lateral position X opt , which is based on bead segmentation within image frames, is more accurate than the procedure used to estimate the vertical position Y opt , which uses a calibration-based defocusing approach . Moreover, the electrical estimate of the vertical position Y, based on the relative prominence P, may exhibit some blurring due to fitting inaccuracy at low relative prominence values (corresponding to high Y values). It is noticed that the electrical method used for the lateral

Fig. 7
Calibration experiment (6-μ m beads, 3 μl/min flow rate, Re p = 0.11 ). a Density plot of relative prominence P versus peak-topeak time ratio R. The 1 and 99% percentile values of R (red lines) and P (green lines) are shown. The corresponding values are labelled as R L , R H (Table 1), and P L , P H (Table 2). See also Supplementary material, Figure S3. b Density plot of electrical velocity V versus relative prominence P. The cubic 95% quantile regression is shown in green. The relative prominence value P M corresponding to the maximum velocity V max = 0.12 m/s is visualized (Table 2). Colour-bar values are normalized with respect to the highest bin count (colour figure online)  Table 2 Values of 1 and 99% percentile of P (Fig. 7a) and value of P at maximum velocity (Fig. 7b). Resulting calibration parameters (Eq. 5) position determination could be replicated in a vertical configuration using five pairs of facing electrodes (Spencer et al. 2016), thus providing an even more accurate estimate of particle vertical position, at the expense of a more complex microfabrication process. The density plot of the electrical velocity V against the optical velocity V opt is shown in Fig. 8c. An excellent agreement is found between the two measurements (0.97 correlation coefficient, 0.01 m/s root-mean-square difference, that is 2.3% of maximum particle velocity). The velocity range measured optically and electrically is (0.20-0.39) m/s, with average particle velocity of 0.30 m/s. By comparison, the velocity distribution in steady state, hydrodynamically fully developed, laminar flow for Newtonian fluids in rectangular channels (Spiga and Morini 1994) has maximum 0.39 m/s and average 0.19 m/s for the present channel cross section and flow rate. In fact, particles tend to avoid trajectories very close to the channel walls, so that their average velocity is greater than average fluid velocity.

Quantitative particle positioning within the channel cross section
The R ↦ X and P ↦ Y mappings (namely Eqs. 3 and 5 with the coefficient values reported in Tables 1 and 2) were applied to the data sets considered in Fig. 6. The resulting (X, Y), (X, V) and (Y, V) density plots are shown in Fig. 9. For each experimental condition, the percentage of particles falling in each bin of a 3 × 3 regular grid in the (X, Y) plane (Fig. 9, first column) is also indicated (bin size W∕3 × H∕3). At Re p = 0.11 , when particles are randomly distributed, 21% of the particles visited the central bin and lower percentages were found in the other bins (especially the corner bins), due to a minimum distance of particle centres from channel walls imposed by finite particle size. By increasing Re p , the particles moved away from the central bin towards peripheral bins. At Re p = 3.14 , only 1% of the particles visited the central bin, and the highest particle concentration was found in the bin corresponding to the centre of the long walls, in agreement with the literature (Zhou and Papautsky 2013;Amini et al. 2014) (see Supplementary material, Figure S5).
The evolution of particle distribution in the channel cross section as Re p increases is reflected in changes of the velocity profiles observed in the (X, V) and (Y, V) planes (Fig. 9, second and third column, respectively).

Conclusion
In this work, we reported the first high-throughput electrical system for the measurement of cross-sectional position and velocity of single particles flowing through a microchannel. Lateral position was computed measuring the ratio of oblique and transverse transit times, whereas vertical position and velocity were extracted from the relative prominence and transit time of a bipolar double-Gaussian signal measured with a suited wiring scheme. Transit-time ratios and relative prominence values supplied valuable topological information on particle position in the microchannel cross section. Moreover, they could be transformed to quantitative estimates of lateral and vertical positions using a simple calibration procedure. Optical validation confirmed the soundness of the proposed approach. As an applicative example, the effects of particle size and flow rate on inertial particle focusing were investigated. Fig. 8 Comparison between electrical estimates X, Y, V and optical estimates X opt , Y opt , V opt (6-μ m beads, 10 μl/min flow rate). Data are shown as density plots (2500 events), and the bisector line is dashed in red. Colour-bar values are normalized with respect to the highest bin count. a Electrical lateral position X (Eq. 3) versus optical lateral position X opt (0.97 correlation coefficient, 2.2 μ m root-mean-square difference, that is 5.6% of channel width). b Electrical vertical position Y (Eq. 5) versus optical vertical position Y opt (0.94 correlation coefficient, 1.7 μ m root-mean-square difference, that is 8.1% of channel height). c Electrical velocity V (Eq. 6) versus optical velocity V opt (0.97 correlation coefficient, 0.01 m/s root-mean-square difference, that is 2.3% of maximum particle velocity) (colour figure online) Effect of particle Reynolds number Re p on particle focusing. Measurements of 6-μ m (odd rows) or 10-μm-diameter (even rows) beads, pumped through the device at 3 μl/min (a-f), 10 μl/min (g-l) or 30 μl/min (m-r) flow rate, are shown. For each experimental condition, density plots of electrical position Y versus electrical position X (first column), electrical velocity V versus electrical position X (second column), electrical velocity V versus electrical position Y (third column) are reported. In the (X, Y) density plot (first column), the percentage of particles falling in each bin of a 3 × 3 grid (white dotted lines) is also indicated. The higher the particle Reynolds number Re p , the more pronounced the hydrodynamic focusing. At least 6500 events are reported in each density plot and colour-bar values are normalized with respect to the highest bin count (colour figure online) The method has been demonstrated with reference to spherical beads. It applies as well to spherical cells, whose behaviour at low frequency (i.e. below the range of the -dispersion) is similar to that of an insulating particle (e.g. Schade-Kampmann et al. (2008), Haandbaek et al. (2014)). The method is also promising for nonspherical particles or cells because, although the latter yield impedance signals depending on their orientation (Jones 1995;Fernandez et al. 2017), the transit times and the relative prominence are not significantly affected by particle orientation.
This fully electrical label-free approach to particle monitoring, implemented with a simple coplanar chip layout and based on current signals instead of massive image data, represents an effective tool for investigating the behaviour of flowing particles at the microscale. Moreover, the device can be simultaneously used to perform single-cell impedance spectroscopy, thus achieving a multi-parametric characterization embracing size, membrane properties and intracellular conductivity in addition to particle position and velocity.