Universal optimal geometry of minimal phoretic pumps

TL;DR

This paper reveals that the optimal design for minimal phoretic pumps, which use chemically active patches to generate flow, is universal and does not depend on specific chemical properties.

Contribution

It identifies the universal optimal arrangement of three chemical patches in minimal phoretic pumps, advancing understanding of their design principles.

Findings

Optimal patch arrangement maximizes flow rate.

Universal design independent of chemical specifics.

Provides guidelines for efficient phoretic pump construction.

Abstract

Unlike pressure-driven flows, surface-mediated phoretic flows provide efficient means to drive fluid motion on very small scales. Colloidal particles covered with chemically-active patches with nonzero phoretic mobility (e.g. Janus particles) swim using self-generated gradients, and similar physics can be exploited to create phoretic pumps. Here we analyse in detail the design principles of phoretic pumps and show that for a minimal phoretic pump, consisting of 3 distinct chemical patches, the optimal arrangement of the patches maximizing the flow rate is universal and independent of chemistry.

Equations

$${D\frac{\partial C}{\partial y}|}_{y=0}=-A(x),$$

$$C(x,y)=C_0-\frac{L}{2\pi D}\sum _{n=-\infty }^{\infty }\frac{a_n}{n}\frac{\sinh \left[\frac{2n\pi (y-H)}{L}\right]{\mathrm{e}}^{\frac{2\mathrm{i}\pi nx}{L}}}{\cosh \left[\frac{2n\pi H}{L}\right]},$$

$$a_n=\frac{1}{L}{\int }_0^{L}A(x)\exp \left(-\frac{2\mathrm{i}n\pi x}{L}\right).$$

$${(\mathbf{u}\cdot {\mathbf{e}}_x)|}_{y=0}={M(x)\frac{\partial C}{\partial x}|}_{y=0},$$

$$Q={\int }_{S_x}\mathbf{u}\cdot \mathrm{d}\mathbf{S}=\frac{H}{2}⟨u_s(x)⟩,$$

$$Q=-\frac{H}{D}\sum _{n=1}^{\infty }\tanh \left(\frac{2n\pi H}{L}\right)\text{mbox}Im[a_n{m}_{-n}],$$

$$A(x)=\frac{A_2{M}_1-A_1{M}_2}{M_1-M_2}+\left(\frac{A_1-A_2}{M_1-M_2}\right)M(x).$$

$$A(x)=\sum _{p=1}^P{A}_p{1}_{S_p}(x),M(x)=\sum _{p=1}^P{M}_p{1}_{S_p}(x),$$

$$Q/L=\sum _{n=1}^{\infty }$$

$$Q/L$$

$$\mathcal{F}(A,M)=\sum _{j=1}^3(A_{j+1}-A_j)(M_j+M_{j+1}),$$

$${\left(\frac{\partial \mathcal{G}}{\partial l_1}\right)}_{l_2}$$

$$C=C_0+\frac{(A_j+A_{j+1})(H-y)}{2}+\frac{(A_j-A_{j+1})}{2}\tilde{C}\left(\frac{x-x_j}{H}\right),$$

$$Q_{j,j+1}$$

$$Q/L\sim \frac{h^2}{4D}\sum _{j=1}^P(A_{j+1}-A_j)(M_j+M_{j+1})=\frac{h^2}{4}\sum _{j=1}^P{\alpha }_{j,j+1}.$$

$$U$$

$$Q/L=\sum _{j=1}^4{\mathcal{F}}_j(A,M){\mathcal{H}}_j(l_1,l_2,l_3,l_4,h),\text{with }$$

$${\mathcal{F}}_1={\alpha }_{23}+{\alpha }_{34}+{\alpha }_{42},{\mathcal{H}}_1=\sum _{n=1}^{\infty }\frac{h(-1{)}^{n+1}\tanh (2\pi nh)}{{\pi }^2{n}^2}\cos (\pi n{l}_1)\prod _{k=2}^4\sin (\pi n{l}_k).$$

Full text

I Introduction

The rapid development of microfluidics, which already has had a deep impact on both biology and chemistry whitesides01 ; beebe02 ; hansen03 ; whitesides2006 , was enabled by key advances in continuum physics. Indeed, it is our understanding of the surface-dominated physics at the micron scale which has allowed the invention of a whole array of small-scale devices to precisely control flow and transport processes in microfluidic devices squires2005 .

