Entrance Effects in Concentration-Gradient-Driven Flow Through an Ultrathin Porous Membrane

Transport of liquid mixtures through porous membranes is central to processes such as desalination, chemical separations and energy harvesting, with ultrathin membranes made from novel 2D nanomaterials showing exceptional promise. Here we derive, for the first time, general equations for the solution and solute fluxes through a circular pore in an ultrathin planar membrane induced by a solute concentration gradient. We show that the equations accurately capture the fluid fluxes measured in finite-element numerical simulations for weak solute-membrane interactions. We also derive scaling laws for these fluxes as a function of the pore size and the strength and range of solute-membrane interactions. These scaling relationships differ markedly from those for concentration-gradient-driven flow through a long cylindrical pore or for flow induced by a pressure gradient or electric field through a pore in an ultrathin membrane. These results have broad implications for transport of liquid mixtures through membranes with a thickness on the order of the characteristic pore size.


I. INTRODUCTION
Fluid transport through pores and porous membranes plays a key role in many processes of fundamental and practical interest, including cellular homeostasis in biological systems, 1 chemical separations, 2 desalination, 3 and energy conversion. 4,5Thus, a general theoretical understanding of the parameters that control these transport phenomena has broad implications for a variety of domains.Many theoretical models of fluid transport in porous membranes have considered flows only within the pores [6][7][8][9] and have neglected the effect of transport between the membrane pores and the fluid outside the membrane.These socalled entrance or access effects can dominate fluid transport processes when the membrane thickness approaches the characteristic pore size 10,11 or when the fluid-solid friction becomes small. 12,13The most extreme examples of this situation are membranes of atomic thickness made from 2D materials such as graphene and its derivatives [14][15][16][17][18][19] or molybdenum sulfide. 20,21ch 2D membranes are of great interest, as they have been shown to exhibit exceptional properties compared with conventional membranes for applications such as desalination 14 and electrical energy harvesting from salinity gradients. 20uid fluxes across a membrane can be induced by a variety of driving forces, including gradients of pressure, electrical potential, or solute concentration.Equations have previously been derived to quantify entrance effects on fluid flow driven by a pressure gradient 22,23 and on fluid flow 24 and ionic electrical currents 25,26 induced by an electric field acting on a electrolyte solution.However, to date, no theory has been developed to describe the entrance effects on fluid fluxes driven by concentration gradients and how they vary with relevant parameters.
Fluid fluxes driven by concentration gradients are of particular relevance in applications such as chemical separations, 2,16 desalination, 3,14,15,17 and salinity-gradient-driven energy harvesting. 4,5,13,20,27The work presented here focuses specificially on entrance effects on the concentration-gradient-driven process of diffusioosmosis, 28 in which flow of a solutecontaining solution is driven by an osmotic-pressure gradient that develops in the inhomogeneous interfacial fluid layer induced by interactions of the fluid with the solid surfaces.
Diffusioosmosis has been shown to play a key role in astonishing energy densities measured for salinity gradient energy harvesting in a nanotube membrane. 27Thus, entrance effects on this phenomenon are of considerable interest.
Here we derive, for the first time, general equations to quantify the diffusioosmotic solution flux and solute flux of a dilute solution through a circular aperture in a 2D membrane as a function of the aperture size and the strength and range of the interactions between the solute and membrane surface.We verify the accuracy of the equations by comparison with finite-element numerical simulations.We go on to compare the scaling behavior predicted for concentration-gradient-driven flow through a circular aperture with those for other membrane geometries and driving forces and discuss the implications of these results for real systems.

