2.1. Problem Formulation
Our study addresses the pressing need for continuous, real-time monitoring of pro-inflammatory biomarkers in human sweat. Such monitoring can provide valuable insights into an individual’s physiological state and early warning signs of inflammation-related conditions. To achieve this, we have designed a novel wearable microfluidic immunosensor patch, illustrated in
Figure 1. This biosensor features a dual-unit structure to optimize the efficiency and accuracy of biomarker detection.
Unit (I), as detailed in
Figure 1(b), serves as the initial stage of sample introduction in the biosensor. It incorporates an array of coils that generate a localized magnetic field. This magnetic field strategically traps and pre-concentrates magnetic nanoparticles (MNPs) functionalized with antibodies specific to the target pro-inflammatory biomarkers within the microfluidic channel. The pre-concentration step significantly enhances the biosensor’s specificity by increasing the probability of successful interactions between the MNPs and target biomarkers present in the sweat sample. In addition, the coils induce microfluidic mixing, ensuring a homogeneous distribution of both biomarkers and MNPs throughout the microfluidic channel. This uniform distribution creates an optimized environment for the efficient formation of immunocomplexes, which occur when the target biomarkers bind to their corresponding antibodies on the MNPs. By promoting efficient immunocomplex formation, the coils within Unit (I) accelerate the overall assay speed, enabling faster detection of biomarkers.
Unit (II), depicted in
Figure 1(c), acting as the selective detection zone of the biosensor, features an upper coil that captures the biomarker-tagged MNPs from Unit (I). The captured MNPs are then precisely guided into a microfluidic channel embedded within a dielectric domain defined by dual parallel electrodes. This configuration forms a capacitive structure, where the presence of biomarker-tagged MNPs induces a measurable change in capacitance. This change occurs because the dielectric material within the domain interacts with the electric field differently when MNPs are present due to their unique electrical properties. By precisely measuring this capacitance shift, the sensor indirectly quantifies the concentration of target biomarkers. The strategic placement of metallic electrodes within the dielectric domain plays a crucial role in maximizing the signal-to-noise ratio, ensuring accurate and reliable detection of even minute quantities of biomarkers.
The compartmentalized design, featuring distinct units for binding and detection, offers notable advantages over conventional single-unit biosensors. This separation of functions provides greater control over the assay conditions, minimizes cross-contamination and non-specific interactions, and ultimately enhances the overall performance and reliability of the biosensor. Consequently, this innovative design is well-suited for the continuous, real-time monitoring of pro-inflammatory biomarkers in wearable applications.
To validate and optimize the design of our proposed biosensor, we developed a comprehensive 3D model (
Figure 1(d)). This model encompasses three distinct domains, each simulating a specific physical process occurring within the device:
Magnetic Field Domain: This domain focuses on calculating the magnetic field distribution generated by the coils. Accurate prediction of the magnetic field is crucial for understanding how MNPs are manipulated and trapped within the microfluidic channels.
Fluid Dynamics Domain: This domain simulates the behavior of the biofluid (sweat) within the microfluidic channels, considering the complex interactions between the fluid, MNPs, and channel walls. This simulation helps optimize channel geometry and flow conditions for efficient biomarker capture.
Electrical Domain: This domain analyzes the electrical behavior of the system, particularly the capacitive response resulting from the presence of MNPs within the dielectric domain. Precise modeling of the electrical properties ensures accurate biomarker quantification.
By utilizing the finite element method (FEM) in COMSOL Multiphysics V 6.0, we numerically solve the governing equations for the magnetic, fluidic, and electric fields. This comprehensive simulation approach offers valuable insights into the biosensor’s operational dynamics. We can use these insights to iteratively refine the design, ultimately improving sensitivity and specificity, and bringing us closer to realizing a wearable biosensor that empowers individuals with real-time health monitoring capabilities.
2.2. Physics and Mathematical Framework of the Biosensor
Expanding upon the introduction and design description, we delve into the fundamental physical and mathematical principles underpinning the microfluidic immunosensor’s operation. This in-depth analysis, focusing on the three critical domains visualized in
Figure 2(a), provides a rigorous understanding of the sensor’s functional mechanisms.
Domain 1: Planar Coils
Domain 1 focuses on characterizing the magnetic field generated by the coils within the biosensor. To optimize magnetic field strength while minimizing power consumption, we evaluated four distinct copper planar coil designs (
Figure 2(b) and
Figure 2(c)), each varying in outer radius and turn count (R500 to R2000). All designs maintained consistent copper wire dimensions (width (w) = 10
m, height (h) = 10
m, and a separation between adjacent wires (s) = 10
m) to ensure a fair comparison.
We harnessed Maxwell’s equations to accurately predict the magnetic field distribution around the coils. The magnetic vector potential (A), which enables the derivation of the magnetic flux density (B) and magnetic field strength (H), is defined as: [
32]
Furthermore, the Ampere-Maxwell equation, incorporating current density (J), provides a comprehensive framework for modeling the electromagnetic behavior of the coils: [
32]
The relationship between magnetic flux density (B), magnetic field strength (H), and the medium’s properties is described by: [
32]
where
is the magnetic permeability in free space, and
is the relative permeability of the medium.
For accurate simulation of electromagnetic fields within the coil devices, we employed quadratic discretization in the Finite Element Method (FEM). This method, utilizing second-degree polynomial approximations, offers higher accuracy and a better representation of complex magnetic field distributions compared to traditional linear discretization methods.
To quantify coil efficiency, we used the power merit factor (), a critical metric relating maximum magnetic field strength () to electrical power consumption ((P)), which is the product of voltage (U) and current (I). By carefully selecting coil dimensions and operational parameters through simulation, we aimed to strike an optimal balance between field intensity and power consumption for long-term device reliability and performance.
Domain 2: Microfluidic Platform
In Domain 2, our focus shifts to the microfluidic platform, where the dynamics of sweat flow and MNP transport are crucial. We model the flow of sweat, primarily composed of water, through the microchannel using the incompressible Navier-Stokes equations
The behavior of MNPs within this flow is of paramount importance. These nanoparticles experience magnetofluidic forces arising from the interaction between their magnetic properties and the external magnetic field gradient. By meticulously accounting for these forces, along with viscous drag forces (modeled using Stokes’ law), we gain valuable insights into the MNP trajectories and their interactions with the channel walls. This knowledge is essential for optimizing capture efficiency and ensuring controlled MNP delivery to the detection zone.
The governing incompressible Navier-Stokes equation for our fluid dynamics simulation, incorporating the velocity vector
, fluid density
, pressure gradient
, dynamic viscosity
, and external forces
, is given by: [
33]
The magnetofluidic force
on the MNPs is quantified by: [
34]
This equation indicates that the magnetofluidic force is proportional to both the cube of the particle diameter
D and the gradient of the squared magnetic field intensity
. For the drag force
, acting on a spherical MNP with diameter
D, we apply Stokes’ law: [
35]
where
denotes the fluid’s viscosity (
),
the mass of the particle (kg),
the particle velocity vector (m/s), and
the density of the particle (
).
To solve the complex equations governing particle motion within the fluid flow, we employ the Generalized alpha method (GAM). This advanced numerical technique excels at managing the temporal evolution of the system, meaning it can accurately capture how the motion of the MNPs changes over time within the dynamic fluidic environment. GAM also effectively suppresses spurious oscillations, which are numerical artifacts that can arise in simulations and compromise the accuracy of the results. By mitigating these oscillations, GAM ensures the stability of our simulations and the reliability of the data they produce.
Domain 3: Capacitive Electrodes
The final domain, Domain 3, focuses on the capacitive sensing mechanism that underlies the biosensor’s detection capabilities. By strategically arranging metallic electrodes within the microfluidic channel, we create a capacitor that is highly sensitive to changes in capacitance induced by the presence of MNPs, as shown in
Figure 2(d). This direct-field capacitive sensing approach offers significant advantages in terms of sensitivity.
The capacitance (
C) of the system is primarily governed by the dielectric properties of the materials between the electrodes, as described by the equation: [
36]
where
represents the vacuum permittivity,
is the relative dielectric constant of the material between the electrodes,
A denotes the electrode area (
), and
d is the electrode separation distance (
).
As MNPs accumulate in the detection zone, they alter the effective dielectric constant
), leading to a measurable change in capacitance (
), quantified by:
This expression highlights the sensor’s sensitivity to changes in biomarker concentrations, as these changes are reflected in MNP-induced variations in dielectric properties.
To further refine our model, we consider the microfluidic channel’s structure and estimate the equivalent capacitance (
) of the biosensor, taking into account the PDMS walls and the variable channel capacitance (
Figure 2(d)):
Here,
and
represent the dielectric constants of PDMS and the channel, respectively. This comprehensive model, integrating both direct-field capacitive sensing and magnetic trapping, provides a robust framework for understanding and optimizing the detection and quantification of MNP-tagged biomarkers.
In practical terms, if the initial capacitance without MNPs is (
) , and the capacitance changes to (
C) after the introduction of MNP-biomarker complexes, the change in capacitance (
) can be calculated using Equation
8. This (
) value serves as a unique fingerprint for the specific biomarker concentration within the sweat sample.
To translate this fingerprint into a quantitative measurement of biomarker concentration, we utilize a pre-established calibration curve. This curve is generated by meticulously measuring the change in capacitance () across a range of known concentrations of the target biomarker during the biosensor’s development phase. By referencing this curve, we can accurately determine the unknown concentration of the biomarker in a new sweat sample. This is achieved by simply matching the measured () value from the new sample to the corresponding concentration on the calibration curve. This approach allows for precise and reliable determination of biomarker levels in real-world samples.