## Abstract

Adsorption-based microfluidic sensors are promising tools for biosensing. Advanced mathematical models of time response and noise of such devices are needed in order to improve the interpretation of measurement results, and to achieve the optimal sensor performance. Here the mathematical models are presented that take into account the coupling of processes that generate the sensor signal: adsorption–desorption (AD) of the target analyte particles on the heterogeneous sensing surface, and mass transfer (MT) in a microfluidic chamber. The response kinetics and AD noise (which determines the ultimate sensing performance) of protein biosensors are analyzed, assuming practically relevant analyte concentrations, sensing surface areas and MT parameters. The condition is determined under which MT significantly influences the sensor characteristics relevant for reliable analyte detection and quantification. It is shown that the development of improved mathematical models of sensor temporal response and noise can be used as one of strategies for achieving better sensing performance.

### Keywords

- microfluidic biosensor
- surface heterogeneity
- mass transfer
- adsorption–desorption noise
- mathematical model

## 1. Introduction

Microfluidic adsorption-based biosensors are promising devices for real-time, in-situ and low-cost analysis of samples taken from the environment, food or living organisms, enabling detection of the presence and measurement of the amount of target biological specimens: biomolecules (such as proteins or DNA fragments), microorganisms, or other biological structures [1, 2, 3]. Such sensors are highly sensitive, capable of operation with small sample quantities, and also small, lightweight and energy efficient, thus being especially suitable for autonomuous, portable and distributed sensing applications [4, 5, 6]. Due to such characteristics, development of microfluidic sensors is very significant for environmental protection, medicine, agriculture, food inspection, public healthcare and security, and other fields, where they can substitute large and expensive laboratory equipment, typically located far from the place where the samples are taken for the analysis. They also enable the development of new fields of biosensor applications, such as the personalized medical point-of-care diagnostics, telemedicine, wearable sensors etc. [7, 8, 9].

In order to utilize the great potential of adsorption-based microfluidic biosensors for practical applications, research is performed aiming to enable optimization of their performance. The increase of sensitivity and response rate of the sensors, better selectivity, lowering of the minimal detectable concentration, higher reliability of measurement results ant their more accurate interpretation, as well as the research and development of new measurement methods, which enable obtaining of more information about one or multiple adsorbed substances at the same time, are of great practical significance. In that sense, of particular interest is to know the dependence of the temporal response and the sensor noise on the parameters of the sensing element, the measurement system, and the experimental conditions, which requires the development and application of mathematical models that take into account physical processes and phenomena relevant for generating the response and its fluctuations. The adsorption and desorption (AD) processes are inevitably taken into account in modeling of the response and noise, since they are fundamental for sensor operation, and are also the source of adsorption–desorption noise, which sets the fundamental limits of detection and quantification of the analyte. In structures of micrometer and nanometer dimensions, AD noise can dominantly determine the values of minimal detectable and quantifiable signal, as well as other limiting sensor performances, especially in the case of low analyte concentration [10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20]. In various cases, apart from the AD process of the target analyte, different additional processes influence the sensor response kinetics and noise, so it is necessary to take them into account in mathematical models.

The objective of this work is to present mathematical models of the temporal response and adsorption–desorption noise of microfluidic adsorption-based biosensors, that take into account the processes responsible for the generation of the sensor signal: adsorption–desorption of the target analyte particles on the heterogeneous sensing surface, and mass transfer (MT) in a microfluidic chamber. Section 2 presents the mathematical models of the sensor temporal response, while Section 3 describes the AD noise models, developed for the cases considered in Section 2. In Section 4, the results will be presented of the temporal response analysis, and of the analysis of AD noise of protein biosensors, both performed by using of the models presented in Sections 2 and 3, assuming practically relevant analyte concentrations, sensing surface areas and MT parameters. The conclusions will be summarized in Section 5.

## 2. Mathematical models of sensor temporal response

In adsorption-based biosensors (e.g. SPR (* Surface Plasmon Resonance*), resistive graphene-based, CNT (

*) or NWFET (*Carbon NanoTube

*), SAW (Surface Acoustic Wave), FBAR (*NanoWire Field Effect Transistor

*), microcantilever sensors) detection of the target analyte and measurement of its concentration are based on the change of a measurable parameter of the sensing element, caused by analyte adsorption on the active surface [21, 22, 23, 24, 25, 26, 27, 28]. Namely, the adsorption leads to the change of some of the sensing element’s physical parameters (e.g. the mechanical strain or mechanical structure’s mass, the effective density of the surface layer or its conductivity, the refraction index, the distribution of electric charges on the surface), which changes at least one measurable parameter of the sensor (e.g. the deflection or the resonant frequency of the mechanical structure; the amplitude, frequency or phase of surface or bulk acoustic waves in mechanical sensors; the resistance or current in electrical sensors; the intensity, reflection, transmission, or absorption of light in optical sensors).*thin Film Bulk Acoustic wave Resonator

Here, the term “temporal response of adsorption-based sensors” will be used to denote the temporal change of a certain physical parameter of the sensing element, which is induced by analyte adsorption. At any given time, that change is determined by the number of analyte particles bound to the sensing surface, and that number depends on the target analyte concentration in the sample. This enables the measurement of the concentration of the target substance to be performed by measuring the change of the physical parameter. Therefore, the mathematical model of the sensor temporal response is based on the model of the time evolution of the number of adsorbed particles as a function of analyte concentration, assuming that there is a known (preferably linear) relation between the two quantities.