A. Diffusioosmotic flow
Consider the flow of a solution containing a single solute type through a circular aperture of radius a in an infinitesimally thin planar wall, as illustrated in Fig. 1.Assuming that the fluid flows can be described by continuum hydrodynamic equations for low-Reynolds-number steady-state flow of a dilute solution of an incompressible Newtonian fluid, the governing equations are 29 where u, j, p, and c are the solution velocity, solute current density, pressure, and solute concentration, respectively, η is the solution shear viscosity, and U is the solute-membrane interaction potential.D and µ are the solute diffusivity and mobility respectively, which we assume are related by the Einstein relation, µ = D k B T , 30 where k B is the Boltzmann constant and T is the temperature.U is the interaction potential per solute molecule and so −c∇U is a body force per unit volume acting on the fluid due to the solute-membrane interactions.
U is assumed to depend on the position in the fluid relative to the membrane surface.For a neutral solute, U typically depends on the distance from the surface 28,29 and has a range on the order of the solute molecular diameter.Further assuming that advection of the solute is negligible compared with diffusion (i.e.low Péclet number flow), Eq. ( 2) for the solute flux simplifies to The solution velocity and solute flux are assumed to satisfy the usual no-slip and zero flux boundary conditions at the membrane surface, i.e. u = 0 and n • j = 0, where n is the unit vector normal to the membrane surface.
We note that a similar approach based on the widely used Poisson-Nernst-Planck-Stokes equations for electrolytes, 30 in which the electric potential energy plays an analogous role for the electrolyte that the interaction potential U does for a neutral solute, could be used to extend this study to concentration-gradient-driven electrolyte transport.However, such an extension is non-trivial, as the electric potential must be determined by solving an additional coupled differential equation (the Poisson equation) that depends on the solute (electrolyte) concentration, rather than being specified.Thus, we leave this extension to electrolytes to future work.
Consider the fluid flow induced by a concentration difference, ∆c = c H − c L , between the two sides of the membrane, with the pressure far from the membrane the same on both sides of the membrane, i.e. p H = p L = p ∞ , as shown in Fig. 1.
Our derivation uses a combination of cylindrical (r, z, φ) and oblate spheroidal (ζ, ν, φ) coordinates, where z = aνζ, r = a (1 + ν 2 )(1 − ζ 2 ), and θ is the angle about the z axis (0 ≤ ν < ∞, −1 ≤ ζ ≤ 1, and 0 ≤ φ < 2π). 24,31Writing the solute concentration in the presence of a concentration gradient as where c s (ζ, ν) is to be determined, and inserting this expression into Eq.( 4) and the boundary condition for the solute current density j gives with the boundary condition n • ∇c s = 0 at the membrane surface.If the solute-membrane interaction potential is small relative to the thermal energy (U k B T ), Eq. ( 6) reduces to Solving this equation, subject to the boundary conditions on c s at the membrane surface and far from the membrane (c s → c H = c ∞ + ∆c 2 and c L = c ∞ − ∆c 2 in the upper and lower half-planes, respectively, where U → 0), gives 31 We have verified using finite-element numerical simulations (see Fig. S3 of the supplementary material) that Eq. ( 8) with Eq. ( 5) accurately describe the solute concentration even when U is several times k B T .A possible reason why Eq. ( 8) appears to be accurate outside the regime for which it was derived is that the second term in Eq. ( 6) that was neglected to arrive at Eq. ( 7) and ( 8) can be small even if the magnitude of the potential U is large: for example, for c s given by Eq. ( 8) this term is zero for a potential that is a function only of ζ or in the pore mouth (at ν = z = 0) for a potential that is a function only of distance from the membrane, due to the orthogonality of ∇U and ∇c s in these cases.
The fluid flow induced by the concentration gradient can be obtained from the reciprocal theorem for steady incompressible creeping flow, 32 which allows the fluid flow due to a body force F acting on the fluid to be related to the pressure-driven flow in the same pore geometry, 24 for which an analytical solution exists for the fluid velocity through a circular aperture. 22,32As shown by Mao et al. 24 for the related problem of electroosmosis through a circular aperture, where ū is the fluid velocity induced by a pressure difference ∆p = p H − p L for ∆c = 0 in the system geometry in Fig. 1 and the integral is over the volume V occupied by the fluid.For concentration-gradient-driven flow described by Eq. ( 1), F = −c∇U .The pressure-driven flow velocity can be obtained from the stream function ψ = − a 3 6πη (1 − ζ 3 ) ∆p for the flow 32 using ū = 1 r φ × ∇ψ, 32 where φ is the unit vector in the φ direction, as Inserting this expression for ū and F = −c∇U into Eq.( 9) (with c given by Eqs. ( 5) and ( 8)) and making use of with Equation ( 12) is the main result of this work.
Furthermore, the solute flux density can be obtained, using Eqs.( 4) and ( 5), as The solute flux across the membrane is where the unit vector normal to the pore mouth is n = ν = ẑ and the surface integral is over the pore aperture.Evaluating the solute flux at the pore mouth (at ν = z = 0) using Eqs ( 8) and ( 13) yields = −2D∆c

