A density-based threshold model for evaluating the separation of particles in heterogeneous mixtures with curvilinear microfluidic channels

Particle separation techniques play an important role in biomedical research. Inertial focusing based microfluidics using nonlinear channels is one of the promising label-free technologies for biological applications. The particle separation is achieved as a result of the combination of inertial lift force (FL) and Dean drag force (FD). Although the mathematical expressions of FL and FD have been well derived in prior studies, they are still complicated, which limits their popularity in practice. Recent studies modified these expressions through experiments and proposed a threshold model, which assumes that only particles larger than the threshold will be well focused. Although this threshold model has been used in recent studies, two varying versions of the threshold model (TM1 and TM2) prevents standardisation in practice. In addition, both models were developed with regular low-density particles and may not be applicable to samples with higher density or samples with irregular shapes. Here, we evaluated the threshold models with samples of different densities. Based on these evaluations, we derived a modified model (TM4), which additionally considers the factor of particle density to improve the accuracy of existing models. Our results demonstrated that TM4 could more reliably predict the sorting efficiency of samples within a wider density range.

Scientific Reports | (2020) 10:18984 | https://doi.org/10.1038/s41598-020-75878-w www.nature.com/scientificreports/ where ρ is the fluid density, U is the maximum fluid velocity, μ is the fluid viscosity, and D is the hydraulic diameter expressed as: where h and w are the height and width (for rectangular cross-section), respectively. In inertial microfluidics, the resultant force due to shear-gradient and wall-induced forces acting on particles in the opposite direction is called F L, which is given by: where C L is lift co-efficient, which is a function of particle position, and a is the particle diameter. Apart from inertial lift force, particles flowing in the curved channel will experience a F D resulting from 2 counter-rotating vortices (Dean vortices), which can be quantified by Dean number (De) as: where R is the radius of curvature of the channel path. The Dean vortices induced at the channel cross-section will cause particles to move back and forth across channel width. This lateral velocity for particles, Dean velocity (U D ), is estimated by: With Dean velocity (U D ), the F D acting on the particles can be expressed as: The focusing effect for given particle dimensions is determined by the ratio of F L and F D (R f = F L /F D ), which varies with particle diameter (R f ∝ a 3 ). The largest particles will be focused on the channel position closest to the inner wall as they experience the larger F L compared to other particles at that position in the system. A threshold model was derived based on the modification of previously established Eqs. (3), (4), (5), and (6), among which only particles larger than the threshold will experience significant F L, thereby balancing F D to facilitate particle focusing. The threshold of particle focusing is proportional to the size of the microchannel, which can satisfy the following inequality, as stated: where x = D( hydraulic diameter) for TM1 and x = h (channel height) for TM2.
Particle recovery rate (PRR) is often used to quantify the particle sorting performance in microfluidics 23 . In this study, the PRR of the device used ( Fig. 1) is defined as: Evaluation of TM1 and TM2 with regular low-density particles. Given that TM1 and TM2 were developed with regular low-density particles ( ∼ 1.05 g/ml), both threshold models were initially analysed by COMSOL, and then polystyrene (PS) microsphere experiments were performed. According to Eq. (7), the estimated thresholds of TM1 and TM2 are 21 μm and 15 μm, respectively. Therefore, 21 μm and 15 μm regular PS particles were first studied. In the COMSOL simulation, it was found that 21 μm particles (threshold of TM1) were collected at outlet 4 ( Fig. 2a), and 15 μm particles (threshold of TM2) were collected at outlets 3 and 4, respectively (Fig. 2b). Using the estimated values provided by COMSOL, we proceeded to conduct actual experiments to validate this hypothesis further.
For the experiments with 15 μm and 21 μm PS microspheres, more than 90% (91.6 ± 7.3%) of 21 μm particles were recovered in outlet 4 (Fig. 2c), and the rest of the particles was found at outlets 3 and 5 (5.2 ± 3% and 3.2 ± 4.4%, respectively). The particles found in outlets 3 and 5 could be a result of unfocused streamlines in the first few seconds of flow (Fig. S1). For 15 μm particles, the PRRs at outlets 3 and 4 were 52.4 ± 4.6% and 47.6 ± 4.7%, respectively (Fig. 2d). The highest PRR of 21 μm microspheres indicated that the threshold estimated by TM1 for PS microspheres was more accurate than TM2 for devices with similar spiral radius and aspect ratio parameters.
Evaluating the overestimation of TM1 threshold. The previous results cannot exclude the possibility of overestimating the threshold. For example, even if the threshold in the system was 21 μm (TM1), particles larger than 21 μm (such as 25 μm) would still be well focused, but 25 μm was definitely not the precise threshold value. To address this issue, we studied the PRRs of microspheres in the range of 21 to 15 μm. If the recovery     (Table S2). All 18 μm particles (T t1 = 18 μm) were well focused on outlet 4 ( Fig. 3a), which was similar to the PRR of 21 μm particles (Fig. 2a). Hence, the threshold value of 16.5 μm (T t2 = (18 + 15)/2 = 16.5 μm) (test 2) was used for the next test. Then, all 16.5 μm particles were well focused on outlet 4 ( Fig. 3b), which was similar to the PRR of 18 μm particles (T t1 ). Hence the threshold value of 15.7 μm (T t3 = (16.5 + 15)/2 ∼ 15.7 μm) was used for the next test (test 3). The simulation of 15.7 μm particles showed that only 66.7% particles (n = 3) were concentrated on outlet 4 (Fig. 3c), significantly lower than the PRR of 16.5 μm particles. As a result, the target threshold range was likely to be between 16.5 μm and 15.7 μm (15.7 μm < T < 16.5 μm), and the theoretical threshold value for test 4 was defined as 16.1 μm. Since the precise threshold range was between 16.1 μm and 16.5 μm (16.1 μm < T < 16.5 μm), the scale of overestimation in practice was negligible. Hence, the revised upper limit of 16.5 μm was used as a proposed threshold for regular low-density samples. For the corresponding threshold model (TM3) modified from TM1 and TM2, the x in Eq. (7) was taken to (D + 3h)/4.
The outcomes of COMSOL simulation were validated by actual experiments. Specifically, we evaluated TM3 with 16.5 μm PS microspheres (Fig. 3d). The highest PRR was found at outlet 3 (77.6 ± 6.2%), while the rest of the microspheres was found at outlet 4 (22.4% ± 6.2%). Although the hPRR of 16.5 μm microspheres at outlet 3 were higher than that of 15 μm microspheres (52.4 ± 4.6%), it was lower than that of 21 μm microspheres (91.6 ± 7.3%) (Fig. 3e). Thus, TM1 remains the most accurate model for regular low-density particles in this study.
Evaluation of all models with irregular low-density particles. One of the main aims of redefining the model used to estimate the separation threshold is to expand the application range of inertial focusing techniques for samples of different densities, shapes, and sizes. Here, we investigated the accuracy of all three threshold models for irregularly shaped samples using polystyrene (iPS) particles of various shapes and sizes (< 100 μm). Through actual experiments, we obtained the recovery rates of 15 μm, 16.5 μm 21 μm iPS samples (± 2 μm) in outlets corresponding to the thresholds of TM1, TM2, and TM3. The hPRRs of 21 μm, 16.5 μm and 15 μm iPS particles were 41 ± 4.8% (outlet 3), 31.2 ± 10% (outlet 4) and 27.6 ± 4.9% (outlet 3) (Fig. 4). Then, we evaluated whether TM1 will underestimate the threshold by obtaining the PRR of the larger iPS particles Evaluation of all models with irregular high-density particles. The focus position of particles within the channels depends on the particle density 12 . However, since inertial microfluidics research is mainly designed for biomedical applications, low-density particles such as polystyrene microbeads ( ∼ 1.05 g/ml) were often used because their densities are similar to that of biological samples. In recent years, scientists are now interested in exploring the use of inertial microfluidics for higher density targets. Hence we used irregular talcum (iTalc) particles ( ∼ 2.75 g/ml) to investigate how particle density could affect the accuracy of the existing threshold models. It could be seen from the results that the hPRRs of 21 μm, 15 μm and 16.5 μm iTalc particles at outlet 3 were 38 ± 11.5%, 31.3 ± 9.1%, and 32.2 ± 9.9%, respectively (Fig. 5). These low PRRs of respective diameters were probably due to the high density of talcum. According to Eqs. (3) and (6), F L and F D are independent of the change in particle density, and thus an increase of particle density might not influence the balance of F L and F D . However, in inertial microfluidics, apart from the F L and F D , the particles migrating in the curved channel were also subjected to the centrifugal force F c which is given by: where ρ p and ρ f are the particle and fluid densities, respectively, V is the tangential velocity of the particles, and r is the radius of the particle trajectory 10 . In most studies, F c is usually ignored for proximate densities of particle and fluid medium ( ρ p − ρ f ) 12 . For regular and non-regular PS samples, since the densities of PS ( ∼ 1.05 g/ml) and water ( ∼ 1 g/ml) are very similar, F c acting on microbeads is negligible, and thus particle focusing depends on the balance of F L and F D (F L = F D ). However, for higher density particles such as talcum, the F c acting on the particles would be higher (keeping other parameters unchanged), and thus F c became non-negligible. When the particles were subjected to Dean vortices in the inner cross-section, F c acted in the direction opposite to F L , and therefore a larger F L was required to balance the additional F c and to focus the particles of higher density (F L = F D + F c ). Hence, the deviation of PRRs values from TM1 estimation strongly suggested the need for a more precise model for high-density targets.
Modification of TM1 for irregular high-density particles. Since larger particles will experience a larger F L in the system, the threshold for iTalc particles was expected to be larger compared to that of iPS particles. In order to derive a more precise model for samples with varying densities, we obtained the PRR of the larger iTalc particles (26 ± 2 μm) at each outlet. The results showed that the highest PRR of 26 ± 2 μm iTalc www.nature.com/scientificreports/ particles (50.2 ± 12.6%) found at outlet 3 was significantly higher than that of 21 μm (38 ± 11.5%) (Fig. 5d). This result was in line with our expectation that the larger F L induced by larger particles would balance the additional F c . Then, by analyzing the PRR of iTalc particles larger than 26 ± 2 μm (28-50 μm), we assessed whether the threshold was saturated at 26 ± 2 μm. The hPRR for larger iTalc particles at outlet 3 was found to be 48.6 ± 15.3%, which was very similar to the PRR of the smaller iTalc particles (26 ± 2 μm) (Fig. 5e). Therefore, the threshold for high-density particles (iTalc) was probably saturated around 26 μm. Similar to the result of iPS particles, a slight increase in size did not increase the highest PRR, despite the larger F L experienced by larger particles. This could be due to the presence of other opposite forces, F c and F D, which would also be increased in the system (F c ∝ a 3 and F D ∝ a ). Therefore, there was a saturation point at which the focusing effect did not increase as the particle dimension increased. To allow estimation of PRRs for samples with varying densities, we thus derived a model (TM4) as: where c was a density coefficient greater than or equal to unity (c ≥ 1). The value of c increased with the density of sample, for example, for polystyrene particles ( c = 1; ∼ 1.05 g/ml) and for talcum particles ( c = 1.215; ∼ 2.75 g/ml). Since the c value is equal to 1 for PS, and 1.215 for Talc, substituting above c values to TM4, the thresholds for iPS and iTalc particles were calculated to be 21 μm and 26 μm, corresponding to the highest PRRs obtained previously (Figs. 4 and 5).
Applicability of TM4 for irregular particles of different densities. In order to study whether TM4 could be applied to particles of different densities, we hypothesized that the density coefficient c of TM4 would increase linearly with particle density, from 1.05 g/ml to 2.75 g/ml (Fig. 6a), whereby the proposed linear correlation can be modelled by: where c is the density coefficient, and ρ is the particle density, respectively. We aimed to verify this linear correlation using irregular polyethylene terephthalate (iPET) ( ∼ 1.38 g/ml) and polyamide (iPA) ( ∼ 1.15 g/ml). According to Eq. (11), the density coefficient (c) of iPET was around 1.04. Based on TM4, the threshold for iPET was approximately 22.5 μm. We then obtained the PRRs of iPET particles of various sizes (22.5 ± 2 μm, 26.5 ± 2 μm, 20.5 ± 2 μm, and 28.5-50 μm). We found that the highest recovery rates of 22.5 μm, 20.5 μm, and 26.5 μm iPET particles were around 46.9 ± 2.3%, 42.1 ± 12%, 43.9 ± 4.3% and 46.2 ± 5.1% respectively (Fig. 6b-e). As the highest PRR of 22.5 μm particles was higher than that of 20.5 μm particles and close to 26.5 μm and 28.5-50 μm iPA particles, 22.5 μm was determined as the threshold for iPET particles. Then, iPA particles were  (Fig. 7a-d) showed that the threshold for iPA was around 21.5 μm (60% hPRR at outlet 3) because the hPRR of 21.5 μm iPA was higher than that of 17.5 ± 2 μm and similar to that of 25.5 ± 2 μm. Therefore, the linearity was validated (Fig. 7e).
Overall, we derived an improved model (TM4) to estimate the separation efficiency of samples with varying densities by incorporating the factor of particle density. Our studies demonstrated that TM1 was only suitable for the low-density particle (1.05 g/ml), while TM4 was applicable for samples with a broader range of densities, including high-density samples (2.75 g/ml). The linear correlation between c and ρ makes TM4 more applicable to samples of different densities. For example, since the value of c is 1 in TM4 for PS particles, TM4 will be exactly equivalent to TM1, which was the most accurate for PS particles, and thus TM4 retained the same accuracy as compared to TM1 for low-density samples. For intermediate-density and high-density particles (e.g., PET and Talc), since the centrifugal force acting on them was non-negligible, the parameter was considered in the derivation of TM4 by introducing a density coefficient c, which was customized for particles of different densities. Specifically, the threshold of suitable density range for TM4 was between 1.05 to 2.75 g/ml. We hope these studies can also provide a systematic approach for future efforts to derive an optimal prediction model for a wider range of biological samples.