In the analyses of microfluidic biosensors time response, which have been published in the literature, it is often assumed that the change of the number of adsorbed particles occurs only due to the AD process of the target substance, and the interpretation of experimental results is performed according to that [29, 30]. This simplified interpretation of the events occurring on the sensing surface is justified under the conditions that ensure a negligible influence of other processes, and assume a homogeneous adsorption surface. Analyses that include some of the additional effects that also influence the current number of adsorbed target particles are less abundant in the existing literature. For example, the mass transfer process of target adsorbate particles in a reaction chamber, toward or away from adsorption sites on the sensing surface, coupled with the AD process, is considered in [31, 32, 33]. Competitive AD processes of target and competitor substances are analyzed in [34, 35]. An example of the analysis encompassing multiple phenomena is the research of adsorption processes of two or three analytes on the surface of the same sensor, coupled with mass transfer processes of corresponding particles [36, 37, 38]. Sensing surface heterogeneity and mass transfer processes in biosensors are considered within the analysis of sensor models response in [30, 36, 39], where the need is emphasized for taking into account both of these factors simultaneously while creating experiments and interpreting the measurement results.

In this section, first the starting system of equations will be presented that model the change of the number of adsorbed particles on a sensing surface with an arbitrary number of different adsorption sites, and take into account mass transfer processes of analyte particles in the sensor chamber. Subsequently, simplified physical models will be defined and presented by suitable equations for certain cases of practical significance, in a similar way as in [39]. Also, conditions will be defined under which the application of the approximate models is justified. The case of adsorbing surface with two types of adsorption sites will be particularly considered.

### 2.1 Starting equations

It is assumed that the sensing element of a sensor is in a microfluidic reaction chamber of a rectangular cross-section, with the sample to be analyzed flowing through it. In Figure 1 (left) a schematic representation is given of the chamber and the sensing element, with the designations of dimensions and coordinate axes. During the laminar fluid flow, a parabolic velocity profile is formed in the chamber (the profile belongs to the * z* = const. plane, and it is constant in the

*-axis direction in the adsorbing surface zone).*z

On the sensor’s active surface the recognition and binding of the target analyte occurs. Adsorption–desorption and mass transfer are the key processes for binding of particles to the sensing surface. The former is the process of binding of the analyte particles to surface adsorption sites due to a certain affinity, and unbinding from them. The latter encompasses the processes (convection and diffusion) by which the analyte particles are transported through the microfluidic chamber to or from the binding sites. The coupling of these processes determines the spatial and temporal dependence of the analyte concentration in the chamber.

In the analysis, it is assumed that the adsorption occurs in a single layer only, that only one analyte particle can be adsorbed on an adsorption site, and that the probability of adsorption on any given site (or desorption from it), does not depend on the occupancy of adjacent sites. In the case of a homogeneous sensing surface (in the sense of the affinity toward the target analyte particles), the AD process is characterized by a single pair of adsorption and desorption rate constants, * k*and

_{a}

*.*k

_{d}

The active surface of affinity-based sensors is often not comprised of uniform adsorption sites to which the particles of the target analyte bind [30, 40]. The reason for this can be a non-uniform surface morphology or chemical composition, or the existence of multiple possible binding sites for analyte particles due to the material structure on the molecular scale. In more detail, the presence of various defects, irregularities, cavities, pores, impurities, contaminants, different functional groups on the surface, etc. are some of numerous features that constitute the surface heterogeneity. Materials in the form of flakes with specific binding sites on edges and basal planes, such as liquid-phase exfoliated graphene [41], are also an example of a heterogeneous sensing surface. When a surface is functionalized with specific capturing probes for the target analyte, adsorption sites heterogeneity can be manifested through nonspecific binding of analyte particles to the surface. An adsorption surface can also be heterogeneous due to uneven binding of functionalizing entities (e.g. non-uniform orientation of capturing proteins when attached to the surface), which influences the affinity toward the analyte, and the efficiency of its binding. Whatever the reason, a heterogeneous sensing surface can be characterized by different affinities of different surface sites for the analyte particles binding. The analyte’s AD process is then characterized by a certain distribution of adsorption and desorption rate constants across the surface, and it is described by a model that takes into account the surface non-uniformity.

If we assume that there is * n*types of adsorption sites on the surface, the surface heterogeneity can be described by a discrete set of

*values of adsorption energies for the given analyte, thus yielding*n

*pairs of adsorption and desorption rate constants (*n

*,*k

_{ai}

*),*k

_{di}

*= 1,2…*i

*. Adsorption of a single analyte on*n

*types of adsorption sites, randomly distributed across the adsorbing surface of a microfluidic sensor, is mathematically described by*n

*+ 1 equations, and by boundary and initial conditions. One of the equations is the convection-diffusion equation (for a microfluidic chamber of a rectangular cross-section, whose illustration is shown in Figure 1 (left))*n

Here, * D*is the diffusion coefficient of analyte particles, and

*is the flow velocity of the sample. The remaining*v

*equations pertain to the processes of reversible binding of particles to adsorption sites of different types. Each of these processes, which are coupled in a general case, can be regarded as one of*n

*components of a complex AD process. It is known that the Langmuir adsorption model assumes the uniformity of adsorption sites. In the case of heterogeneous surface sites, it can be assumed that the adsorbent surface consists of a collection of locally homogeneous surfaces, and that the adsorption on each of them can be considered as Langmuir adsorption. All the parts of the sensing surface*n

*with the adsorption sites of the type*S

*, constitute the surface*i

*of area*S

_{i}

*. The above-mentioned set of*A

_{i}

*equations can thus be written as*n

Here * η*is the surface density of the analyte particles adsorbed on sites of the

_{i}

i

^{th}type,

η

_{max,i}is the surface density of adsorption sites of the matching kind, and

