A Reaction-Advection-Diffusion Model to describe Non-Uniformities in Colorimetric Sensing using Thin Porous Substrates

This study develops and validates a reaction-advection-diffusion model to explain non-uniform color distributions and ring-like patterns in paper-based colorimetric sensors, demonstrating that mass transport and reaction dynamics alone can drive spatial variations without evaporation, thereby providing critical insights for optimizing sensor design and protocols.

The Gist

The Big Picture: Why Your Paper Sensor Looks "Messy"

Imagine you are using a piece of special paper to test water for lead or nitrate. You put a drop of water on the paper, and it's supposed to turn a uniform red color if the chemical is there. But often, instead of a nice, even red circle, you get a weird pattern: maybe a dark ring in the middle, a pale center, or even multiple rings like a target.

This paper asks: "Why does the color spread out unevenly, and how can we fix it?"

Most people think this happens because of the "Coffee Ring Effect" (like when coffee dries and leaves a ring of dirt at the edge). But the authors say: "No, that's not the main culprit here." Even if the liquid soaks into the paper instantly (before it can evaporate and form a coffee ring), the colors still get messy.

The real reason is a complex dance between flow (how the liquid moves), reaction (how the chemicals mix), and stickiness (how much the chemicals stick to the paper).

---

The Two Main Characters: The "Drop" and the "Paper"

To understand the experiment, imagine the paper sensor as a stage with two actors:

1. The Embedded Actor (The Paper): One chemical is already baked into the paper fibers. Let's call this the "Paper Guest."
2. The Droplet Actor (The Water): The other chemical is in the liquid drop you place on top. Let's call this the "Visitor."

When the Visitor soaks into the paper, they meet the Paper Guest and react to create a colored "Baby" (the product). The problem is figuring out where these Babies are born.

The Two Ways to Set the Stage

The researchers studied two different setups:

• Setup A (Reagent-Embedded): The Paper Guest is the chemical you want to detect. The Visitor is the test liquid. (Think of a pregnancy test strip where the test line is pre-printed).
• Setup B (Analyte-Embedded): The Paper Guest is the test liquid (pre-concentrated on the paper), and the Visitor is the chemical added to trigger the color. (Think of catching a fish in a net, then adding bait to see if it's the right kind of fish).

The Two-Stage Dance

The researchers realized the process happens in two distinct phases, like a two-act play:

Act 1: The Rush (Imbibition)
When you drop the liquid, it doesn't just sit there. It gets sucked into the paper like a sponge.

• The drop shrinks (the Visitor gets smaller).
• The wet spot on the paper grows (the "Paper Guest" gets visited).
• The Magic: As the liquid rushes outward, it carries the chemicals with it. If the chemicals react too fast, they make color right where they meet. If they react slowly, they get carried further out before making color. This creates the rings.

Act 2: The Calm (Diffusion)
Once the drop is fully soaked in, the liquid stops moving. Now, the chemicals just slowly drift around (diffuse) and finish reacting. The pattern created in Act 1 is mostly locked in by this point.

Why Do the Rings Form? (The "Traffic Jam" Analogy)

Imagine a highway (the paper) where cars (chemicals) are driving.

• The Visitor enters from the center and drives outward.
• The Paper Guest is parked all along the highway.

Scenario 1: The Fast Reaction (The Traffic Jam)
If the cars react instantly when they see a parked car, they crash (react) right at the center. The center gets crowded with "accidents" (color), and the outer roads stay empty. You get a solid dark center.

Scenario 2: The Slow Reaction (The Overtake)
If the cars drive past the parked cars without reacting immediately, they travel further out. But eventually, they run out of fuel (reactants). They might all crash at the same specific distance from the center. This creates a perfect ring in the middle of the paper, with a pale center and pale edges.

Scenario 3: The "Coffee Ring" Myth
Usually, we think rings form because the liquid evaporates at the edge, dragging dirt there. But this paper proves that even if the liquid soaks in faster than it evaporates, you still get rings. The ring is caused by the competition between how fast the liquid moves and how fast the chemicals react.

The "Stickiness" Factor

The paper also looked at how "sticky" the chemicals are.

• Mobile Chemicals: They flow freely with the water.
• Immobile Chemicals: They get stuck to the paper fibers (like Velcro).

The Discovery:

• If you want a uniform color in Setup A, you should make the "Paper Guest" sticky (immobile). This stops it from running away to the edges.
• If you want a uniform color in Setup B, it depends on whether the product (the color) is sticky or mobile.

The "Goldilocks" Solution: Paper Thickness and Porosity

The researchers found that the type of paper matters a lot:

• Thin, dense paper: The liquid moves slowly. The reaction happens in a small spot. You get intense color, but it's very uneven (messy rings).
• Thick, porous paper: The liquid moves fast and spreads out. The reaction is spread over a larger area. You get a very uniform color, but it is fainter (less intense).

The Trade-off: You have to choose between a bright, messy ring or a faint, perfect circle.

The Real-World Test

To prove their theory, they did real experiments:

1. Lead Detection: They put lead on the paper and added a reagent drop. The model predicted exactly where the dark rings would appear, matching their photos perfectly.
2. Nitrite Detection: They did the reverse (reagent on paper, nitrite drop). Again, the model predicted the weird patterns (sometimes multiple rings) that they saw in the lab.

The Takeaway

This paper is like a recipe book for better sensors. It tells engineers:

• Don't just blame evaporation for messy colors; look at how fast the liquid soaks in and how fast the chemicals react.
• If you want a uniform sensor, use thicker, more porous paper and maybe "glue" (immobilize) one of the chemicals to the paper so it doesn't run away.
• You can actually predict where the rings will form by doing the math before you even build the sensor.

In short: To get a perfect color on paper, you have to control the traffic flow of the chemicals, not just wait for the water to dry.

Technical Summary

1. Problem Statement

Colorimetric paper-based sensors (µPADs) are widely used for point-of-care diagnostics due to their low cost and simplicity. However, a major limitation is the non-uniform distribution of the colored reaction product, which compromises measurement accuracy and reliability.

• The Gap: While the "coffee-ring effect" (caused by evaporation-driven flow on impermeable surfaces) is well-known, it does not fully explain non-uniformities in porous substrates where imbibition (liquid absorption) dominates over evaporation.
• Specific Challenge: Existing models fail to explain the emergence of ring-like patterns at intermediate radial positions (not just at the periphery) or the formation of multiple rings. The interplay between hydrodynamics, mass transport, and reaction kinetics within the porous substrate remains unclear.

2. Methodology

The authors developed a theoretical framework based on reaction-advection-diffusion equations to model the spatial distribution of a colored product formed by the reaction between two species:

• Species S: Embedded uniformly in the porous substrate.
• Species D: Delivered via a sessile droplet that imbibes into the substrate.
• Product P: The colored result of the reaction ($D + S \to P$).

Model Stages

The process is divided into two distinct stages:

1. Stage 1 (Droplet Imbibition):
   • The droplet penetrates the substrate, creating a moving boundary problem.
   • Domain I ($0 \le r < R_d$): The region under the droplet acts as a Mixed Flow Reactor (MFR) with uniform concentration.
   • Domain II ($R_d \le r \le R_p$): The wetted region outside the droplet base where radial flow occurs, driven by capillary pressure.
   • Assumptions: Imbibition is much faster than evaporation (neglecting evaporation effects); vertical saturation is instantaneous due to thin substrate; flow is axisymmetric.
2. Stage 2 (Post-Imbibition):
   • Occurs after the droplet is fully absorbed.
   • Liquid is stationary; transport is governed solely by diffusion and reaction.

Key Features

• Mobility Factor ($m_i$): A parameter ($0 \le m_i \le 1$) is introduced to account for species-substrate interactions (e.g., adsorption, precipitation, or solubility). $m_i=1$ implies a fully mobile species; $m_i=0$ implies a fully immobile species.
• Configurations: The model analyzes two setups:
  • Reagent-Embedded (RE): Reagent is in the substrate; Analyte is in the droplet.
  • Analyte-Embedded (AE): Analyte is preconcentrated in the substrate; Reagent is in the droplet.
• Numerical Solution: The governing non-linear PDEs are solved using the Method of Lines (discretizing space with finite differences) in MATLAB.

3. Key Contributions

1. Mechanism of Intermediate Rings: The study demonstrates that ring-like patterns can form at intermediate radii even in the absence of evaporation. This challenges the conventional view that coffee rings are solely evaporation-driven.
2. Role of Transport-Reaction Competition: The paper establishes that the spatial distribution is governed by the competition between the reaction rate and advective/diffusive transport within the porous medium.
3. Generalized Framework: The model incorporates a mobility factor, allowing it to predict outcomes for various chemical scenarios (e.g., precipitating products vs. soluble products) and substrate types.
4. Validation: The model is validated against experimental data for Lead (Pb²⁺) detection (AE configuration) and Nitrite (NO₂⁻) detection (RE configuration), successfully capturing complex patterns including single rings, double rings, and core-dominant profiles.

4. Key Results and Findings

Hydrodynamics

• Imbibition occurs significantly faster than evaporation (timescale ratio $\approx 4.5 \times 10^3$), justifying the neglect of evaporation in the model.
• The droplet radius decreases while the wetting radius increases until the droplet is fully absorbed.

Influence of Concentration Ratio ($M_A$)

• Stoichiometric Ratio ($M_A \approx 1$): Leads to incomplete reaction and the formation of a ring-like pattern at an intermediate radius in both configurations.
• Excess Reagent ($M_A < 1$):
  • AE Configuration: The reaction peak shifts inward (towards the center) as reagent availability increases.
  • RE Configuration: The profile becomes flatter and more uniform as reagent excess increases.

Influence of Substrate Properties

• Thickness & Porosity: Thicker and more porous substrates lead to greater uniformity in color distribution but result in reduced color intensity (lower peak concentration). This suggests a trade-off between homogeneity and sensitivity.
• Multiple Rings: The model explains the emergence of multiple rings under specific conditions where local maxima in reaction rates shift during the process.

Effect of Species Immobilization

• Immobile Product (e.g., Precipitate):
  • AE: Immobilizing the embedded analyte (S) significantly improves uniformity.
  • RE: Immobilizing the embedded reagent (S) is less effective and can lead to radially decreasing profiles.
• Mobile Product:
  • RE: Immobilizing the reagent (S) effectively reduces color gradients.
  • AE: Immobilization has little effect if the product is mobile.

5. Significance and Implications

• Design Optimization: The study provides actionable guidelines for optimizing µPADs. For instance, selecting thicker, more porous papers can improve uniformity, while adjusting the reagent-to-analyte ratio can control the location of the color signal.
• Protocol Development: Understanding that non-uniformity arises from internal transport-reaction dynamics (not just evaporation) allows researchers to design better protocols, such as preconcentration steps or specific immobilization strategies.
• Broad Applicability: While focused on circular sensors, the framework is extendable to lateral flow assays and other complex geometries, offering a fundamental understanding of reactive transport in porous media.

In conclusion, this work shifts the paradigm from viewing non-uniformity as an evaporation artifact to understanding it as a complex interplay of advection, diffusion, and reaction kinetics within the porous substrate, providing a robust tool for the next generation of reliable colorimetric sensors.