B. Limiting cases and scaling behavior
The diffusioosmotic mobility predicted in Eq. ( 12) depends crucially on the range of the interaction potential U , which we will call λ.However, the term exp(−U/k B T ) − 1 in Eq. ( 12) is averaged spatially with a geometry-dependent weight, with a complicated dependence on the (oblate-spheroidal) coordinates ζ and ν.As a consequence, the mobility, and its scaling with the pore radius a and interaction range λ, may depend on specific details of the geometry dependence of the interaction.Therefore, we consider various limiting cases for the geometry of the potential and the consequences for the dependence of the scaling with the pore radius and interaction range.
1. Case of potential that depends only on the distance to the membrane surface and/or pore edge In most circumstances, the potential U is expected to be a function of the distance d from the membrane surface.However, the integral in Eq. ( 12) for the mobility cannot be simplified due to the geometrical interplay between the various variables.
Let us assume for simplicity that the solute excess/depletion at the membrane surface can be represented by a step function as a function of the distance d from the surface, i.e.
where α characterizes the solute excess close to the membrane surface (α > 1 for solute adsorption and α < 1 for solute depletion).While the mobility still cannot in general be evaluated analytically in this case, analytical solutions exists in certain limits.In particular, for a solute-wall interaction range much larger than the pore radius (λ a), e −U (d)/k B T can be approximated as a constant independent of the coordinates, and we find from Eq. ( 12) in this limit that On the other hand, if λ a, the integral in Eq. ( 12) can be approximated as Details of the derivation of this equation can be found in the supplementary material.
A similar calculation can be performed for the case in which the interaction originates only from the pore edge, and U depends on the distance to the edge.As detailed in the supplementary material, the mobility in this case is Solute flux.Similarly, using Eq. ( 17) in Eq. ( 16) and noting that d = a − r in the pore mouth at z = 0, the solute flux across the membrane for the step-function potential can be shown to be for λ ≤ a and for λ a.

Case of potential that depends only on ζ
In the case in which the potential U is a function of the ζ coordinate only (see Fig. 1), the diffusioosmotic mobility in Eq. ( 12) simplifies to Assuming that U only depends on ζ is a stringent condition in terms of symmetry, but interestingly such a result is expected for a circular aperture in a planar membrane at a fixed electrostatic potential (however, in the absence of screening).
To determine how Q scales with λ, particularly in the limit λ a, we can use the In addition, U is only non-zero when Eq. ( 23) is non-zero.Thus, in this limit, U (ζ) ∼ U [aζ 2 /(2λ)] and the integral in Eq. ( 23) becomes where x ≡ [a/(2λ)] 1/2 ζ.So the mobility in this case scales as κ DO ∼ a 3 (λ/a) 3/2 ∼ (aλ) 3/2 .Equation ( 23) can also be rewritten in radial coordinates, by focusing on the pore mouth (ν = z = 0), as where r is the distance to the center of the pore in the membrane plane (z = 0).As an alternative approach to predicting the dependence of the mobility on λ, we again restrict ourselves for simplicity to a step-function interaction versus distance d = a − r from the pore edge, i.e.
Inserting Eq. ( 26) into Eq.( 25) gives the diffusioosmotic mobility for the step-function potential as For a small interaction range λ compared with the pore radius, the predicted scaling of the mobility is which is identical to the scaling predicted directly from Eq. (23).