One of the standard issues in small devices is the difficulty of driving flows. In a straight channel, or pipe, the rate at which a Newtonian liquid flows from one side of the channel to the other scales as the fourth power of a relevant cross-sectional channel length scale times the applied external pressure gradient stonereview . When length scales become tens of microns or less, the external pressures required to drive flows become prohibitively large and as a result the community has turned to surface-driven methods where a flow is induced locally squires2005 ; stonereview .

In the biological world, surface flows are often created along tissues, or groups of cells, by the time-varying beating of short cilia sleigh1988 resulting in effective slip boundary conditions for the neighbouring flow Blake1971a ; brennen1977 . Although artificial cilia have been realised in the lab, the dynamics and performance of biological ciliary arrays has proven difficult to reproduce experimentally fahrni2009 ; babataheri2011 ; coq2011 ,

Instead, a popular method to generate flows near surfaces in the lab consists in taking advantage of phoretic mechanisms where externally-applied physico-chemical gradients (such as charge, temperature, composition…) create local body forces on the fluid in thin layers near surfaces which, through the action of viscous stresses, entrain a bulk flow anderson1989 . A famous example of such methods is electrophoresis wherein an electric field applied along a channel filled with an electrolyte drives a flow due to charge imbalance near the electrical double layer at the junction between the fluid and surfaces squires2005 .

While externally-applied gradients are able to drive flows, gradients which are instead generated locally on the surface of colloidal particles can be used to generate locomotion Howse2007 ; Ebbens2011 ; julicher2009 ; moran2017 . Self-propulsion of such phoretic swimmers can result either from chemical gradients directly patterned on the particles themselves via coated catalysts golestanian2007 or from transport instabilities for chemically-homogenous particles thutupalli2011 ; michelin2013c ; izri2014 , and have proven popular model systems in the field of active matter marchetti_review . A canonical example of such catalytic reactions is the decomposition of hydrogene peroxide on platinuum-coated surfaces Howse2007 or iron oxide catalysts palacci2013 , but many other chemical reactions have also been considered duan2015 ; yadav2015 .

The physico-chemical principles used for phoretic swimmers can in principle also be exploited to induce flow transport in confined devices such as microchannels, and therefore to create pumps michelin2015b ; shen2016 ; yang2016 ; tan2017 . Yet, the existing literature has only so far provided limited insight on the fundamental design principles of such pumps, and we propose here a detailed analysis of the link between pump design and performance. In particular, with a view toward experimental realisation, an important practical question is that of minimal geometrical design. What type of surface chemistry would be simple to fabricate yet effective at creating transport?

In the case of swimmers, the minimal design is that of a Janus particle whose surface is covered by two distinct, homogeneous patches of which at least one is phoretically active. By symmetry, a Janus channel cannot be used to pump flows, and the simplest design has three patches. In this paper, we solve theoretically the $P$-patch problem. We demonstrate that in the minimal case of $P=3$ patches, the optimal pump design, i.e. the geometrical arrangement of chemical patches leading to maximum phoretic flow rate, is universal and independent of chemistry, in stark contrast with Janus particles.

II Results

II.1 Model and performance of a generic phoretic pump

We consider an infinite, straight two-dimensional channel of width $H$ (Fig. 1). One of the channel walls, located at $y=0$, is chemically-coated with a catalyst along a repeated pattern of period $L$. The catalyst allows a chemical reactant in the liquid to produce a new solute species of concentration $C(\mathbf{x})$. In the limit of large reactant concentration, we may assume that the solute release occurs at a fixed rate, or activity, $A(x)$,

[TABLE]

where $D$ is the molecular diffusivity of the solute. At sufficiently small length scales, both advective and unsteady transports are negligible and the dynamics of the solute concentration is purely diffusive, $D\nabla^{2}C=0$. For simplicity, we assume that the upper wall allows for free exchanges of solute with a chemical reservoir so that the product concentration along it is homogeneous, $C(y=H)=C_{0}$. These conditions uniquely determine the solute concentration within the channel as

[TABLE]

with $a_{n}$ the Fourier coefficients of $A(x)$ given by

[TABLE]