*is the analyte concentration adjacent to the sensing surface. It is assumed that the adsorption sites are uniformly distributed across the surface*C

_{s}

*, so*S

_{i}

η

_{max,i}does not depend on the coordinates

*and*x

*.*z

The boundary and initial conditions are as follows: 1. at the entrance of the chamber (* x* = 0) the concentration is equal to that in the analyzed sample:

*(*C

*,0,*t

*,*y

*) =*z

C

_{0}; 2. at the exit from the chamber the continuity condition assumes free convection: ∂

*/∂*C

*= 0 for*x

*=*x

*; 3. the initial adsorbate particle distribution in the chamber is uniform:*L

_{c}

*(0,*C

*,*x

*,*y

*) =*z

C

_{0}, 4. in the zone of the sensing element, which is defined by the coordinates

*= 0,*y

x

_{1}≤

*≤*x

x

_{1}+

*,*L

_{s}

z

_{1}≤

*≤*z

z

_{1}+

*(*w

_{s}

*is the adsorption zone length, and*L

_{s}

*is its width), there is a balance between the diffusion flux in the direction perpendicular to the surface, and the net adsorption rate per unit surface*w

_{s}

The ^{th}Eq. (3) is valid on the parts of adsorbing surface with the adsorption sites of the ^{th} type; on the remaining part of the chamber’s surface where * y* = 0, and on the whole chamber’s surface where

*=*y

*, the zero flux condition is valid; 5. on the edges of every adsorption zone the flux equals zero (the adsorbed particles leave the adsorbing zone only by desorption); 6. the adsorption begins at the moment*h

_{c}

*= 0:*t

*(0,*η

*,*x

*) = 0.*z

For the determination of * C*(

*,*t

*,*x

*,*y

*) and*z

*(*η

_{i}

*,*t

*,*x

*), it is necessary to solve Eqs. (1) and (2) with the mentioned initial and boundary conditions, by using numerical methods. The number of adsorbed particles on the sites of the*z

i

^{th}type,

*(*N

_{i}

*), can be determined by integration of*t

*(*η

_{i}

*,*t

*,*x

*) over the sensing element surface. If the contribution of a particle adsorption to the senor signal does not depend on the place on the sensing surface where the particle is bound, but only on the type of the adsorption site, the sensor response will depend on the numbers of adsorbed particles per each site type.*z

When the total sensor response is the sum of contributions of * n*components, of which each is the product of the number of adsorbed particles of a given component,

*(*N

_{i}

*), and the corresponding weight factor*t

*(*w

_{i}

*= 1, 2, …,*i

*), the temporal response is*n

The weight factor * w*equals the average contribution of a single adsorbed particle of the

_{i}

i

^{th}component to the sensor response. However, for some types of adsorption-based sensors, the response can be represented by the simpler expression

For instance, in the case of resonant micro/nanocantilevers and acoustic wave mechanical sensors, particle adsorption changes the mass of the mechanical sensing structure, so that * w*is determined by the mass of a single analyte particle, thus it is independent on the type of the site where the particle was adsorbed. Also, in the case of plasmonic sensors, it is reasonable to assume that the mean refractive index change when an analyte particle is adsorbed is

*= (*w

*-*n

_{a}

*)/*n

_{e}

*, where*N

_{max}

*is the refractive index value of the analyte,*n

_{a}

*is the refractive index of the surrounding medium, and*n

_{e}

*is the total number of adsorption sites on the surface [42]. In these cases, it is justified to use Eq. (5).*N

_{max}

Numerical solving of Eqs. (1) and (2) requires the surface distributions of different adsorption sites to be known, which is rarely the case. Instead of that, approximations can be introduced in order to simplify the equations, which can even enable obtaining of the analytical solution valid for the transient regime and/or the steady state. The approximations are introduced based on the comparison of adsorption, convection and diffusion time scales.

### 2.2 Adsorption limited response kinetics

A significant reduction of mathematical complexity of the problem is possible when it is justified to assume that the concentration of the analyte in the reaction chamber is spatially uniform, constant in time, and equal to the concentration in the sample injected in the chamber. Eq. (3) is then reduced to * C*(

*,*t

*,*x

*,*y

*) =*z

C

_{0}= const. The spatial independence of the concentration

*o and*C

η

_{max,i}implies the uniformity of the surface density of adsorbed particles

*on the surface*η

_{i}

*, so the number of particles adsorbed on sites of the*S

_{i}

i

^{th}type is

*=*N

_{i}

*. The number of adsorption sites of that type is*η

_{i}A

_{i}

N

_{max,i}=

η

_{max,i}

*=*A

_{i}

*, where*ν

_{i}N

_{max}

*is introduced as a measure of abundance of these sites in the total number of sites on the sensing surface,*ν

_{i}

*. The model of the multicomponent AD process of a single analyte on the surface with*N

_{max}

*types of adsorption sites, derived from Eqs. (2), is therefore represented by the set of*n

*mutually independent equations*n

The analytical solution of each of Eqs. (6), with the initial condition * N*(0) = 0, is

_{i}

where * N*is the number of adsorbed particles in the steady state, which establishes with the time constant of the AD process,

_{ie}

*. For the steady state d*τ

_{AL,i}

*/d*N

_{i}

*= 0, i.e.*t

*=*a

_{i}

*. The sensor response is determined by Eq. (4) or Eq. (5).*d

_{i}

This model is applicable in the cases where the transfer flux toward the sites of each type is sufficiently greater than the adsorption flux, so that the binding kinetics is adsorption limited [31, 32, 37]. This means that the number of adsorbed particles on the given surface and at the given analyte concentration is determined only by the AD process parameters at any given time. This is the ideal case for obtaining the data about the AD process and its kinetics from the measured time response of the sensor. The transport-adsorption regime of this kind is typical for fast diffusing particles (such as gas molecules and certain biomolecules of small mass).