Case of potential that depends only on ν
In the case in which the potential U is a function of the ν coordinate only (see Fig. 1), the diffusioosmotic mobility in Eq. ( 12) simplifies to Following similar reasoning to the previous section, in terms of the distances, d 1 and d 2 , from a point with a given ν value to the foci of the ellipses or hyperbolae in Fig. 1, . So a potential U (ν) that depends only on ν depends only on the average distance, (d 1 + d 2 )/2.Therefore, assuming that the potential has a distance range λ entails Thus, in this limit, U (ν) ∼ U [aν 2 /(2λ)] and the integral in Eq. ( 29) becomes where x ≡ [a/(2λ)] 1/2 ν.So the mobility in this case scales as κ DO ∼ a 3 (λ/a) 1/2 ∼ a 5/2 λ 1/2 .

III. NUMERICAL RESULTS
To validate the theory, finite-element method (FEM) simulations of concentrationgradient-driven flow were carried out using Comsol Multiphysics (version 4.3a) 33 through pores with various radii and solute-membrane interactions.Here we consider that the solute interacts with the membrane via a potential that depends on the (shortest) distance d to the membrane surface.Accordingly, the solute-membrane interaction potential was modelled using a hyperbolic tangent function, defined by parameters and λ describing the strength and range of the potential.In all simulations, the average solute concentration c ∞ was 10 −3 σ −3 and the solute diffusivity D was σ 2 /τ , and unless otherwise stated the aperture radius a was 10 σ and λ was σ, where σ is the unit of length (σ can be regarded as the diameter of a fluid molecule) and τ = ησ 3 /(k B T ) is the unit of time.Details of the FEM simulations, which all correspond to low-Péclet number flow, are given in the supplementary material.
We have quantified the concentration-gradient-driven solution and solute fluxes measured in the numerical simulations and predicted by our theory in terms of the diffusioosmotic mobility κ DO defined by and the solute permeance P s defined by The equations that we have derived for the solution and solute fluxes (Eqs.( 12) and ( 15)) predict that the fluxes are linearly related to the concentration difference ∆c and thus that the conductances and resistances defined in Eqs. ( 33) and ( 34) are independent of ∆c.We have verified that this is indeed the case for the range of concentration differences studied in the numerical simulations (∆c = 10 −6 to 3 × 10 −4 σ −3 = 10 −3 to 0.3 c ∞ ), as shown in Fig. S5 of the supplementary material.
Figure 2 shows the diffusioosmotic mobility κ DO from the simulations and theory of a circular aperture as a function of aperture radius a and solute-membrane interaction range λ for two different values of the solute-membrane interaction strength , k B T /10 and k B T , with all other parameters kept constant.The theory curves were calculated by numerically evaluating the integral in Eq. ( 12) with the solute-membrane potential U in Eq. ( 32).The sign of κ DO has been defined so that a positive and negative values correspond to fluid flow in opposite and same direction, respectively, to the applied concentration gradient.
Hence, for solute depletion at the membrane surface ( > 0) the flow is towards higher solute concentration (κ DO < 0), whereas for solute adsorption the flow is towards lower concentration (κ DO > 0). 34gure 2 shows good quantitative agreement between the theory and simulations for the variation of κ DO with all relevant parameters for k B T .Since the theory assumes a weak potential in deriving Eq. ( 8) for the solute concentration, we indeed find deviations between the prediction and the simulations as the magnitude of the solute-membrane potential increases.Nevertheless, the agreement is reasonable well beyond the regime of validity of this approximation.Note that for the values of in Fig. 2 and so from Eq. ( 18) or (19) κ DO is approximately proportional to , but this scaling is not expected in general and already starts to break down for = k B T .
Figure 2 also compares the simulation results with the approximate scaling with the pore radius a and solute-membrane interaction range λ predicted by the theory.For a λ, Eq. ( 18) predicts that κ DO is proportional to a 3 and independent of λ, which is evident in the scaling for small a in Fig. 2(a) and in the saturation at large λ in Fig. 2(b).On the other hand, for a λ, Eq. ( 19) predicts scaling of κ DO with aλ 2 , which is seen to hold in the large-a regime in Fig. 2(a) and in the small-λ regime in Fig. 2(b).
Figure 3 shows the analogous comparison between the FEM simulations and theory for the solute permeance P s .The theory curves were calculated by numerically evaluating the integral in Eq. ( 16) with the solute-membrane potential U in Eq. (32).As for the diffusioosmotic conductance, the theory accurately captures the simulated solute permeance for k B T , with deviations between the theory and simulations becoming evident for magnitudes of the solute-membrane potential greater than k B T .For the parameters used in Fig. 3(a), the first term in Eq. ( 21) dominates and so the permeance P s shows the expected linear scaling with the pore radius a.For λ a, Eq. ( 21) predicts that P s varies from its value at λ = 0 with a scaling as √ λ, which is evident in Fig. 3(b).