Due to the differential affinity of the chemically-patterned wall with the reactant and product molecules, local surface gradients in solute concentration result in a net slip velocity outside a thin interaction layer providing an effective slip boundary condition for the flow velocity, $\mathbf{u}$, as anderson1989

[TABLE]

with $M(x)$ the local diffusiophoretic mobility along the wall. This mobility stems from the difference in affinity with the wall surface (or short-range interaction potential) between the solute and solvent molecules within a thin interaction layer anderson1989 . This simple framework can be easily generalized to other phoretic mechanism, such as thermophoresis anderson1989 ; bickel2013 ; baraban2013 , or other geometries (e.g. axisymmetric channels or patterning of top and bottom walls).

Using the fundamental properties of Stokes’ flow, the resulting flow rate induced by the phoretic pump through any cross-section $\mathcal{S}_{x}$ of the channel is given by michelin2015c

[TABLE]

where $\langle\cdot\rangle$ is the average in $x$ over a period $[0,L]$. Using Eqs (2), (4) and (5), the pumping rate $Q$ is then obtained as

[TABLE]

where $m_{n}$ are the Fourier coefficients of the mobility $M(x)$.

As expected for phoretic problems in the diffusive limit, the pumping rate is a bilinear function of the activity, $A(x)$, and mobility, $M(x)$, and no pumping is possible if either is constant along the channel, nor if $M(x)=\lambda A(x)+\mu$ where $\lambda$ and $\mu$ are two arbitrary constants. For a Janus-type channel patterning consisting of the repetition of two patches with properties $(A_{1},M_{1})$ and $(A_{2},M_{2})$, one can write

[TABLE]

Consequently, two-patch patterns are unable to pump, a result which was expected since such systems are left-right symmetric (i.e. $x\leftrightarrow-x$) with respect to the midpoint of any of the patches. This is of course a fundamental difference with phoretic propulsion of microparticles for which a minimal two-patched Janus patterning leads in general to locomotion golestanian2007 .

II.2 Pumping rate of a $P$-patch channel

While realizing continuous variations of the chemical properties of the wall is experimentally difficult, a simple patterning consists of the periodic repetition of $P\geq 3$ patches: on each patch $S_{p}$ of length $L_{p}$ (with $\sum_{j=1}^{P}L_{j}=L$), both $A(x)$ and $M(x)$ are constant and take values $A_{p}$ and $M_{p}$, i.e.

[TABLE]

where $\mathbf{1}_{S_{p}}(x)=1$ for $x\in S_{p}$ and {\color[rgb]{0,0,0}\mathbf{1}_{S_{p}}(x)}=0 otherwise. The Fourier coefficients $a_{n}$ and $m_{n}$ can be obtained from Eq. (3), and the flow rate of the channel is computed from Eq. (6) as

[TABLE]

with $l_{p}=L_{p}/L$ the reduced length of $S_{p}$, $h=H/L$ the channel aspect ratio and $\alpha_{pq}=(M_{p}A_{q}-M_{q}A_{p})/D$. This generic form expresses the pumping rate in the channel as the sum of pair interactions between patches, whose intensity depends on their lengths and the distance between their centers. Note that the flow rate $Q$ depends on the $P(P-1)/2$ coefficients $\alpha_{pq}=-\alpha_{qp}$ rather than the $2P$ chemical characteristics $(A_{j},M_{j})$ but $\alpha_{pq}$ may not be defined independently from each other. These coefficients also set the characteristic velocity scales generated in such pumps, which are similar to those for the flows generated by phoretic swimmers golestanian2007 .

II.3 The optimal and minimal phoretic pump

Since channels with $P=2$ can never pump, the minimal phoretic pump has $P=3$ patches. In that case, using $\sum l_{j}=1$, the pumping rate in Eq. (9) becomes SM

[TABLE]

and is the product of two functions: (i) $\mathcal{F}(A,M)=\alpha_{12}+\alpha_{23}+\alpha_{31}$ which depends exclusively on the chemical properties of the patches and (ii) $\mathcal{G}(l_{1},l_{2},l_{3},h)$ which depends only on the geometry of both channel and patches. Note that the function $\mathcal{G}$, written here in a symmetric form with respect to $(l_{i})_{1\leq i\leq 3}$, is effectively a function of $l_{1}$, $l_{2}$ and $h$ only (since $l_{1}+l_{2}+l_{3}=1$).