### 2.3 Mass transfer influenced response kinetics

When the transport flux of analyte particles is lower than the adsorption flux, or comparable with it, the time evolution of the numbers of adsorbed particles is affected by mass transfer processes [31, 32, 37]. The temporal and spatial change of the analyte concentration in the sensor chamber depends on the parameters of AD and mass transfer processes in the system of a given geometry. As the adsorption flux becomes more dominant, the change of the analyte concentration in the chamber becomes more pronounced due to the slow compensating influx of particles carried by transport processes in the space that is being depleted of particles due to their fast binding to the adsorbing surface. Here, the analysis will be focused on the situation in which a thin depleted layer is formed adjacent to the sensing surface, because it is common in many types of microfluidic sensors (especially those whose reaction chamber height is of the order of ≥10 μm [43]). It also enables approximations to be used, which lead to the simplified time response mathematical model, and also yield the analytical expression for the AD noise spectral density.

When the diffusion time scale is greater than the convection time scale (expressed by _{c}^{2}/* D* >

*/*L

_{c}

*, where*v

_{m}

v

_{m}is the mean convection velocity), not all the particles from the chamber volume can participate in the adsorption process, but only those from the layer of a certain thickness adjacent to the adsorbing surface. If the thickness of that layer is small compared to the chamber height and the length of the adsorbing zone, the temporally and spatially variable analyte concentration in the chamber can be approximated by the two-compartment model (TCM), whose applicability is experimentally confirmed [31, 44, 45]. According to that model, the chamber can be divided into two parts (compartments). One of them is the inner compartment, which is adjacent to the adsorbing surface, and contains a variable analyte concentration due to the depletion of analyte particles. The other, outer compartment, contains the analyte at the same concentration as it is in the sample injected in the chamber,

C

_{0}. Also, according to the model, all quantities are averaged over the adsorbing surface, and the transport between the two compartments is described by the mass transfer coefficient,

*. It is assumed that the zones*k

_{m}

*(with binding sites of the type*S

_{i}

*) consist of a multitude of smaller areas scattered over the sensing surface (as illustrated in Figure 1 (right)), so that different adsorption sites are mixed. The averaging of quantities over the sensing surface can thus approximate the actual conditions with sufficient accuracy, making the use of TCM justified. Furthermore, the model equates the rate of change of the total number of adsorbed particles with the net number of particles that enter the inner compartment in unit time, i.e.*i

where the mass transfer coefficient is given as * k* = 1.467(

_{m}

D

^{2}

*/(*v

_{m}

*))*L

_{s}h

_{c}

^{1/3}[31]. Since the total number of particles adsorbed on the entire sensing surface equals

and Eqs. (2) after averaging of the quantities over the surface yield

Eqs. (8)–(10) imply that the analyte concentration adjacent to the adsorbing surface is determined by the expression

The time evolution of the numbers of adsorbed particles on different types of sites is determined by Eqs. (10) and (11), i.e. by * n* + 1 coupled equations. This system of equations is significantly simpler than the starting system (Eqs. (1) and (2)), and can be efficiently solved by using numerical methods for the given initial conditions. The sensor response is then determined by Eq. (4) or (5). The steady-state values of the numbers of adsorbed particles are obtained from Eqs. (10) for d

*/d*N

_{i}

*= 0, and they are the same as in the case of adsorption-limited regime*t

The presented mathematical model is applicable for the thin depleted zone adjacent to the sensing surface, which is more likely to exist at lower * D*values, higher

*, and faster convection. The analysis of Eq. (11) shows that at*h

_{c}

_{s,TCM}≈_{0}, so the model given by Eq. (10) reduces to the model that is valid in the case of adsorption-limited binding (Eq. (6)). Therefore, Eq. (13) is the condition for the transfer flux to dominate over the adsorption flux (while the diffusion is slow compared to the convection). When the condition (13) is not satisfied, it is necessary to use the mathematical model given by Eqs. (10) and (11) for the response analysis of sensors in which a thin depleted zone is formed. Based on these considerations, it can be concluded that the equations obtained by using TCM have broader applicability than expected: they are valid both for mass-transfer influenced and for adsorption-limited (the case of sufficiently high * k*, according to Eq. (13)) binding kinetics.

_{m}

The simplest heterogeneous sensing surface contains two types of adsorption sites. The corresponding equations are

where * i*equals 1 or 2.

## 3. Mathematical models of adsorption-desorption noise

Due to the inherently stochastic nature of processes involved in the analyte particle binding-unbinding events on the sensing surface, the number of adsorbed particles randomly fluctuates, even after reaching the steady state, so it can be expressed as

where Δ* N*denotes fluctuations. The fluctuations of the number of adsorbed particles, Δ

*(*N

*), result in the fluctuations of the sensor’s time response, Δ*t

*(*R

*), which constitute the inevitable adsorption–desorption (AD) noise. Based on Eq. (4), these fluctuations are related as*t

assuming * n*AD processes on the sensing surface, where Δ

*denote the fluctuations of the number of particles that participate in the*N

_{i}

i

^{th}AD process.

Here, the goal is to obtain the analytical expression for the spectral density of sensor AD noise, when all the transient processes are finished, i.e. when an steady state is established. According to Eq. (16), the basis of this analysis is the analysis of fluctuations of the number of adsorbed particles around the steady-state values. Before presenting the theoretical models of AD noise of sensors with heterogeneous sensing surface, a short overview will be given of the already published results that include mathematical modeling of AD noise of sensors and other micro/nanodevices.

