3.1. Electrode Dimension Optimization with FEM
The optimal electrode layout and dimensions were determent using FEM simulations shown in
Figure 3 based on the constrictions discussed in
Section 2.1.
Three-dimensional finite element modeling in COMSOL was used to predict the magnitude and uniformity of the electric field between the electrodes when varying their width. This was done by evaluating the normalized current density through a cut plane in the narrow area between the electrodes as indicated in red on
Figure 3A. A 3 V potential at 231 kHz was applied to the central excitation electrode while the two outer measurement electrodes were defined as ground (0 V). The medium between the electrodes was defined as water with a conductivity of 1.6 S/m and a relative permittivity of 80.
Figure 3.
FEM simulation of the current density. (A) 3D sketch of the channel with wide electrodes. The current density was evaluated through the detection cross-section (red surface) between the electrodes using FEM. (B) Normalized current density through the detection cross-section as a function of electrode length the optimized chip. (C) Contour plot of the current density in the detection cross-section for the conventional electrode layout. (D) Contour plot of the current density in the detection cross-section for the improved electrode layout. (E) Line plot of the x-component of the current density for z = 2.5 μm as a function of the y-coordinate across the detection cross section for varying electrode widths. (F) Line plot of the x-component of the current density for y = 0 μm as a function of the z-coordinate across the detection cross-section for varying electrode widths.
Figure 3.
FEM simulation of the current density. (A) 3D sketch of the channel with wide electrodes. The current density was evaluated through the detection cross-section (red surface) between the electrodes using FEM. (B) Normalized current density through the detection cross-section as a function of electrode length the optimized chip. (C) Contour plot of the current density in the detection cross-section for the conventional electrode layout. (D) Contour plot of the current density in the detection cross-section for the improved electrode layout. (E) Line plot of the x-component of the current density for z = 2.5 μm as a function of the y-coordinate across the detection cross section for varying electrode widths. (F) Line plot of the x-component of the current density for y = 0 μm as a function of the z-coordinate across the detection cross-section for varying electrode widths.
The current density through the cross-sectional cut-plane was evaluated by surface integration for channels with electrode widths ranging from 10 μm (flush with the channel wall) to 50 μm (
Figure 3B). The simulation predicts that the current between the electrodes will saturate when the electrode width is in the order of 20 μm, indicating that the design of the chip will not gain any more sensitivity from increasing the width beyond this point. While a larger current is expected to result in better sensor sensitivity it is not necessarily beneficial to extent the electrode width beyond 20 μm. The key disadvantage of widening the electrodes relates to the fluid dynamic properties of the channel. A significant enlargement of the electrode area can give rise to problems with turbulent flow around the edges of the enlargement; which is highly undesirable in a system relying on detection of transitions. Furthermore, issues can arise with the flow in the enlarged channel regions where particle trapping can occur due to areas with dead flow. However, particle trapping was not observed during experiments with the channel layout proposed in this paper.
Figure 3C,D show the contour plots of the electric current density through the cross-sectional detection plane for the original design with an electrode width of 10 μm and the improved design with an electrode width of 20 μm. It is seen that the current density varies more for the channel with the improved layout compared to the channel with the original layout; however the overall current density is 68% higher in the improved channel. The larger current density of the optimized chip should allow for higher signal and thus better SNR.
The current density across the detection cross-section in the
y- and
z-directions for different electrode widths is plotted in
Figure 3E,F. An electrode width of 10 μm corresponds to a conventional, straight channel. The variation of current density in the
z-direction is not influenced much by the enlargement of the electrode width, while in the y-direction there is a significant increase near the channel wall when the electrode is widened. As the flow rate used in this these experiment are sufficiently low, we do not expect that the particles will be influenced by inertial forces from the liquid, which means that the
z- and
y-position of the particle during the transition will be random. Some will transition near the edge and some in the center of the channel. It is therefore expected that while the improved electrode layout will provide an overall higher SNR it will also lead to a larger spread in measured current.
3.2. Chip Characterization
Figure 4A,B shows the impedance characterization of a single set of electrodes from 200 kHz to 10 MHz, of both the conventional and improved chips. The shape of the curves follows the typical trend for a fluidic channel [
3].
Figure 4.
Frequency sweep of chips used for the experimental work. (A) Impedance magnitude as a function of frequency for the conventional and improved chips. (B) Phase shift/angle as a function of frequency for the conventional and improved chips.
Figure 4.
Frequency sweep of chips used for the experimental work. (A) Impedance magnitude as a function of frequency for the conventional and improved chips. (B) Phase shift/angle as a function of frequency for the conventional and improved chips.
The frequency dependant characterization of the impedance magnitude is plotted in
Figure 4A for the conventional and improved chips. The impedance is seen to be higher in the conventional chip compared to the improved chip, meaning that the current is larger between the electrodes in the improved design, which is in accordance with the FEM simulations. The frequency dependent phase angle is plotted in
Figure 4B and follows the expected trend for the channel design [
3]. The small discrepancy is expected to be due to variations in the fabrication of the two chips.
Figure 5A shows the time dependant noise level of the two different electrode layouts. The root mean square of the noise from the different signals are with an
n value of 2000, conventional real = 2.4585 × 10
−10 A, conventional imaginary = 1.90681 × 10
−10 A, optimized real = 2.9216 × 10
−10 A, and optimized imaginary = 4.1071 × 10
−10 A.
Figure 5B–F show single transitions of the different detected particles for the two systems. It is seen that the shape of the time dependant transition is slightly elongated and more defined in the optimized chip design due to the different geometrical layout of the channel and electrodes. This is due to the lower flow velocity in the enlarged regions of the channel and the concentrated electric field in the constricted area between the electrodes, respectively.
Figure 5.
Signal plots in time of the transitions of the different particles in the conventional and optimized chip design. (A) Plot of the noise level of the different chip systems. (B) 1 μm bead in the conventional chip. (C) 2 μm bead in the conventional chip. (D) 0.5 μm bead in the optimized chip. (E) 1 μm bead in the optimized chip. (F) 2 μm bead in the optimized chip.
Figure 5.
Signal plots in time of the transitions of the different particles in the conventional and optimized chip design. (A) Plot of the noise level of the different chip systems. (B) 1 μm bead in the conventional chip. (C) 2 μm bead in the conventional chip. (D) 0.5 μm bead in the optimized chip. (E) 1 μm bead in the optimized chip. (F) 2 μm bead in the optimized chip.
Figure 6 shows a histogram of the current response of samples with polystyrene beads with a diameter of 2 μm and samples of mixed polystyrene beads with diameters of 0.5, 1, and 2 μm, while they pass the electrodes in the conventional chip (green and red) and improved chip (blue and magenta) with an excitation signal on the center electrode of 3 V AC at 231 kHz.
As can be seen in
Figure 6, the data obtain from the conventional chip yields a well-defined distribution of 2 μm beads (green) with a peak around 60 nA. It is also seen that the mixed beads sample injected through the conventional chip (red) split into two distinct distributions; one with a peak around 60 nA, and one with a peak closer to 10 nA. Since one of the populations coincides with the distribution of 2 μm beads it is assumed the beads with diameters of 0.5 μm and 1 μm make up the population at the low end of the current spectrum. It is also possible, however, that the 0.5 μm beads are not detected at all. Based on this, it is concluded that it is not possible to differentiate between polystyrene beads with diameters of 0.5 μm and 1 μm using the conventional electrode design. The histogram distributions of the recorded current response from the optimized chip (blue and magenta) are also shown in
Figure 5. For the optimized chip it is seen that the mixed beads sample (magenta) split into three distributions indicating the system’s ability resolve 0.5, 1, and 2 μm beads. The two histograms also show, as expected from the FEM simulations, that the current response from identical particles is larger for the improved chip. Further, the larger spread of the 2 μm beads in the histogram from the improved chip arise either from the larger variation in the current density or comes from the chip’s better size resolution. This is in agreement with the simulations shown in
Figure 3. The analysis of the distributions seen in
Figure 6 is shown in
Table 1. The measured current is expected to scale linearly with particle volume. This means that the current is expected to increase by a factor of eight when moving from beads with a radius 0.5 μm to beads with a radius of 1 μm and, similarly, another factor of eight when moving from a radius of 1 μm to a radius of 2 μm. However, the values in
Table 1 do not support this assumption. It is seen that the mean, median and mode for the analyzed population scale with a factor of five to seven when the radius of the detected beads is doubled. Since the measurements were performed on mixed samples (
i.e., the measurement on all bead sizes were performed under identical conditions) for both chip designs, it is very unlikely that this discrepancy is due to measurement error. It is speculated that this is due to the vertical variation of the electric field in the conventional chip and the vertical and horizontal variation in the optimized chip.
Figure 6.
Histogram of the maximum current response of samples with beads with a diameter of 2 μm for the optimized (blue) and conventional (green) chip designs and samples with mixed beads with diameters of 0.5, 1, and 2 μm for the optimized (magenta) and conventional (red).
Figure 6.
Histogram of the maximum current response of samples with beads with a diameter of 2 μm for the optimized (blue) and conventional (green) chip designs and samples with mixed beads with diameters of 0.5, 1, and 2 μm for the optimized (magenta) and conventional (red).
Table 1.
Calculated mean, median, mode and standard deviation of the two populations obtained from the measurements of the mixed sample in the conventional chip (red) and the three populations obtained from the measurements of the mixed sample in the optimized chip (magenta) shown in
Figure 6.
Table 1.
Calculated mean, median, mode and standard deviation of the two populations obtained from the measurements of the mixed sample in the conventional chip (red) and the three populations obtained from the measurements of the mixed sample in the optimized chip (magenta) shown in Figure 6.
Electrode layout | Sample | Mean (nA) | Median (nA) | Mode (nA) | Standard deviation (nA) |
---|
Conventional | Mixed—2 μm PS beads (red) | 58.8 | 58.5 | 58 | 3.8 |
Conventional | Mixed—1 μm PS beads (red) | 8.4 | 8.5 | 10 | 3.1 |
Optimized | Mixed—2 μm PS beads (magenta) | 84.3 | 83.3 | 81 | 6.3 |
Optimized | Mixed—1 μm PS beads (magenta) | 16.4 | 15.7 | 15 | 3.6 |
Optimized | Mixed—0.5 μm PS beads (magenta) | 3.4 | 3.0 | 3 | 1.4 |
In comparison with the FEM simulations it seen that the improved sensitivity of the system corresponds well with what is expected.