Discussion
Some parameters essential to channel design, such as C L and U D in Eq. (3) and (5) can only be estimated through experimental data 12 . These equations can be modified by utilizing the size-dependent nature of R f (R f = F L /F D and R f ∝ a 3 ). However, there is a lack of standardisation of the threshold model used, leading to a significant discrepancy in threshold prediction. In addition, current models were developed using low-density samples ( ∼ 1.05 g/ml) and may not be applicable to samples of higher density. We believe that the applications of inertial microfluidics will be extended to the use of high-density particles in the foreseeable future. One potential application is to facilitate the monitoring of heavy metals in drinking water. These high-density metals (such as Ag, Cu, Ni, and Cr) are carcinogens present in very low concentrations. Nonetheless, current point-of-care technologies for monitoring metals are subjected to high limits of detection (LOD), of which the LOD values are higher than the maximum allowable limit suggested by WHO 24 . Therefore, an improved threshold model Figure 6. The proposed linear correlation between density coefficient, particle density, and the recovery rate of irregular polyethylene terephthalate (iPET) of various diameters. (a) The graph illustrates that c linearly increased with ρ within the range of 1.05 g/ml to 2.75 g/ml. The PRRs of (b) 22.5 ± 2 μm, (c) 20.5 ± 2 μm, (d) 26.5 ± 2 μm, and (e) 28.5-50 μm were indicated respectively. Since the hPRR of 22.5 μm was similar to 26.5 μm but higher than that of 20.5 μm, 22.5 μm was determined as the threshold for iPET. Data were shown as mean ± SD of three independent experiments. *p < 0.01, **p < 0.001, ***p < 0.0001. Scientific Reports | (2020) 10:18984 | https://doi.org/10.1038/s41598-020-75878-w www.nature.com/scientificreports/ incorporating sample density for point-of-care inertial microfluidic biosensors can be useful for the detection or removal of toxic metals from drinking water. Overestimation of the threshold, even as small as a few microns, will also significantly affect the design of a microfluidic device and prevent the optimal separation of the subpopulations of cells. For example, the diameter of white blood cells ( ∼ 9.5 μm) is similar to that of red blood cells ( ∼ 7 μm). Here, we evaluated the accuracy of the existing models (TM1 and TM2) and derived an improved model (TM4) to estimate the separation efficiency of samples with varying densities by incorporating the factor of particle density. TM4 provides a more precise estimation of the separation efficiency of samples with varied densities and will aid researchers in their design of microchannels for particle separation based on the inertial focusing effect.
Similar to TM1 and TM2, it should be highlighted that TM4 is only applicable to channels with a rectangular cross-section. Interestingly, the previous threshold models did not consider geometric considerations, and both TM1 and TM2 were derived based on rectangular cross-section channels. Since the focus of this manuscript is to expand on existing threshold models (TM1 and TM2) and incorporate the component of sample density for inertial microfluidics, the standard rectangular cross-section design was used for our simulation and experimentations. Our preliminary studies (Fig. S2) suggested that all current threshold models (including TM4) were indeed not applicable for microfluidic channels with varied geometric designs (e.g., trapezoidal). Future efforts should be directed towards the modification of existing models for different geometric shapes.
Other factors can also influence the sorting efficiency of particles. Previous studies have also emphasized that unless the asymmetry is large, the shape of the particle has only a small effect on inertial focusing, although the relevant mechanism was not clarified in the literature 25 . In order to show whether there were some PS powders with great asymmetry, a scatter diagram of the circularity distribution of iPS powders was generated (Fig. S3). Since the scattering points were uniformly distributed from 0.1 to 1, the shape of the particles was highly diverse. This implied that a certain portion of particles could have a higher degree of asymmetry. As previous research on the correlation between the particle shape and the focus position was carried out with straight channels, further studies will be required to consider how Dean resistance is affected by the shape of the particles.
The experimental results of the microbeads were different from that of the simulated runs. For example, experimental results showed that outlet 4 had a higher PRR of 21 μm PS microbeads, which could be due to the stronger F L experienced. On the contrary, the 16.5 μm microbeads were less focused in outlet 3, possibly due to a moderate F L acting on the microbeads 26 . The shift of the equilibrium position from outlet 4 to 3 in experimental results could also be due to differential positioning depending on the nature of F L . Particles circulating along with Dean vortices at the inner section would experience larger F L at the proximate center location than at proximate inner channel wall 7 . Compared with 21 μm microbeads, 15 μm and 16.5 μm microbeads could experience Figure 7. Recovery rate of irregular polyamide (iPA) particles of (a) 17.5 ± 2 μm, (b) 21.5 ± 2 μm, (c) 25.5 ± 2 μm, (d) 27.5-50 μm diameters. As the highest recovery rate of 21.5 ± 2 μm iPA (59.8%) was higher than that of 17.5 ± 2 μm (52.1%), and very close to that of 25.5 ± 2 μm (55.7%), 21.5 μm was determined as the threshold for iPA. Data were shown as mean ± SD of three independent experiments. *p < 0.01, **p < 0.001, ***p < 0.0001. (e) The linearity of density coefficient and particle density. Points A, B, C, and D corresponded to the experimental data obtained with iPS, iPA, iPET, and iTalc particles, respectively.