The chemical function $\mathcal{F}$ can be rewritten (using the convention $A_{j+3}=A_{j}$),

[TABLE]

and can thus be interpreted as the sum of contributions of adjacent pairs of patches which each induces a net flow proportional to the mean mobility multiplied by the difference in activity. A similar result is at the heart of the self-propulsion of Janus microswimmers golestanian2007 .

The explicit separation of the chemical and geometric dependences of the pumping rate in Eq. (10) confers a universality to the three-patch configuration: The variation of the flow rate with the geometric patterning of the channel is not affected by the values of the chemical activities or mobilities. In particular, this means that the optimal pump, found by maximising the function $\mathcal{G}$, is unique and identical for all chemistry.

The variation of $\mathcal{G}$ within the $2D$ parameter space $\mathcal{I}_{3}=\left\{0\leq l_{1},l_{2},l_{3}\leq 1,\,\sum l_{i}=1\right\}$ is shown in Fig. 2. For all aspect ratios $h$, $\mathcal{G}$ vanishes if any $l_{j}=0$ (Janus limit), which are the boundary points on $\mathcal{I}_{3}$. The gradient of $\mathcal{G}$ with respect to $(l_{1},l_{2})$ is given by

[TABLE]

with $l_{3}=1-l_{1}-l_{2}$ and $(\partial{\cal G}/\partial l_{2})$ is obtained similarly; the only point within $\mathcal{I}_{3}$ where $|\cal G|$ has a maximum is $l_{1}=l_{2}=l_{3}=1/3$, which confirms the results of Fig. 2. The optimal minimal (3-patch) pump is therefore unique and, independently of the chemistry, is the one where all patches have equal lengths.

The dependence of the pumping ability of the channel on its geometry can be further understood by examining the impact of the channel aspect ratio, $h=H/L$, on the optimal flow rate, ${\cal G}_{\rm max}$ (Fig. 3). For large $h$, the flow rate varies linearly with $h$. In that case, the concentration distribution at the lower walls is independent of $h$ at leading order, except for its mean value which does not contribute to pumping. The resulting phoretic slip forcing is therefore independent of $h$, and similarly to Couette (shear) flow, the net pumping is linear in $h$.

In the opposite limit, $h\ll 1$, the maximum flow rate scales quadratically with $h$ (Fig. 3). In that case, the concentration profile is almost piecewise constant in $x$ (i.e. away from the junctions between patches). Zooming in on the boundary between patches $j$ and $j+1$ for $x\approx x_{j}$, the leading-order concentration $c$ can be rewritten as

[TABLE]

where $\tilde{C}(s)$ is an odd function of $s$ with $\tilde{C}(\pm\infty)\pm 1$. Since $M(x)$ is piecewise constant near $x_{j}$, the resulting contribution of this junction to the pumping flow rate, $Q_{j,j+1}$, is obtained at leading-order for $h\ll 1$ as

[TABLE]

The total flow rate thus depends only on the junction between adjacent patches and is then obtained (for $P$ patches) as

[TABLE]

Comparing with Eq. (10) for $P=3$ shows that $\mathcal{G}_{\textrm{max}}=h^{2}/4$, in excellent agreement with the full solution (Fig. 3), and that $\mathcal{G}\approx\mathcal{G}_{\textrm{max}}$ when $h\ll 1$, as also observed on Fig. 2, demonstrating the robustness of the optimal design in that limit.

The universal nature of the optimal geometry for minimal ($P=3$) phoretic pumps is intimately linked to the number of independent chemical properties setting the flow rate. For $P$ patches, $2P$ different properties come into play, ($A_{i},M_{i}$). Denoting by $\mathcal{A}$ and $\mathcal{M}$ a characteristic magnitude of activity and mobility, dimensional analysis imposes that $Q=\mathcal{AM}\times\tilde{Q}$, and $\tilde{Q}$ depends on only $2(P-1)$ parameters. No net pumping is obtained if either $(A_{i})_{i}$ or $(M_{i})_{i}$ are all identical, or if both sets are linearly correlated, providing three additional constraints, such that the pumping rate effectively only depends on $2P-5$ independent chemical parameters. For $P=3$ patches, this confirms that a single chemical function controls the pumping rate, conferring its universality to the minimal pump.

II.4 Minimal phoretic swimmers vs. minimal phoretic pumps