An insight into the existing literature shows that AD noise analyses have been usually limited to the consideration of individual phenomena pertinent to fluctuations of the number of adsorbed particles – e.g. fluctuations originating from a stochastic single-analyte AD process [46] or fluctuations due to surface diffusion of adsorbed particles [47].

First papers on AD noise in micro/nanodevices were focused on resonant mechanical structures [12, 48, 49]. Dating from the same period is the first paper on AD fluctuations in micro-biosensors [13]. In subsequent publications, the analysis of AD fluctuations was performed for various types of sensors (mass sensors with micro/nanocantilevers [15, 16, 50], semiconductor resistive gas sensors [51], plasmonic sensors [52], quartz crystal microbalance gas sensors [19]), which operate in a single-gas environment, assuming Langmuir adsorption. In [53] a theoretical AD noise model is presented for the Wolkenstein adsorption of particles of a single gas, applicable in the case of chemical adsorption in semiconductor resistive sensors.

In real situations, certain processes (e.g. various cases of non-specific adsorption, or mass transfer processes), which are coupled with the AD process of target particles, affect the sensor response, and thus influence the response fluctuations. Therefore, it is necessary for the AD noise analysis to include the coupling of multiple processes, depending on the considered practical case. In [54] the expression is derived for the PSD of adsorbed mass fluctuations due to the coupled AD processes of an arbitrary number of gases, by using the Langevin approach. The AD noise model in the case of multilayer adsorption according to the BET model is presented in [55, 56]. Analyte diffusion within the sensor chamber is considered together with the binding of particles to the surface sites, in order to analyze the fluctuations and noise figures of merit of biosensors in [14, 57] assuming the adsorption surface of infinite capacity. The coupling of mass transfer (convection and diffusion) and AD processes of one or multiple substances on the sensing surface is taken into account in the development of the AD noise model of microfluidic sensors, which is presented in [58, 59, 60]. The combined effect of the AD process, the mass transfer in the sensor chamber and the surface diffusion on the fluctuations of the number of adsorbed particles is analyzed in [61], and a good match is shown between the derived PSD of AD noise and the experimental results obtained by using a graphene gas sensor [62]. The influence of the analyte depletion from the sample on the AD noise is modeled and analyzed in [63]. In [64] the analysis is presented of the signal-to-noise ratio of a nanowire biosensor, based on stochastic simulations of the AD process coupled with diffusion.

According to the current trends in micro- and nanosensor development (the decrease of the sensing surface area, and the decrease of detectable concentrations), the analysis of AD fluctuations becomes increasingly significant for the estimation of limiting performances of such devices, and for optimization of sensor design and experimental methods. In spite of that, the topic of AD noise is scarce in the literature, compared to the total number of papers on chemical and biological sensors. Also, the published experimental results pertinent to AD noise are very scarce, and can be found for gas sensors [17, 18, 19, 20].

In the following part of this Section, mathematical models will be presented that take into account the existence of different types of adsorption sites on the sensing surface.

### 3.1 Langevin method for multicomponent stochastic processes

There are two approaches that are commonly used for the analysis of the fluctuations of the number of adsorbed particles: the first is based on the master equation, and the second on the Langevin equation, with the use of Wiener-Khinchin theorem. Both the approaches are described in detail in Supplementary data of Ref. [59] for a competitive AD process of two analytes, and all the given expressions can be generalized in a simple manner in order to be valid for * n*-component AD processes, where

*≥ 2. Here, the Langevin method will be presented, which enables efficient determination of the analytical expression for the power spectral density (PSD) of the sensor response fluctuation, starting from the kinetic macroscopic equations that describe a multicomponent stochastic process.*n

When adsorption and desorption processes occur on the heterogeneous active surface of a sensor, the number of particles adsorbed on each type of adsorption sites stochastically fluctuates in time. The fluctuations of the number of particles adsorbed on sites belonging to the ^{th} type are denoted with Δ* N*(

_{i}

*= 1, 2 …*i

*). All these fluctuating processes, when observed together, constitute a single complex random process, which belongs to the class of multicomponent Markov “gain and loss” processes [65]. Fluctuations on one type of sites (Δ*n

*) are a single component of that process. A complex process that has*N

_{i}

*(*n

*≥ 2) components can be represented by an*n

*-dimensional column vector*n

**= [Δ**ΔN

N

_{1}Δ

N

_{2}… Δ

*]*N

_{n}

^{T}(the superscript “T” denotes matrix transposition). The Langevin equation expressed in the matrix form is then (Supplementary data of Ref. [59])

where ** ξ**(

*) = [*t

ξ

_{1}(

*)*t

ξ

_{2}(

*)…*t

*(*ξ

_{n}

*)]*t

^{T}is the vector of the Langevin source functions, and

**is the square**K

*×*n

*matrix whose elements are determined by the process parameters that influence the dynamics of analyte particles binding-unbinding random events. Elements of this matrix will later be derived for each of the analyzed cases of adsorption on heterogeneous surfaces. The symbol “·” is the matrix multiplication operator.*n

The Langevin equation in the complex domain (* ω* = 2π

*,*f

*being the Fourier frequency) yields*f

(where Δ** N**(j

*) = [Δ*ω

N

_{1}(j

*) Δ*ω

N

_{2}(j

*) … Δ*ω

*(j*N

_{n}

*)]*ω

^{T},

**(j**ξ

*) = [*ω

ξ

_{1}(j

*)*ω

ξ

_{2}(j

*)…*ω

*(j*ξ

_{n}

*)]*ω

^{T}, and

**is the**I

*×*n

*unity matrix), and then the*n