Scientific Reports
| (2020) 10:18984 | https://doi.org/10.1038/s41598-020-75878-w www.nature.com/scientificreports/ sufficient F L to balance F D when they were located in a lateral position closer to the center of the channel. In this study, COMSOL was only a means to provide a rough estimation, as the accuracy of COMSOL was also influenced by many factors, such as the number of mesh elements. Increasing the number of elements generally led to higher accuracy, but the demand for computation power was also higher. In this study, the number of mesh elements was standardised to around 37,000. Despite these limitations, we demonstrated that the outcomes of the numerical simulation were close to our actual experiments of 15 μm and 21 μm microspheres. Biological samples are often heterogeneous, especially when clinical samples are involved. Prior studies on the detection of mesenchymal cancer cells from bladder wash urine and blood samples had demonstrated the fluctuations in target cell retrieval in clinical samples 23 . Here, we evaluated the performance of different models using particles of different size, shape, and density. Our result demonstrated that TM1 was a more accurate model for low-density particles. However, PPRs of higher density samples deviated from expected values. We thus modified TM1 by considering the effect of the additional centrifugal force F c acting on high-density particles in the curvilinear channel. The derived model (TM4) comprises of a density coefficient c such that the threshold obtained fits experimental values. The robustness of TM4 for samples of varying densities (1.05 g/ml ≤ ρ ≤ 2.75 g/ml) was confirmed with intermediate-density particles (PET and PA). In addition, TM4 was also as accurate as TM1 for low-density particles (Table S3).
Future analysis should cover other parameters that could influence the threshold, such as changes in channel aspect ratio, channel radius (R), and cross-section dimensions (same aspect ratio), which were reported to affect inertial focusing 27 . For example, apart from shape and density, particle deformability is one of the factors affecting the inertial effect owing to an additional life force experienced by non-rigid particles 28 . Similar to other inertial studies using polymeric particles (usually PS), all particles were assumed to be rigid because of their relatively high elasticity ( ∼ 2 − 4GPa ). However, the mechanical characterization of cells is similar to a drop of liquid enclosed by a membrane exhibiting low elasticity 29 . For example, the elastic modulus of mammary cells (MCF10A) and breast cancer cells (MCF7) widely used for studying breast cancer were around 440 Pa and 300 Pa, respectively 30 , while the stiffness of healthy leukocytes was even lower (10-60 Pa) 31 . Based on differences in deformability, isolation of more deformable MCF7 and healthy leukocytes from less deformable MCF10A and immature leukocytes had been reported in previous studies despite their similarity in sizes 31,32 .
Overall, we demonstrated that TM1 remained accurate for regular and non-regular low-density particles (e.g., PS). However, a modified model TM4, which incorporated the density coefficient c (TM4), was more suitable for samples with varied densities, especially those with high-density particles (e.g., Talc). The linear correlation of density coefficient c and particle density ρ was also validated. These findings are highly applicable to applications where the density of targeted particles is high or heterogeneous.