While the minimal phoretic pump must include three different patches, minimal phoretic swimmers are able to break symmetries using only two. However, in contrast to the results obtained above and showing universality of the three-patch pump, the optimal minimal (Janus) swimmer is not universal but its geometry depends on the surface chemistry. This can be seen by evaluating the swimming velocity of an unit-radius axisymmetric Janus sphere coated with two different materials $(A_{1},M_{1})$ on the portion $\mu\leq z\leq 1$ of its surface (front side) and $(A_{2},M_{2})$ for $-1\leq z\leq\mu$ (back). The result is

[TABLE]

where $L_{n}^{\prime}(\mu)$ is the derivative of the $n$-th Legendre polynomial golestanian2007 . The non-universality of Janus swimmers can then be demonstrated by highlighting a few examples. When $M_{1}=M_{2}$, the Janus swimmer with maximum speed is hemispheric, and thus the optimal value is $\mu_{\textrm{opt}}=0$. In contrast, when $M_{1}=-M_{2}$, the hemispheric particle with $\mu=0$ does not swim and instead maximizes its swimming speed for $\mu_{\textrm{opt}}\approx\pm 0.61$. The optimal Janus swimmer, i.e. the value of $\mu$ maximizing $|U|$, is therefore not universal and optimizing the patterning of the surface of the swimmer requires a detailed knowledge of its chemical properties, in contrast with minimal phoretic pumps which are always optimal for $l_{i}=1/3$.

II.5 Optimal pumps beyond $3$-patch patterns

The universality for pumps is lost for $P>3$ as the pumping rate now depends on $2P-5>1$ independent chemical parameters. In the case of $P=4$ patches, the net flow rate, Eq. (9), becomes

[TABLE]

with $\mathcal{F}_{j}$ and $\mathcal{H}_{j}$ obtained by circular permutation for $j\geq 2$ SM . The $\mathcal{F}_{j}$ contribution is essentially a modulation of the $3$-patch pump obtained for $l_{j}=0$. The pumping rate nevertheless depends on only three independent parameters since these four contributions are not independent ($\mathcal{F}_{1}+\mathcal{F}_{3}=\mathcal{F}_{2}+\mathcal{F}_{4}$). All possible geometries now span the three-dimensional parameter space $\mathcal{I}_{4}=\left\{0\leq l_{1},l_{2},l_{3},l_{4}\leq 1,\,\sum l_{i}=1\right\}$. Depending on the surface chemistry, the optimal pumping rate is reached (i) within $\mathcal{I}_{4}$ if $\mathcal{F}_{1}\mathcal{F}_{3}$ and $\mathcal{F}_{2}\mathcal{F}_{4}$ are both positive (non-trivial $4$-patch pump) or (ii) on its boundary if either quantity is negative (degenerated $3$-patch pump), in which case the universal optimal pump with three equal-length patches is recovered SM . These two possibilities are illustrated on Fig. 4 where the dependence of the pumping efficiency on the geometry of patterning is represented over $\mathcal{I}_{4}$.

III Discussion

In summary, this work proposes a generic mathematical framework to evaluate and optimize the phoretic pumping performance of a straight microchannel periodically-coated with active surfaces. Focusing on patterns which are well suited for experimental realization, namely a succession of materials with uniform chemical properties (patches), we show that the minimal pump features three different patches and is optimal for three patches of equal lengths regardless of their chemical properties. Although we focused on diffusiophoresis, our results are also applicable to thermophoresis and electrophoresis (at least in the weak gradient limit when surface slip is proportional to the concentration gradient anderson1989 ) and could be extended to more complex geometries using numerical computations. For clarity and generality, we purposely focused here on the simplest chemical formulation of the problem, i.e. a prescribed fixed-flux of a single chemical component (reactant or product). Our framework could nevertheless be extended to account for a more detailed description of the chemical reaction, for example by including several chemical components or multi-step reactions to describe the wall activity.

The most important result of our study is the universality of the optimal geometric design. This is a unique feature of the phoretic pumping problem that does not have an equivalent for phoretic swimmers. Furthermore, this universality is likely to be critical for experimental development since determining independently the chemical and phoretic properties of active materials is challenging experimentally. There is therefore no need for a trial-and-error experimental approach to phoretic pumps.

Acknowledgements.

This project has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreements 714027 to SM and 682754 to EL).