*×*n

*matrix of single-sided power spectral and cross-spectral densities of the numbers of adsorbed particles,*n

S

^{2}

_{ΔN}(

*), whose elements are determined by the expression*ω

The expressions for the elements of the matrix S_{ξ}^{2} can be derived by using the formal statistical approach, as presented in [66]. Here only the final result is shown (* d*are the effective probabilities of increase or decrease of the number of adsorbed particles at sites of the

_{ie}

i

^{th}type in unit time)

By using the Wiener-Khinchin theorem [67] and Eqs. (16) and (19) the expression is obtained for the PSD of the sensor response fluctuations [59], i.e. for the PSD of AD noise

where the row vector of weight factors ** W** = [

w

_{1}

w

_{2}…

w

_{n}] is introduced.

### 3.2 Adsorption limited binding

In the case of the transport-adsorption regime, known as the rapid mixing regime, the transfer flux toward the adsorption sites dominates over the adsorption flux, so the dependence of the concentration on the spatial coordinates is negligible, and it results in adsorption-limited response kinetics, determined by the kinetic Eqs. (6). These equations are mutually independent, and show the change of the number of adsorbed particles on any of * n*sites, in the form of the difference of the actual adsorption and desorption rates,

*and*a

_{i}

*, which are linear functions of the number of adsorbed particles*d

_{i}

*(*N

_{i}

*= 1,2 …*i

*). The fluctuations can then be analyzed by directly applying the Langevin method. This approach is usually applicable in the case of AD process of gas particles. In other cases (e.g. in biological sensors, where the analyte particles are typically large macromolecules with slow diffusion in liquid samples) it is necessary to take into account the influence of transfer processes on the fluctuations of the number of adsorbed particles.*n

Eqs. (6) and (15) directly yield the system of * n*independent Langevin equations

after a random source function is added on the right side of each equation.

When this system of equations is written in the matrix form, as shown in SubSection 3.1, the quadratic diagonal matrix ** K**is obtained, and its elements are

The elements of the matrix ^{2}_{ξ}are (Eq. (20))

where * N*is given by Eq. (7).

_{ie}

Since all the quadratic matrices on the right side of Eq. (19) are diagonal, the matrix ^{2}_{ΔN}(* ω*) is also diagonal, i.e. the cross-spectral densities

S

^{2}

_{ΔN,il}are equal to zero (the random processes Δ

*and Δ*N

_{i}

*are statistically independent for every pair of*N

_{l}

*and*i

*,*l

*≠*i

*). The spectral densitities of fluctuations of the numbers of adsorbed particles are of the Lorenzian type*l

with the characteristic frequency _{c,AL,i} = 1/(2π_{AL,i}), where _{AL,i} = * k* +

_{ai}C

_{0}

*, and the AD noise PSD of a sensor with*k

_{di}

*types of adsorption sites is determined by the sum of Lorenzians*n

### 3.3 Mass transfer influenced binding

The use of TCM for approximation of spatially and temporally dependent analyte concentration in a sensor chamber (see SubSection 2.3) has enabled the presentation of the macroscopic kinetic equations (Eqs. (10), (11) and (14)) in the form in which the change of the number of adsorbed particles in unit time equals the difference between the instantaneous effective rates of adsorption and desorption, _{eff,i}and _{eff,i}, which are explicitely dependent only on the instantaneous numbers of adsorbed particles on adsorption sites of all types. Such form of the equations is suitable for fluctuation analysis. In order for the Langevin method to be applicable for obtaining the PSD of fluctuations of a sensor response when coupling of complex (* n*-component) AD process and mass transfer is considered, a linear approximation of the functions

a

_{eff,i}and

d

_{eff,i}around the equilibrium values of the numbers of adsorbed particles is used (assuming small fluctuations relative to the equilibrium values, Δ

*<<*N

_{i}

*) [59, 66]*N

_{ie}

All the derivatives are calculated for _{1} = _{1e}, _{2} = _{2e}… * N* =

_{n}

*(the equilibrium values are determined by Eq. (12)), which is in the above expressions denoted with the subscript “e” within the derivatives. It can be noticed that*N

_{ne}

a

_{eff,i}(

N

_{1e},

N

_{2e}…

*) =*N

_{ne}

d

_{eff,i}(

*), since the adsorption and desorption rates are pertinent to the equilibrium state, according to the condition d*N

_{ie}

*/d*N

_{i}

*= 0. By substituting Eqs. (27) and (28) into Eqs. (10), linearized kinetic equations are obtained, which, after the addition of Langevin source functions on their right side assume the form of Langevin equations, shown by the matrix Eq. (17). The obtained expressions imply that the elements of the quadratic matrix*t

**are**K

and after differentiation

Based on Eq. (20)

Now, all the quantities are known that enable the determination of the spectral and cross-spectral densities of fluctuations of the numbers of adsorbed particles on different types of adsorption sites, according to Eq. (19), as well as the spectral density of the sensor response fluctuations, based on Eq. (21)

When there are two types of adsorption sites on the sensing surface, the elements of the matrices ** K**and

S

^{2}

_{ξ}are, according to Eqs. (30)–(32)

so the PSD of the AD noise, calculated based on Eq. (21), is

where

The characteristic frequencies of the AD noise spectrum are _{c,MT,i} = 1/(2π_{MT,i}), where * i*is 1, 2 or 3.

The expressions for * K*and

_{ii}

*for a fast mass transfer (at suffuciently high*K

_{ij}

*) become approximately equal to the expressions given by Eq. (23), which are valid for the case of adsorption-limited binding. Also, in that case, Eq. (33) reduces to Eq. (26) for*k