IV. DISCUSSION
An interesting outcome of the previous results is that the diffusioosmotic mobility κ DO = −Q/∆c across a circular aperture in ultrathin membrane is strongly dependent on the details of the interaction of the solute with the membrane.We have shown in particular that the mobility scales with the pore radius a and interaction range λ as a 3−γ λ γ , with an exponent γ that depends on the underlying symmetries of the potential: when λ a, for a potential that depends only on the ζ coordinate (in the oblate-spheroidal system, see Fig. 1) an exponent γ = 3/2 is predicted, for a potential that depends only on the ν coordinate an exponent γ = 1/2 is predicted, whereas for a potential depending on the distance to the membrane a value γ = 2 is found.
It is also interesting to compare the results for the circular aperture with those obtained in long cylindrical pores (e.g. as a model for nanotubes).As derived in detail in the supplementary material, the diffusioosmotic mobility of a long cylindrical pore of length L is proportional to (aλ) 2 /L for λ a.When compared to the case leading to an exponent γ = 2, the scaling of κ DO with pore size and interaction range λ is therefore recovered by replacing the length of the nanopore L with the pore size a, which is indeed expected for entrance effects.However, as shown with the case leading to an exponent γ = 3/2 or 1/2 for the diffusioosmotic mobility, this situation is not universal.The different scaling relationships derived for κ DO for a circular aperture and a long cylindrical pore are summarized in Table I.
These result are relevant for transport through finite-length pores, for which the total resistance to flow can often be accurately given by the sum of the resistance due to the pore interior and that due to the pore ends, which can be approximated by that of a circular aperture. 13,23,26e predicted scaling behavior of the diffusioosmotic mobility and solute permeance of a circular aperture for concentration-gradient-driven flow differ markedly from the scaling behavior derived previously for other types of flows in the same system geometry.For example, the hydraulic conductance (solution flux per unit pressure difference) in pressure-driven fluid flow through a circular aperture has been shown 22,23,32 to be proportional to a 3 in contrast to the proportionalty with a 3−γ where γ > 0 for the diffusioosmotic mobility in concentration-gradient-driven flow in the limit λ a.
On the other hand, for λ a, the diffusioosmotic mobility shows the same a 3 scaling as the hydraulic conductance.The equivalent scaling of the hydraulic conductance and diffusioosmotic mobility in this limit is because concentration-driven-flow in which the solute-membrane interaction range is larger than the aperture radius is equivalent to osmotic transport through a semipermeable membrane, with the osmotic pressure gradient due to the concentration gradient playing an equivalent role to the pressure gradient in pressure-driven flow. 35Likewise, the scaling behavior predicted here for the diffusioosmotic mobility differs from the electroosmotic conductance for electric-field-driven fluid flow of an electrolyte, which has been shown to be proportional to a 3 /λ D for a λ D , 24 where λ D is the Debye length characterizing the electric double layer width, equivalent to λ here; in the same limit, scaling of the diffusioosmotic mobility with a 3 is predicted here.
The electrical conductance across a circular aperture in electric-field-driven transport of an electrolyte has been shown to be proportional to a for an uncharged membrane, 25 with the addition of surface charge to the membrane only changing the length scale in the scaling relationship to an effective radius a eff that is the sum of the aperture radius and the Dukhin length characterizing the ratio of the surface electrical conductivity to the bulk electrical conductivity, 26 without changing the scaling exponent.This scaling differs from that derived for the solute permeance for λ ≤ a, for which the equivalent effective radius is , in which the second term in the sum depends on both the solute-membrane interaction range λ and pore radius a.
We can consider the implications of our theory for realistic systems, in particular for diffusioosmotic transport of electrolyte solutions.Extending the theory to electrolytes is desirable for applications, such as salinity-gradient-driven energy conversion 4,20,27 and desalination, 3,14,15,17 but it is technically difficult, so we leave this derivation for future work.
However, as a rule of thumb, one may expect that this would amount to replacing the solutemembrane interaction range λ by the salt-concentration-dependent Debye length λ D ∝ c −1/2 salt .A counterintuitve outcome of the non-universal dependence of the diffusioosmotic mobility as a function of the interaction range λ is a possible impact on the dependence of the mobility on the salt concentration.In a long pore under a salinity difference ∆c salt , the solvent flux is predicted to behave as Q = K DO ∆ log c salt , with the mobility K DO ∝ c 0 salt πa 2 /L for a long pore of length L, 27 i.e. no dependence on salt concentration.This is due to the dependence of the diffusioosmotic mobility κ DO ∼ λ 2 D ∼ c −1 salt , so that Q ∼ ∆ log c salt .Now, for a circular aperture, we have shown that the dependence of the diffusioosmotic mobility κ DO on the interaction range scales as κ DO ∼ λ γ with a non-universal exponent γ.Thus, the case of γ = 2 (occurring for potentials that depend on the distance to the membrane) will lead to K DO ∝ c 0 salt as for the long-pore case; but cases with γ = 2 (as highlighted above in two cases), will lead to K DO ∝ c 1−γ/2 salt , exhibiting therefore a curious dependence on c salt .