Methods
Device feature and fabrication. The spiral device with a rectangular cross-section (the main device in this study) consisted of a 10-loop curvilinear channel with 300 μm spacing between loop interval. The width and height of the channel were 500 μm and 220 μm, respectively. The device had 1 inlet and 5 outlets. The end of the 500 μm wide channel was evenly split into five outlets, each of which had 100 μm width. Spiral devices with trapezoidal cross-section consisted of an 8-loop curvilinear channel with 300 μm spacing between loop interval. The width of the channel was 600 μm, while the inner and outer heights were 80 μm and 130 μm, respectively. The device had 1 inlet and 2 outlets. The 600 μm wide channel was evenly split into two outlets.
Microdevices were fabricated using standard soft lithography from an aluminum master mold. Briefly, polydimethylsiloxane (PDMS) was prepared by mixing a 10:1 ratio of base and curing agent (Dow, Germany). The prepared PDMS was degassed in degasser and poured onto the master mold. The mold filled with prepared PDMS was baked at 70 °C oven for 2 h. After baking, the solidified PDMS was peeled off from the mold carefully. The locations for inlet and outlets were punched using a 1.5 mm diameter puncher (Integra, USA). The punched PDMS device was bonded to a flat solidified PDMS layer to complete the channels in a plasma machine. After plasma treatment for 3 min, the bonded PDMS device was placed in a 70 °C oven for 30 min to strengthen the bonding. Finally, 15 cm and 6 cm long sections of flexible plastic tubing (Tygon, USA) were inserted into the device inlet and outlets, respectively.
COMSOL flow simulation. To evaluate the particle focusing of different sizes for respective threshold models. COMSOL Multiphysics software was used to simulate the inertial focusing on our existing microdevice. The three-dimensional (3D) image of the microdevice was created by AutoCAD in STL format (Fig. 1a). All dimensions were the same, except for the shorter length of the branch at all outlets (to save computation time). The properties of water were applied to systematic fluid for simulation (ρ = 1000 kgm −3 ; µ = 10 −3 kgm −1 s −1 ). The laminar flow (spf) module with incompressible flow property was used in the system. F low velocity was set at the inlet, while the fixed pressure boundary condition was set at all outlets. Using these perimeters, a flow magnitude profile for the system was simulated (Fig. 1b).
Based on the simulation of the flow profile, particle tracing for the fluid flow (fpt) module was chosen to analyse the motions of particle trajectories. The drag force, wall-induced life force, and shear gradient life force parameters were applied to the system. Interactions between particles were ignored as the diluted sample was often used in the inertial microfluidic application. Solid particles ( 1.05 g/ml) was used to model spherical polystyrene (PS) microbeads. Each set of samples consisted of 3 particles of the same size evenly distributed at the inlet to initiate simulation (Fig. 1a). The simulation time for each set of particles was 2 s with 0.01 s time step for tracing particle trajectories. The 'freeze' condition was applied at the end of outlets to ensure all particles within the boundaries. The number of particles freezing at the end of boundaries was counted.
Scientific Reports | (2020) 10:18984 | https://doi.org/10.1038/s41598-020-75878-w www.nature.com/scientificreports/ Flow rate optimization. For the inertial microfluidic system, the Reynolds number (Re) is required to be within the range of 1 to 100 (~ 1 < Re < ~ 100), which corresponds to the volumetric flow rates approximately from 0.021 ml/min to 2.1 ml/min. Since a higher flow rate was preferable, different flow rates (1.5 ml/min, 1.7 ml/min, 1.9 ml/min, and 2.1 ml/min) were used to process 21 μm PS regular microbeads in the system. 21 μm spherical microbeads were used as the particle size was above the thresholds of both TM1 and TM2, which indicated that 21 μm PS microbeads should be well-focused in the device. We confirmed that based on the highest PRR (91.7% PRR) obtained under this parameter, the optimal flow rate is 1.7 ml/min (Fig. S4). Although the PRR derived from the flow rate of 1.9 ml/min suggested that this parameter could also be used, it was close to the limit of Reynold's number applicable to inertial focusing based microfluidics, and thus only 1.7 ml/min was used for subsequent experiments.
Preparation of irregular particles. iTalc (baby powder) was commercially available (Johnson & Johnson China Ltd, China). For iPS, iPET, and iPA powders, they were prepared using a Retsch CryoMill cryogenic grinder (Haan, Germany). After a 7-min precooling stage with five shakes per second, the grinding stage was used for 1.5 min with 25 shakes per second in the presence of liquid nitrogen. The crushed microplastic powders (< 100 μm) were then collected through a 100 μm filter. The powders were mixed with DI water for experimentation in a biological safety cabinet to minimize potential contamination.
Particle sorting experiments. 1 μl microbead suspension was mixed well with 1 ml deionized water and transferred to 10 cc syringe. Another 10 cc syringe was prepared by injecting 2 ml deionized water and was loaded to a syringe pump. The inlet and outlets of the microdevice were connected to the syringe and five centrifugal tubes using plastic tubing. The device and centrifugal tubes were arranged accordingly (Fig. 1c). The syringe pump was run at 1.7 ml/min for 10 s. The loaded syringe was replaced by another syringe with microbead suspension. The syringe pump (New Era Pump Systems Inc., USA) was run at 1.7 ml/min for 20 s. The sample collected at centrifugal tubes was then transferred to a 24-well plate (SPL Life Science, Korea). Each well with microbead suspension was mixed well with an additional 100 μl deionized water. The 24-well plate was placed at rest for 30 min to settle down the particle suspension. The number of regular particles in each well was counted under an optical microscope (Nikon, Japan) (Fig. S5). All experiments were carried out in triplicates.
For irregular particles, all images captured by the optical microscope (Fig. S6) were analysed by ImageJ. Briefly, images were converted to 8 bits with suitable threshold levels to customize the particle counting function. The adjusted images were carefully compared to their original images to ensure that the threshold selected was optimal. Their circular diameters were measured with ImageJ and Excel. Briefly, we first obtained the boundary area (A) of each particle by ImageJ. Then, the circular diameter (d) was calculated by d = √ (4*A/π) using Excel. Due to the high variation of sample sizes, the upper and lower ranges of two microns (± 2 μm) were included for the measurement of the specific dimension.