_{m}

*= 2. This demonstrates a wider applicability of the expressions derived by using TCM. Namely, when the zone adjacent to the sensing surface, which contains the analyte particles that can participate in the adsorption is narrow, the expressions derived by using TCM for mass transfer influenced binding are valid even in the case of high*n

*values, when the binding is adsorption limited.*k

_{m}

## 4. Results and discussion

In order to analyze the effects of mass transfer and surface heterogeneity on the sensor temporal response and AD noise, numerical calculations are performed for the case of protein biosensor, on whose active surface two types of adsorption sites exist with different affinities toward the target analyte. The results are presented in terms of the adsorbed mass on the sensing surface, assuming the values of adsorption and desorption rate constants from the ranges corresponding to biomolecules and biosensors [68]: _{a1} = 1.3^{.}10^{−18} m^{3}/s, _{d1} = 0.4 1/s, _{a2} = 1.3^{.}10^{−20} m^{3}/s and _{d2} = 0.02 1/s, the mass of a single analyte particle * w* = 20 kDa, and the total of

*= 10*N

_{max}

^{8}adsorption sites on the sensing surface of area

*= 10*A

^{−9}m

^{2}. The analyte concentration is

C

_{0}= 5

^{.}10

^{17}1/m

^{3}.

Figure 2 shows the temporal change of adsorbed mass, representing the temporal response of the sensor whose measured parameter is a function of the total bound mass of the target protein. The total mass adsorbed on the sensing surface is shown (black lines), as well as the adsorbed amounts on each of the two types of adsorption sites (red lines for the type 1 sites, blue lines for the type 2 sites). The five diagrams correspond to different shares of the types of adsorption sites in the overall number of sites, which is expressed by different values of the parameter * ν*(1; 0.8; 0.5; 0.2, 0), where

N

_{max,1}=

*and*νN

_{max}

N

_{max,2}= (1-

*)*ν

*. Different cases are illustrated: from the presence of only the high-affinity sites (type 1 sites), to various ratios of the numbers of sites belonging to the two types (*N

_{max}

N

_{max,1}:

N

_{max,2}that equals 4:1, 1:1 and 1:4), and, finally, to the presence of only the low-affinity sites (type 2 sites) on the sensing surface. As the measure of the affinity, the affinity constant is used, which is defined by the ratio

*/*k

_{ai}

*. The curves are obtained by computer simulation, based on the model that takes into account the mass transfer effects (Eqs. (14)), for two values of the mass transfer coefficient:*k

_{di}

k

_{mI}= 2

^{.}10

^{−3}m/s (solid line) and

k

_{mII}= 9

^{.}10

^{−1}m/s (dashed line). For

*>*k

_{m}

k

_{mII}, the obtained curves overlap with those for

k

_{mII}, and they also overlap with the curves obtained by the use of the model that neglects mass transfer (Eqs. (6) and (7)). This means that the value of

k

_{mII}is high enough for the influence of mass transfer on the sensor’s temporal response to be considered as negligible, so that conclusions about the mass transfer influence can be made by comparing the responses for

k

_{mI}and

k

_{mII}.

In the absence of the mass transfer influence (dashed-line curves shown in the diagram), the existence of two types of adsorption sites is clearly noticeable based on the sharp transition from fast to slow transient regime (for * ν*equal to 0.8, 0.5 or 0.2). The time evolution of the total adsorbed mass is rapid at first, as it is dominantly determined by adsorption on the higher affinity sites, with the time constant of approximately 1 s, but it then slowly approaches the steady state, with the time constant of approximately 38 s (AD process on sites of lower affinity). Even with a lower share of low affinity sites, the AD process occurring on them determines the sensor response rate. Different shares of the two types of adsorption sites influence the sensor response value in the steady state.

When the mass transfer is characterized by _{mI}, a pronounced MT influence on the sensor response kinetics is evident at all values of * ν*. Such results are in accord with Eq. (13) by which the analytic criterion is defined for a negligible influence of MT on the response. For two types of adsorption sites, that criterion is

*>>*k

_{m}

*, where*k

_{mg}

*= (*k

_{mg}

k

_{a1}·

*+*ν

k

_{a2}·(1-

*))·*ν

*/*N

_{max}

*, and for the given parameters*A

*equals 1.3*k

_{mg}

^{.}10

^{−1}m/s, 1

^{.}10

^{−1}m/s, 6.6

^{.}10

^{−2}m/s, 2.7

^{.}10

^{−2}m/s, and 1.3

^{.}10

^{−3}m/s for

*equal to 1, 0.8, 0.5, 0.2 and 0, respectively. Hence, when*ν

*=*k

_{m}

k

_{mI}, the above-mentioned condition is not satisfied for any

*. It can be seen that the minimal MT coefficient value at which the response kinetics can be considered as adsorption limited increases with the increase of the type 1 sites share. The influence of mass transfer characterized by the parameter*ν

k

_{mI}on the time evolution of the number of adsorbed particles is more pronounced for the sites with a higher adsorption rate constant. Thus, the initial evolution of the total number of adsorbed particles is much slower compared to the case of adsorption limited binding.

Figure 3 shows the spectral density (SD) of adsorbed mass fluctuations, (_{ΔR}^{2}(* f*))

^{1/2}, i.e. the spectral density of sensor AD noise, for the same parameter values for which the response shown in Figure 2 was obtained. It is obtained as a square root of the PSD, given by Eq. (33), which, apart from surface heterogeneity, takes into account mass transfer. SDs of the total adsorbed mass are shown for two values of the MT coefficient:

k