V. CONCLUSIONS
In summary, we have derived general equations and scaling relationships as a function of pore radius and solute-membrane interaction strength and range for the solution and solute fluxes induced by a solute concentration gradient through a circular aperture in an ultrathin planar membrane.We have shown, by comparing with finite-element numerical simulations, that the equations accurately quantify the fluid fluxes when the solute-membrane interaction strength is small compared with the thermal energy k B T .In the limit of a solute-membrane interaction range much smaller than the pore radius, the theory predicts a non-universal dependence of the fluid fluxes on the pore radius and interaction range.These results have significant implications for applications involving concencentration-gradient-driven flow in membranes in which the thickness is on the order of the pore size, such as those made from 2D nanomaterials, notably in the context of blue energy harvesting.

SUPPLEMENTARY MATERIAL
The supplementary material contains derivations of scaling laws for diffusioosmosis through a circular aperture in an ultrathin planar membrane, the theory of concentrationgradient-driven flow in a long cylindrical pore, and further details and supplementary results of finite-element numerical simulations.

B. Solute interaction with pore mouth only
Assuming that the solute interacts only with the part of the membrane surface at the pore mouth, the solute-membrane potential U (d) will only depend on the distance from the surface at radial coordinate r = a.In this case, the distance is d = (a − r) 2 + z 2 , ∀r, z.Thus, following analogous steps to those in the previous section, if d ≤ λ a, d a 2 (ν 2 + ζ 2 ) and Inserting Eq. (S15) into Eq.(S1) yields Thus, the scaling of Q with the pore radius a and solute-membrane interaction range λ and strength (measured by α) is the same as in the case of solute interactions with the entire membrane.

S2. THEORY OF CONCENTRATION-GRADIENT-DRIVEN FLOW IN A LONG CYLINDRICAL PORE
Here we derive equations for concentration-gradient-driven flow of a dilute solution in a long cylindrical pore (Fig. S1) assuming the same governing equations used in the main text for flow through a circular aperture in an ultrathin planar membrane.We assume that the pore length L is much larger than the pore radius (L a), so that entrance effects can be ignored.The solute-wall interaction potential U (r) in this case is only a function of the radial coordinate r.
Since L a, the radial solute flux will be small relative to the axial solute flux (0 ≈ j r j z ), so from Eq. ( 4) of the main paper, where σ was the unit of length.We verified that the measured solution and solute fluxes did not change significantly with a finer mesh.Table S1 lists the boundary conditions that were used to solve the equations and Table S2 lists the parameters used in the simulations.
Several supplementary results from the FEM simulations are given in the figures below.
Figure S3 shows the simulated solute concentration c and the ratio of the simulated to S7

FIG. 1 .
FIG. 1. Schematic of flow of a solution through a circular aperture of radius a in an infinitesimally thin planar wall.The origin in cylindrical (r, z) coordinates is at the center of the aperture, as indicated, and the axis of symmetry is the depicted z axis.The solute concentration c and solution pressure p far from the membrane are c H and p H , respectively, in the upper half-plane and c L and p L , respectively, in the lower half-plane.Q and J are the solution and solute fluxes, respectively.Contours of constant ζ and ν in oblate-spheroidal (ζ, ν) coordinates are also shown as dashed lines, with unit vectors shown at one point in space.
relationship between ζ and the distances, d 1 and d 2 , from a point with this ζ value to the two foci of the hyperbolae or ellipses of constant ζ or ν that are shown in Fig. 1.For an infinitesimally thin membrane, these foci are located at the pore edge (at r = a), and thus 1 − ζ 2 = [(d 1 − d 2 )/(2a)] 2 .So a potential U (ζ) that depends only on ζ depends only on the relative distance, d 1 − d 2 .Furthermore, assuming that the potential has a distance range λ implies that U

FIG. 2 .
FIG. 2. Diffusioosmotic mobility κ DO versus (a) pore radius a (with λ = σ) and (b) solutemembrane interaction range λ (with a = 10σ) for solute-membrane interaction strength = k B T /10 or k B T from FEM simulations (points) and theory (solid lines).The dashed lines show scaling with various powers of a and λ.
FIG. S2: Schematic (not to scale) of the two-dimensional axisymmetric computational domain used in the FEM calculations.Solid lines denote solid-liquid boundaries and dashed lines liquid boundaries.The geometry has rotational symmetry about the boundary AH.

∆c σ −3 10 − 6 - 3 × 10 − 4 reservoir radius/length w σ 1000
= n • u = n • ∇u = 0 AB c = c H , p = p ∞ GH c = c L , p = p ∞ BC and FG n • j = n • u = n • ∇u = 0 CD,DE, and EF n • j = u = 0 TABLE S2: Parameters used for FEM calculations.The units of length, energy, and time are σ, k B T , and τ = ησ 3 /(k B T ), respectively, where T is the temperature and η is the solution shear viscosity (σ can be considered the diameter of a fluid molecule).Where a range of values is given, the parameter was fixed at the value in parentheses unless otherwise indicated.theoretical solute concentration, c/c theory , where c theory is the solute concentration derived in the theory (given by Eqs.(5) and (8) of the main paper), near the pore entrance for an intermediate solute concentration difference between the reservoirs (similiar results were obtained for the range of parameters simulated, with the maximum deviation between c and c theory of around ±4% observed for = −3k B T and ∆c = 3 × 10 −4 /σ 3 ).These results indicate that the theory accurately predicts the nonequilibrium solute concentration in this FIG. S3: 2D plots of the ((a), (b)) solute concentration c (in units of σ −3 ) from FEM simulations and ((c), (d)) c/c theory , where c theory is the solute concentration derived in the theory (given by Eqs.(5) and (8) of the main paper), as a function of the r and z coordinates near the pore entrance for a solute-membrane interaction strength of ((a),(c)) k B T /10 and ((b),(d)) 3k B T (for a = 10σ and λ = σ).The solute concentration difference ∆c between the reservoirs was 10 −5 σ −3 in both cases.
FIG. S4: (a) Axial velocity u z and (b) axial solute current density j z from FEM simulations versus the radial coordinate r at the pore mouth at z = 0 for varying solute-membrane interaction strength (for a = 10σ and λ = σ).The solute concentration difference ∆c between the reservoirs was 10 −5 σ −3 in all cases.(The arrows indicate the direction of increasing parameter values.)

TABLE I .
Scaling of the diffusioomotic mobility κ DO with pore radius a and solute-membrane interaction range λ for a circular aperture in a 2D membrane, as well as with pore length L for a long cylindrical pore, for different functional forms of the solute-membrane potential U in the

TABLE S1 :
Boundary conditions used to solve the continuum equations in the FEM simulations.The vector n is the surface normal.