_{mI}= 2

^{.}10

^{−3}m/s (solid red lines) and

k

_{mII}= 9

^{.}10

^{−1}m/s (solid blue lines), for each of five values of

*. As in the case of response, the curves obtained for*ν

*>*k

_{m}

k

_{mII}match those for

k

_{mII}. Also matching with them are the curves obtained by using Eq. (26), i.e. the model that neglects mass transfer. According to that model, PSD of AD noise equals the sum of PSDs of fluctuations of adsorbed masses on two types of sites. The SDs corresponding to these components are also shown in the diagrams (dotted lines denote type 1 sites, and dashed lines denote type 2 sites). Hence, the influence of mass transfer on the AD noise is negligible for

*=*k

_{m}

k

_{mII}, and the comparison of the results for the shown SDs of total fluctuations for a given

*enables making conclusions on the influence of mass transfer on the sensor AD noise.*ν

It can be noticed that at all the values of * ν*, mass transfer leads to the increase of the AD noise magnitude in the part of the spectrum where the characteristic frequencies belong. It also increases the AD noise (obtained by integration of the AD noise SD in the frequency range of interest), and thus influences the fundamental detection and quantification limits of a biosensor. For a given

*, mass transfer also causes a shift of the characteristic frequencies of the AD noise spectrum. Thus, the AD noise spectral analysis can yield information not only on adsorption and desorption rate constants and analyte concentration, but also on mass transfer parameters.*ν

By comparing the diagrams for the case of a surface containing a single type of adsorption sites (i.e. for * ν* = 1 or

*= 0) with the diagrams for the sensors with a heterogeneous adsorbing surface (i.e. when*ν

*equals 0.8, 0.5 or 0.2), it is obvious that the latter diagrams may exhibit a greater number of characteristic frequencies (three in the case of two type of adsorption sites), which can indicate the existence of different types of adsorption sites.*ν

The diagrams also show that when the mass transfer influence is negligible, the two characteristic frequencies of the AD spectrum do not change at different shares of adsorption site types on the sensing surface. They are determined by the characteristic frequencies that correspond to the fluctuation spectra of adsorbed masses on each of site types, which are given by expressions _{c,AL,i} = 1/(2π_{AL,i}), where _{AL,i}= * k* +

_{ai}C

_{0}

*, as given in Section 3.2. Contrary to that, when the mass transfer influence is pronounced, all the three AD noise spectrum characteristic frequencies change with*k

_{di}

*, so that their values also contain the data on the share of a specific site type on the sensing surface. Therefore, as the parameters of all the processes and effects that influence the sensor temporal response are contained in the characteristic features of the AD noise spectrum, noise spectrum analysis can be used as an additional source of data in biosensing.*ν

## 5. Conclusions

In this chapter, the mathematical modeling has been performed of time response and adsorption–desorption (AD) noise in microfluidic adsorption-based biosensors whose active surface is heterogeneous in the sense that it contains adsorption sites of different affinity toward the target analyte. The adsorption–desorption processes of analyte particles on different types of adsorption sites, as well as mass transfer in a microfluidic chamber, have been taken into account as the coupled processes that generate the sensor signal. The devised model of AD noise is the first that simultaneously takes into account surface heterogeneity and mass transfer through both the convection and diffusion of analyte particles. The analytical expression of the spectral density of AD noise has been derived. The models of the time response and AD noise of a sensor with heterogeneous adsorbing surface in the case of negligible mass transfer influence have also been presented. The criterion is given based on which it can be discerned whether the mass transfer influence is significant or negligible, so that the appropriate mathematical model of the response and noise can be chosen.

The derived mathematical models have been used for the analysis of the response and AD noise of protein biosensors that have two types of adsorption sites on the sensing surface. The comparison of the results obtained by the use of the models that take into account mass transfer effects, and those that neglect them, for different shares of two types of binding sites, has enabled drawing conclusions on both separate and coupled influences of surface heterogeneity and transport processes on the sensor response and noise.

While a slow mass transfer increases the sensor response time, the existence of different adsorption site types affects both the transient regime and the sensor response magnitude in the steady state. The minimal mass transfer coefficient value at which the response kinetics can be considered as adsorption limited increases with the increase of the share of high affinity binding sites. The influence of mass transfer (of given parameters) on the time evolution of the number of adsorbed particles is more pronounced on sites with a higher adsorption rate constant.

At any share of the two types of binding sites on the sensing surface, mass transfer causes the increase of the AD noise, and thus increases the sensor’s fundamental limits of analyte detection and quantification. For a given ratio of the numbers of adsorption sites of the two types, mass transfer causes a shift of the characteristic frequencies in the AD noise spectrum. When the mass transfer influence is pronounced, the characteristic frequencies shift with the change of the ratio of the numbers of sites of different types. Therefore, the AD noise spectrum analysis can yield information not only on adsorption and desorption rate constants and analyte concentration, but also on mass transfer parameters and on the share of a certain type of binding sites on the sensing surface. It can be used as an additional source of data in biosensing.

As the results of the performed analysis have shown a potentially significant influence of sensing surface heterogeneity and mass transfer processes on the sensor temporal response, as well as on AD noise, which is inevitable in adsorption-based sensors and determines their ultimate sensing performance, the presented mathematical models can enable better interpretation of measurement results, and give guidelines for ensuring lower noise levels and improved detection limits in microfluidic biosensors.

## Acknowledgments

This research was funded by the Ministry of Education, Science and Technological Development of the Republic of Serbia, grant number 451-03-9/2021-2114/200026. It was also supported by the Science Fund of the Republic of Serbia, PROMIS, #6057070, Gramulsen.