留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

Anomalous diffusion in branched elliptical structure

SuleimanKheder ZhangXuelan WangErhui LiuShengna ZhengLiancun

SuleimanKheder, ZhangXuelan, WangErhui, LiuShengna, ZhengLiancun. Anomalous diffusion in branched elliptical structure[J]. 机械工程学报, 2023, 32(1): 010202. doi: 10.1088/1674-1056/ac5c39
引用本文: SuleimanKheder, ZhangXuelan, WangErhui, LiuShengna, ZhengLiancun. Anomalous diffusion in branched elliptical structure[J]. 机械工程学报, 2023, 32(1): 010202. doi: 10.1088/1674-1056/ac5c39

Anomalous diffusion in branched elliptical structure

doi: 10.1088/1674-1056/ac5c39
  • Diffusion is an essential means of transportation in biological, physical, and chemical systems that is used by individuals to achieve their functions, for example, but not limited to, cellular processes, metabolism, and conduction. [ 15] Realizing the features of diffusion, which depend on its microscopic dynamics and environmental space, will provide an important tool to control these processes, besides its possible applications in nanotechnology and nanomedicine. According to the theory of Brownian motion, the probability distribution is Gaussian and mean square displacement (MSD) of the particle position grows linearly with time 〈 r 2( t)〉 = 2 Dt where D is the diffusion coefficient. However, deviations from these usual behaviors have been observed theoretically and experimentally, anomalous (non-Brownian) diffusion, [ 610] Brownian yet non-Gaussian diffusion, [ 1113] and anomalous yet non-Gaussian diffusion. [ 14, 15] The nonlinear dependence of time characterizes the MSD in the anomalous dynamics, namely, 〈 r 2( t)〉] ∼ t α with α ≠ 1. The diffusion dynamics with transport exponent 0 < α < 1 correspond to sub-diffusion dynamics which are observed in artificially crowded liquids, [ 16] lipid bilayer membranes, [ 15, 17] the cytoplasm of biological cells, [ 18] extracellular space, [ 19] and in hydrology, [ 20] to mention a few. Whereas, it defines super-diffusion when α > 1, which is encountered in active systems such as molecular motor transport in cells. [ 21, 22] Different hypotheses have been suggested to understand the mechanism underlying sub-diffusion dynamics, which depends on the situation at hand, among them, trapping models of energetic (binding sites) or geometric nature, labyrinthine environment (geometrical disorder), temporal correlation due to slow mode (viscoelastic environments), space-dependent diffusivity due to the local porosity of the medium or temperature field in the sense of weak gradients, which can be modelled by continuous time random walk (CTRW), random walk on fractals, fractional Brownian motion (FBM) or fractional Langevin equation (FLE), and heterogeneous diffusion process (HDP), respectively. [*,23-27]

    Disordered materials in nature are abundant, including to mention a few the heterogeneous porous media, extracellular space in brain tissue, intra-tissue in muscles, cell cytoplasm, and glasses. [ 2835] Diffusion in disordered media is quite affected by its geometrical structure, for instance, dead ends, bottlenecks, and backends. Since deterministic fractals and percolation clusters are elegant models for the geometrical structure of disorder systems, Pierre-Gilles de Gennes has suggested modelling diffusion in these systems through performing random walks on percolation clusters and fractals, and coined the term “ant in the labyrinth”. [ 26] Random walks in these structures lead to anomalous diffusion. At the percolation threshold, a percolation structure can be idealized as a single infinite cluster, consisting of a conducting path, which corresponds to a backbone, and side branches, or fingers with dangling bonds. At some level of idealization, this structure corresponds to a comb. [ 36] Later, the comb-like structure consisting of the one-dimensional axis (backbone) with perpendicular sides branches (fingers) has been used to understand the underlying mechanism of anomalous diffusion in percolation clusters: while the particle moving randomly along the backbone may trap inside the dead-ends (fingers), performing a Brownian motion, until it returns by chance to the backbone where it can escape from the finger. Along the y-direction (the so-called fingers), the particles perform normal diffusion, while the mean square displacement (MSD) along the backbone grows non-linearly with a time x 2 ( t ) t , and sub-diffusion behavior dominates the process. [ 26, 37, 38]

    To study the transport properties of diffusion process on comb-like structure, different approaches have been applied, e.g., a mesoscopic approach used by Weiss and Havlin, [ 37] a macroscopic approach used by Arkhincheev and Baskin, [ 39] and a microscopic approach used by Méndez et al. [ 40] Arkhincheev and Baskin have expressed the diffusion process in the comb by using the following Fokker–Planck equation: [ 39] P ( x , y , t ) t = δ ( y ) D 2 P ( x , y , t ) x 2 + D 2 P ( x , y , t ) y 2 . Here, D x = δ ( y ) D and D y = D are the diffusion coefficients in the x and y directions, respectively, and the δ-function in the Fokker–Planck operator L FP = δ ( y ) D 2 x 2 + D 2 y 2 means that the diffusion along the x direction (the so called backbone) is allowed only at y = 0.

    Recently, comb-like model has been used to describe anomalous diffusion in spiny dendrites, [ 41] the mechanism of super-diffusion of ultracold atoms in a one-dimensional as a phenomenology of experimental study, [ 42] as well as anomalous diffusion in porous materials. [ 43] Moreover, several complex extensions of comb-like structure, through introducing a modification on the geometry of the backbone or finger shape of the structure, is introduced to describe geometrically induced complex diffusion in nature, for example, comb with a finite finger length, [ 44, 45, 40] cylindrical comb, [ 46, 47] random comb models, [ 48, 49] comb with ramified teeth, [ 40, 41] fractal mesh and grid structures, [ 5052] and more complex structure. [ 5357] Diffusion along the backbone in both two and three-dimensional comb-like structures with a finite finger is a transient sub-diffusion followed by normal diffusion at long times, and the same result is obtained in using a trap model of energetic nature. [ 58] Whereas, diffusion along the backbone in a three-dimensional comb-like structure with an infinite two-dimensional branch is ultra-slow diffusion/enhanced sub-diffusion if normal/sub-diffusion occurs inside branches. [ 46, 47] In addition, recent theoretical studies have been demonstrating that performing different microscopic dynamics in a comb-like model leads to different types of anomalous diffusion. [ 5963] In most previous studies, the geometric properties of the backbone have not been considered. However, the interior of disordered systems are narrow, tortuous, curved channels. [ 26, 64, 65]

    In the last decade, the problem of diffusion in curved confined tubes has attracted remarkable attention of some researchers, as it was found that the behavior of the diffusion process is affected by its geometrical properties, like curvature and torsion, of these tubes. [ 6669] This problem is related to diverse kinds of phenomena, for example not limited to, diffusion of biomolecules in biological cell, [ 70] transport of solute in aquifers and porous medium, and fluxon dynamics in Josephson junction. [ 71] The aim of this paper is to study the diffusive dynamics in branched channels with gradient curvature. Branched elliptical-like channels are common in biology, including to mention a few, the extracellular spaces in the biological tissues (brain, tumors, muscles), intracellular environments, and other porous media, [ 19, 31, 7278] in which diffusion process in these structures are important for other biological processes like extracellular ionic buffering, and delivery of drugs and metabolites, etc. [ 79, 80] Moreover, they are simple models of curved confined sub-diffusive systems with gradient curvature. Motivated by the foregoing, we analyze the dynamics of the diffusion process in an elliptical comb-like structure (see Fig. 1) with reflected boundary conditions. The elliptical motion takes place only for a fixed radius r = R and is interspersed with a radial motion inward and outward of the ellipse. The Brownian radial motion corresponds to diffusion of particles in dead-end spaces, while elliptical backbone describes the dynamics along elliptical channels. We investigate by numerical simulation how probability distribution function (PDF) and mean square displacement (MSD) are affected by geometrical constraints in the elliptical comb-like structure.

    1  Schematic picture of an elliptical comb-like structure, where the radii are continuously distributed over the ellipse of radius R = ρ( θ) with θ ∈ [0,2 π].

    The rest of the paper is organized as follows. We firstly formulate the mathematical description of the problem in Section 2. The numerical method and example on the verification of the numerical algorithm are collected in Section 3. Then results and discussion are illustrated in Section 4. The conclusions are given in Section 5.

    The structure we study here consists of elliptical channel with radii dead ends distributed along it and directed in/outside (see Fig. 1), in which an unbiased diffusion process occurs inside it. Mathematically, suppose an ellipse whose boundary has the polar equation ρ ( θ ) = b / 1 e 2 cos 2 θ , where 0 < e < 1 is the eccentricity of ellipse, b is the shortest diameter with constant value, and θ ∈ [0,2 π]. Then, its domain can be obtained through stretched radius r such that ϒ = ( θ) with θ ∈ [0,2 π], and r ∈ [0,1], where we can move from Cartesian coordinates to stretched polar coordinates through the transformation [ 81] x = r ρ ( θ ) cos θ , y = r ρ ( θ ) sin θ . The distinctive feature of diffusion process in such structure is that the particle motion along θ direction is possible only when r = R and R is a constant. We mean that the tangent diffusion coefficient D θ is different from zero only at r = R, i.e., D θ = D δ ( r R ) with D is a constant. Whereas, the diffusion is free in the radial direction situation and its diffusion coefficient can be written as D r = D. By adjusting the parameter e the elliptical channel (backbone) has different shapes, in which the curvature of the ellipse, κ = a b / ( a 2 sin 2 θ + b 2 cos 2 θ ) 3 with a 2 = b 2/(1 – e 2), decreases as the value of e increases (see Fig. 1). In order to assess how particles diffuse in the structure probability distribution function (PDF) P( r, θ, t), marginal PDF P( r, t), and mean square displacement (MSD) along θ are a useful quantities. In the classical comb model, PDF P( r, θ, t) can be described by Fokker–Planck equation [ 39] t P = δ ( y ) D x 2 P + D y 2 P . The δ function in the Fokker–Planck operator L FP = δ ( y ) D x 2 + D y 2 means that the diffusion along the x direction is allowed only at y = 0. In what follows, we derive the diffusion equation in a stretched polar coordinate system, then find the differential equation that expresses the diffusion process in our structure.

    The gradient operator can be written as follows: = u ^ r ( 1 ρ ( θ ) r ) + u ^ θ ( ρ ( θ ) ρ 2 ( θ ) r + 1 r ρ ( θ ) θ ) . Then, the balance equation and the first Fick’s law can be written respectively as follows: P t = J = [ 1 ρ ( θ ) J r r + J r r ρ ( θ ) ρ ( θ ) ρ 2 ( θ ) J θ r + 1 r ρ ( θ ) J θ θ ] , J = D P = [ u ^ r ( D r ρ ( θ ) P r ) + u ^ θ ( D θ ρ ( θ ) ρ 2 ( θ ) P r + D θ r ρ ( θ ) P θ ) ] . Then the second Fick’s equation reads P t = 1 ρ 2 ( θ ) ( D r + D θ ρ 2 ( θ ) ρ 2 ( θ ) ) 2 P r 2 + 1 r ρ 2 ( θ ) ( D r D θ ρ ( θ ) ρ ( θ ) + 2 D θ ρ 2 ( θ ) ρ 2 ( θ ) ) P r 2 D θ ρ ( θ ) r ρ 3 ( θ ) 2 P r θ + D θ r 2 ρ 2 ( θ ) 2 P θ 2 .

    The dynamic in the structure can be described by the following diffusion equation (substitute ρ ( θ ) = b 1 e 2 cos 2 θ , D r = D, and D θ = D δ ( r R ) into Eq. ( 5)): P t = D [ ( 1 e 2 cos 2 θ b 2 ) 2 P r 2 + ( 1 e 2 cos 2 θ b 2 r ) P r ] + D [ ( 4 e 2 cos 2 θ ( 1 e 2 cos 2 θ ) e 4 sin 2 2 θ 4 b 2 r ( 1 e 2 cos 2 θ ) ) P r + ( e 4 sin 2 2 θ 4 b 2 ( 1 e 2 cos 2 θ ) ) 2 P r 2 + ( e 2 sin 2 θ b 2 r ) 2 P r θ + ( 1 e 2 cos 2 θ b 2 r 2 ) 2 P θ 2 ] × δ ( r R ) , with initial condition P ( r , θ , 0 ) = δ ( r R ) δ ( θ ) , and boundary conditions θ P ( R , + π , t ) = θ P ( R , π , t ) = 0 , r P ( L , θ , t ) = r P ( R , θ , t ) = 0.

    The initial–boundary value problem ( 6)–( 8) can be solved numerically by using the following procedure: Firstly, in order to divide the spatial and temporal domain into grid, we define r i = ih r , θ j = jh θ , and t l = where i = 1 , N r + 1 ¯ , j = 1 , N θ + 1 ¯ , and l = 1 , N T + 1 ¯ , besides h r = R′ / N r , and h θ = π/ N θ are space steps, τ = T/ N T is time step. Let P i , j l be the numerical solution of Eq.( 6) at point ( r i , θ j , t l ), then the first and second order derivatives in the governing equation can be discretized as follows: t P ( r i , θ j , t l ) = ( P i , j l + 1 P i , j l ) / τ , r P ( r i , θ j , t l ) = ( P i + 1 , j l P i , j l ) / h r , θ P ( r i , θ j , t l ) = ( P i , j + 1 l P i , j l ) / h θ , θ 2 P ( r i , θ j , t l ) = ( P i , j + 1 l 2 P i , j l + P i , j 1 l ) / h θ 2 , r 2 P ( r i , θ j , t l ) = ( P i + 1 , j l 2 P i , j l + P i 1 , j l ) / h r 2 , r θ 2 P ( r i , θ j , t l ) = ( P i + 1 , j + 1 l P i + 1 , j l P i , j + 1 l + P i , j l ) / h r h θ . Now, applying the previous approximate scheme to Eq. ( 6) yields the following iterative equations: P i , j l + 1 = [ λ j 1 + λ i , j 3 + ( λ j 2 + λ i , j 4 λ i , j 5 ) δ ( r i R ) ] × P i + 1 , j l + [ λ j 1 + λ j 2 δ ( r i R ) P i 1 , j l + [ 1 2 λ j 1 λ i , j 3 ( 2 λ j 2 + λ i , j 4 λ i , j 5 + 2 λ i , j 6 ) × δ ( r i R ) ] P i , j l + [ ( λ i , j 6 λ i , j 5 ) δ ( r i R ) ] P i , j + 1 l + [ λ i , j 6 δ ( r i R ) ] P i , j 1 l + [ λ i , j 5 δ ( r i R ) ] P i + 1 , j + 1 l , where λ j 1 = D ( 1 e 2 cos 2 θ j ) τ b 2 h r 2 , λ j 2 = D e 4 sin 2 ( 2 θ j ) τ 4 b 2 h r 2 ( 1 e 2 cos 2 θ j ) , λ i , j 3 = D ( 1 e 2 cos 2 θ j ) τ b 2 r i h r , λ i , j 4 = 4 D e 2 cos 2 θ j ( 1 e 2 cos 2 θ j ) τ D e 4 sin 2 ( 2 θ j ) τ 4 b 2 r i h r ( 1 e 2 cos 2 θ j ) , λ i , j 5 = D e 2 sin 2 ( 2 θ j ) τ r i b 2 h r h θ , λ i , j 6 = D ( 1 e 2 cos 2 θ j ) τ b 2 r i 2 h θ 2 .

    In order to test the correctness of numerical solution, we introduce a source function f( r, θ, t) in Eq. ( 6), and this yields the following equation: P t = D [ ( 1 e 2 cos 2 θ b 2 ) 2 P r 2 + ( 1 e 2 cos 2 θ b 2 r ) P r ] + D [ ( 4 e 2 cos 2 θ ( 1 e 2 cos 2 θ ) e 4 sin 2 2 θ 4 b 2 r ( 1 e 2 cos 2 θ ) ) P r + ( e 4 sin 2 2 θ 4 b 2 ( 1 e 2 cos 2 θ ) ) 2 P r 2 + ( e 2 sin 2 θ b 2 r ) 2 P r θ + ( 1 e 2 cos 2 θ b 2 r 2 ) 2 P θ 2 ] × δ ( r R ) + f ( r , θ , t ) , with initial condition P ( r , θ , 0 ) = ( r 1 ) 2 ( r L ) 2 ( θ + π ) 2 ( θ π ) 2 , and boundary conditions r P ( L , θ , t ) = r P ( 1 , θ , t ) = 0 , θ P ( r , π , t ) = θ P ( r , π , t ) = 0 , where f ( r , θ , t ) = ( r 1 ) 2 ( r L ) 2 ( θ + π ) 2 ( θ π ) 2 2 A ( t + 1 ) ( r 1 ) ( r L ) × ( 2 r 1 L ) ( θ + π ) 2 ( θ π ) 2 8 B ( t + 1 ) ( r 1 ) ( r L ) × 2 r 1 L ) ( θ + π ) ( θ π ) C ( t + 1 ) ( 12 r 2 12 r L 12 r + 8 L + 2 L 2 + 2 ) ( θ + π ) 2 ( θ π ) 2 4 E ( t + 1 ) ( r 1 ) 2 ( r L ) 2 ( 3 θ 2 π 2 ) , with A = D ( 1 e 2 cos 2 θ ) b 2 r + D [ 4 e 2 cos 2 θ ( 1 e 2 cos 2 θ ) e 4 sin 2 ( 2 θ ) ] δ ( r R ) 4 b 2 r ( 1 e 2 cos 2 θ ) , B = D e 2 sin 2 ( 2 θ ) δ ( r R ) r b 2 , C = D ( 1 e 2 cos 2 θ ) b 2 + D e 4 sin 2 ( 2 θ ) δ ( r R ) 4 b 2 ( 1 e 2 cos 2 θ ) , E = D ( 1 e 2 cos 2 θ ) δ ( r R ) b 2 r 2 . The exact solution of the equation is given by P ( r , θ , t ) = ( t + 1 ) ( r 1 ) 2 ( r L ) 2 ( θ + π ) 2 ( θ π ) 2 . The comparison between the exact and numerical solution of Eq. ( 13) subject to the initial condition and boundary conditions Eqs. ( 12)–( 15) is shown in Fig. 2. We can see that the curves in a good agreement which indicates the correctness of numerical results.

    1
    2  The comparison of the exact solution via Eq.(16) with numerical solution for b = 1, e = 0.5, D = D = 1 , τ = 10 −4, h r = 10 −2, and h θ = π/20 at T = 1.
    The probability distribution function (PDF) P( x, y, t) is a useful quantity to assess how particles are distributed in the structure. In case of classical comb with a reflected boundary conditions, the evolution of PDF P( x, y, t) can be described by the comb equation t P ( x , y , t ) = δ ( y y 0 ) D x 2 P ( x , y , t ) + D y 2 P ( x , y , t ) with the following initial and boundary conditions: P ( x , y , 0 ) = δ ( x x 0 ) δ ( y y 0 ) , P ( x , y 0 , t ) x | x ± L = 0 , P ( x , y , t ) y | y ± R = 0. The temporal evolution of marginal PDF P 1 ( x , t ) = R R P ( x , y , t ) d y , which characterizes the distribution of particle along the backbone, is detailed in Figs. 3 and 4. Figure 3 shows that, for fixed L, P 1 ( x, t) converges to the same non-Gaussian function P 1 ( x , t ) as time increases, and the speed of convergence is inversely proportional to R (the length of the lateral finger 2 R). Moreover, the function P 1 ( x , t ) is affected by L, in which the concentration reduces as the length of the backbone increases, seen in Fig. 4. This means that the particles diffuse slowly, then they saturate after a period due to limited of the backbone. Moreover, the stay time in the branch until they return to the backbone determines the nature of its diffusion and reaches a saturated state. This will be demonstrated in the next subsection. Now, we show the effects of the parameters b and e on the distribution of particles P( r, θ, t) in the whole elliptical comb-like structure, and the marginal PDF P 1 ( θ , t ) . along θ direction of the structure. The temporal evolution of the distribution in the structure is drawn from the iterative Eq. ( 9) by using Matlab. The parameters b and e reflect the geometric properties of the ellipse. The relationship among the shortest diameter b, the longest diameter a, and the eccentricity of ellipse e reads as a 2 = b 2/(1 – e 2) with 0 ≤ e < 1. Moreover, the curvature of an ellipse is given as κ = a b / ( a 2 sin 2 θ + b 2 cos 2 θ ) 3 . For fixed b, the curvature of the ellipse decreases as the value of e increases. Figure 5 depicts the distribution of particles in the structure with different parameter effects at T = 1. It shows that the concentration at the initial position increases as the shortest diameter b and the eccentricity e of ellipse increases. Namely, the diffusion of particles slows down as the length of lateral branches and the curvature of the elliptical channel increase, in which these increments affect the waiting time and obstruction in the diffusion process.

    3  Time evolution of PDF in the backbone of classical comb, along x direction, for D = D = 1 , L = 5, and (a) R = 1.5, (b) R = 6.
    4  Time evolution of PDF in the backbone of classical comb, along x direction, for D = D = 1 , R = 3, and (a) L = 5, (b) L = 10.
    5  Time evolution of probability distribution in the structure for different value of the parameter b at T = 1: (a) D = D = 1 , e = 0.5, and b = 2, (b) D = D = 1 , e = 0.5, and b = 3, (c) D = D = 1 , e = 0.5, and b = 4, (d) D = D = 1 , e = 0.5, and b = 5.

    Figures 6 and 7 depict the temporal evolution of marginal PDF along θ direction when R = 1/2, which show a symmetrical and non-Gaussian behavior according to different parameters. Figure 6 shows that the diffusion process slows down as the length of lateral branches increases, which can be attributed to the elapsed time in radial dead ends until the particle returns to the backbone. Besides, as in classical comb, it converges to a non-Gaussian function after a period that varies according to the length of the finger b. However, unlike classical comb, this function also alters according to b. The effects of the eccentricity of elliptical channel e, for fixed b, on marginal PDF along θ direction are detailed in Fig. 7. It shows that the concentration at the initial position increases as the eccentricity increases. According to PDF quantity, the results can be summarized as following: The diffusion in a finite elliptical comb-like structure is abnormal and reaches a state of saturation after a period, which varies according to the stay time in lateral dead ends. Moreover, the concentration increases with the increase of the curvature of the elliptical channel, in which the increase of the curvature obstructs particle diffusion.

    6  Time evolution of PDF in the backbone of elliptical comb, along direction θ, for D = D = 1 , e = 0.5, and (a) b = 2, (b) b = 8.
    7  Time evolution of PDF in the backbone of elliptical comb, along direction θ, for D = D = 1 , b = 4, and (a) e = 0.3, (b) e = 0.8.
    For more statistical information about this process and type of diffusion, the temporal evolution of the second moment (the mean square displacement (MSD) of the particle’s positions), which describes the speed of particle motion, in the backbone is calculated as MSD = ( x ( t ) x 0 ) 2 = ( x ( t ) x 0 ) 2 d x ( P ( x , y , t ) d y ) . Firstly, we revisit MSD along x direction in finite classical comb with respect to parameters R and L. We use numerical simulation to obtain 〈 ( x( t) – x 0) 2〉, in which we assume x 0 = 0. As shown in Fig. 8, for fixed L, 〈 x 2( t) 〉 grows non-linearly for a short time, then it converges to a fixed value. However, the speed of convergence decreases as R increases, which agrees with the results about PDF in the previous section. Besides, 〈 x 2( t) 〉 ∼ t α with α < 1 varies with respect to R and the transient sub-diffusion regimes dominate the process. Namely, one sub-diffusion regime will dominate the process as R tends to infinity. It should be noted that it has been established that x 2 ( t ) ~ t for infinite comb model. The effects of L on 〈 x 2( t) 〉 is represented in Fig. 9. It shows that the saturation threshold increases as L increases. Namely, the state of saturation will go away when L tends infinity. It should be mentioned that the saturation of the MSD is also found by considering diffusion in comb structures with stochastic resetting. [ 82, 63]

    1
    8  Temporal evolution of MSD in the backbone of classical comb, along x direction for D = D = 1 , R = {1.5, 3, 6}, and (a) L = 5, (b) L = 10.
    9  Temporal evolution of MSD in the backbone of classical comb, along x direction for D = D = 1 , L = {5, 7.5, 10}, and (a) R = 1.5, (b) R = 6.

    Now, we consider the behavior of MSD along θ direction in elliptical comb. As in the case of diffusion in classical comb, Fig. 10 shows a non-linear behavior of MSD along θ direction, then it converges to constant value and the speed of convergence decreases as b increases. However, for fixed e, it does not converge to the same value, except in case of circle comb, i.e., when e = 0. Moreover, the saturating threshold is reduced as e increases. This can be attributed to increase the perimeter of the elliptical channel, p 2 π ( a 2 + b 2 ) / 2 , and the time of stay in the fingers. For short time, Fig. 10 shows that MSD ∼ t α with 0 < α < 1, and a sub-diffusion process occurs in this direction. Besides, the value of exponent α decreases as e increases. Such transient sub-diffusion behavior was indeed demonstrated, for example, diffusion of fluorophore-labelled dextran (MW 3000) in granular layers (GL) of rat and turtle cerebella, [ 19] also diffusion in crowded media of non-inert obstacles. [ 83]

    10  Temporal evolution of MSD in the backbone of classical comb, along θ direction for D = D = 1 , b = {2, 4, 5, 6, 8, 10}, and (a) e = 0, (b) e = 0.2, (c) e = 0.5, (d) e = 0.8.
    11  MSD along θ direction for short time for D = D = 1 , e = {0, 0.3, 0.5, 0.8}, and (a) b = 4, (b) b = 7.

    The phenomenon of slow diffusion can be attributed here to the time the particles spend in the dead ends until they return to the backbone and the obstruction due to curvature of the ellipse. For more about the role of eccentricity or curvature of the ellipse, we consider the diffusion in the structure with the absence of dead ends. MSD along the elliptical channel is plotted in Fig. 12. In contrast to diffusion in elliptical comb, MSD grows linearly with time, and the diffusion is normal, see Fig. 12(a). However, MSD grows non-linearly with time as the curvature of the channel increases, see Fig. 12(b), and sub-diffusion process is dominated.

    12  MSD along elliptical channel without constrained geometries for short time for D = D = 1 , and (a) e = 0, 0.5, 0.8, and b = 10, (b) e = 0.5 and b = 3, 6, 8.

    In this paper, we have considered the features of diffusion process in an elliptical channel with dead ends and discussed the impact of geometry of the structure on this process. To manipulate this problem, a modification of comb-like model has been suggested, in which the backbone and fingers of the comb correspond to the elliptical channel and radial dead ends distributed along it, respectively. The Fokker–Planck equation that expresses the diffusion process in the structure process has been derived and solved numerically. The results show a transient sub-diffusion behavior dominates the process followed by a saturating state. The sub-diffusion regime and saturation threshold are affected by the length of the elliptical channel lateral branch and its curvature. Besides, the concentration of particles at the initial position increases with the increase in the curvature of the elliptical channel. These established results can be added to the results of previous works that refer to the effects of curvature on diffusion processes in curved channels.

    Acknowledgement Project supported by the National Natural Science Foundation of China (Grant Nos.~11772046 and 81870345).
  • Schematic picture of an elliptical comb-like structure, where the radii are continuously distributed over the ellipse of radius R = ρ( θ) with θ ∈ [0,2 π].

    The comparison of the exact solution via Eq.(16) with numerical solution for b = 1, e = 0.5, D = D = 1 , τ = 10 −4, h r = 10 −2, and h θ = π/20 at T = 1.

    Time evolution of PDF in the backbone of classical comb, along x direction, for D = D = 1 , L = 5, and (a) R = 1.5, (b) R = 6.

    Time evolution of PDF in the backbone of classical comb, along x direction, for D = D = 1 , R = 3, and (a) L = 5, (b) L = 10.

    Time evolution of probability distribution in the structure for different value of the parameter b at T = 1: (a) D = D = 1 , e = 0.5, and b = 2, (b) D = D = 1 , e = 0.5, and b = 3, (c) D = D = 1 , e = 0.5, and b = 4, (d) D = D = 1 , e = 0.5, and b = 5.

    Time evolution of PDF in the backbone of elliptical comb, along direction θ, for D = D = 1 , e = 0.5, and (a) b = 2, (b) b = 8.

    Time evolution of PDF in the backbone of elliptical comb, along direction θ, for D = D = 1 , b = 4, and (a) e = 0.3, (b) e = 0.8.

    Temporal evolution of MSD in the backbone of classical comb, along x direction for D = D = 1 , R = {1.5, 3, 6}, and (a) L = 5, (b) L = 10.

    Temporal evolution of MSD in the backbone of classical comb, along x direction for D = D = 1 , L = {5, 7.5, 10}, and (a) R = 1.5, (b) R = 6.

    10  Temporal evolution of MSD in the backbone of classical comb, along θ direction for D = D = 1 , b = {2, 4, 5, 6, 8, 10}, and (a) e = 0, (b) e = 0.2, (c) e = 0.5, (d) e = 0.8.

    11  MSD along θ direction for short time for D = D = 1 , e = {0, 0.3, 0.5, 0.8}, and (a) b = 4, (b) b = 7.

    12  MSD along elliptical channel without constrained geometries for short time for D = D = 1 , and (a) e = 0, 0.5, 0.8, and b = 10, (b) e = 0.5 and b = 3, 6, 8.

  • [1] Stein W D 1986 Transport and diffusion across cell membranes (Academic Press)
    [2] Bressloff P C 2014 Stochastic Processes in Cell Biology (Springer)
    [3] Sibatov R T, Uchaikin V V 2015 J. Comput. Phys. 293 409 10.1016/j.jcp.2015.01.022 doi: 10.1016/j.jcp.2015.01.022
    [4] Mendez V, Fedotov S, Horsthemke W 2010 Reaction-Transport Systems (Springer)
    [5] Viswanathan G M, Luz M, Raposo E P, Stanley H E 2011 The Physics of Foraging (Cambridge University Press)
    [6] Dix J A, Verkman A S 2008 Annu. Rev. Biophys. 37 247 10.1146/annurev.biophys.37.032807.125824 doi: 10.1146/annurev.biophys.37.032807.125824
    [7] Weiss M, Hashimoto H, Nilsson T 2003 Biophys. J. 84 4043 10.1016/S0006-3495(03)75130-3 doi: 10.1016/S0006-3495(03)75130-3
    [8] Höfling F, Franosch T 2013 Rep. Prog. Phys. 76 46602 10.1088/0034-4885/76/4/046602 doi: 10.1088/0034-4885/76/4/046602
    [9] Sun H G, Sheng H, Chen Y Q, Chen W, Yu Z B 2013 Chin. Phys. Lett. 30 46601 10.1088/0256-307X/30/4/046601 doi: 10.1088/0256-307X/30/4/046601
    [10] Wu W X, Zheng Z G, Song Y L, Han Y R, Sun Z C, Li C P 2020 Chin. Phys.B 29 90503 10.1088/1674-1056/ab99b7 doi: 10.1088/1674-1056/ab99b7
    [11] Wang B, Kuo J, Bae S C, Granick S 2015 Nat. Mater. 11 481 10.1038/nmat3308 doi: 10.1038/nmat3308
    [12] Metzler R, Jeon J H, Cherstvy A G 2016 Biochim. Biophys. Acta - Biomembr. 1858 2451 10.1016/j.bbamem.2016.01.022 doi: 10.1016/j.bbamem.2016.01.022
    [13] Metzler R 2020 Eur. Phys. J. Spec. Top. 229 711 10.1140/epjst/e2020-900210-x doi: 10.1140/epjst/e2020-900210-x
    [14] Metzler R 2017 Biophys. J. 112 413 10.1016/j.bpj.2016.12.019 doi: 10.1016/j.bpj.2016.12.019
    [15] Jeon J H, Javanainen M, Martinez-Seara H, Metzler R, Vattulainen I 2016 Phys. Rev.X 6 21006
    [16] Szymanski J, Weiss M 2009 Phys. Rev. Lett. 103 38102 10.1103/PhysRevLett.103.038102 doi: 10.1103/PhysRevLett.103.038102
    [17] Gupta S, Mel J, Perera R M, Zolnierczuk P, Schneider G J 2018 J. Phys. Chem. Lett. 9 2956 10.1021/acs.jpclett.8b01008 doi: 10.1021/acs.jpclett.8b01008
    [18] Tabei S, Burov S, Kim H Y, Kuznetsov A, Huynh T, Jureller J, Philipson L H, Dinner A R, Scherer N F 2013 Proc. Natl. Acad. Sci. USA 110 4911 10.1073/pnas.1221962110 doi: 10.1073/pnas.1221962110
    [19] Xiao F, Hrabe J, Hrabetova S 2015 Biophys. J. 108 2384 10.1016/j.bpj.2015.02.034 doi: 10.1016/j.bpj.2015.02.034
    [20] Goeppert N, Goldscheider N, Berkowitz B 2020 Water Res. 178 115755 10.1016/j.watres.2020.115755 doi: 10.1016/j.watres.2020.115755
    [21] Caspi A, Granek R, Elbaum M 2000 Phys. Rev. Lett. 85 5655 10.1103/PhysRevLett.85.5655 doi: 10.1103/PhysRevLett.85.5655
    [22] Reverey J F, Jeon J H, Bao H, Leippe M, Metzler R, Selhuber-Unkel C 2015 Sci. Rep. 5 11690
    [23] Rammal R, Toulouse G 1983 J. Phys. Lett. 44 13 10.1051/jphyslet:0198300440101300 doi: 10.1051/jphyslet:0198300440101300
    [24] Sokolov I M 2012 Soft Matter 8 9043 10.1039/c2sm25701g doi: 10.1039/c2sm25701g
    [25] Metzler R, Jeon J H, Cherstvy A G, Barkai E 2014 Phys. Chem. Chem. Phys. 16 24128 10.1039/C4CP03465A doi: 10.1039/C4CP03465A
    [26] Ben-Avraham D, Havlin S 2000 Diffusion and reactions in fractals and disordered systems (Cambridge university press)
    [27] Gao S L, Zhong S C, Wei K, Ma H 2012 Acta Phys. Sin. 61 6 (in Chinese)
    [28] Goodknight R C, Klikoff W A, Fatt I 1960 J. Phys. Chem. 64 1162 10.1021/j100838a014 doi: 10.1021/j100838a014
    [29] Pinner A, Nye P H 1982 J. Soil Sci. 33 25 10.1111/ejs.1982.33.issue-1 doi: 10.1111/ejs.1982.33.issue-1
    [30] Rao P, Jessup R E, Addiscott T M 1982 Soil Sci. 133 154
    [31] Tao A, Tao L, Nicholson C 2005 J. Theor. Biol. 234 525 10.1016/j.jtbi.2004.12.009 doi: 10.1016/j.jtbi.2004.12.009
    [32] Hrabe J, HrabeTová S, Segeth K 2004 Biophys. J. 87 1606 10.1529/biophysj.103.039495 doi: 10.1529/biophysj.103.039495
    [33] Ernest P, Si B R 1964 J. Gen. Physiol. 47 1129 10.1085/jgp.47.6.1129 doi: 10.1085/jgp.47.6.1129
    [34] Safford R E, Bassingthwaighte E A, Bassingthwaighte J B 1978 J. Gen. Physiol. 72 513 10.1085/jgp.72.4.513 doi: 10.1085/jgp.72.4.513
    [35] Santamaria F, Wils S, De Schutter E, Augustine G J 2006 Neuron 52 635 10.1016/j.neuron.2006.10.025 doi: 10.1016/j.neuron.2006.10.025
    [36] Iomin A, Mendéz V, Horsthemke W 2018 Fractional dynamics in comb-like structures (World Scientific)
    [37] Weiss G H, Havlin S 1986 Phys.A 134 474 10.1016/0378-4371(86)90060-9 doi: 10.1016/0378-4371(86)90060-9
    [38] Weiss G H, Havlin S 1987 Philos. Mag.B 56 941 10.1080/13642818708215329 doi: 10.1080/13642818708215329
    [39] Arkhincheev V E, Baskin E M 1991 Sov. Phys. JETP 73 161
    [40] Méndez V, Iomin A, Horsthemke W, Campos D 2017 J. Stat. Mech. Theory Exp. 2017
    [41] Méndez V, Iomin A 2013 Chaos Soliton. Fract. 53 46
    [42] Iomin A 2012 Phys. Rev.E 86 32101 10.1103/PhysRevE.86.032101 doi: 10.1103/PhysRevE.86.032101
    [43] Arkhincheev V E, Kunnen E, Baklanov M R 2011 Microelectron. Eng. 88 694 10.1016/j.mee.2010.08.028 doi: 10.1016/j.mee.2010.08.028
    [44] Liu L 2018 J. Comput. Math. 36 563 10.4208/jcm doi: 10.4208/jcm
    [45] Liu L, Zheng L, Liu F 2017 J. Mol. Liq. 233 326 10.1016/j.molliq.2017.03.034 doi: 10.1016/j.molliq.2017.03.034
    [46] Iomin A, Méndez V 2016 Chaos Soliton. Fract. 82 142
    [47] Dzhanoev A R, Sokolov I M 2018 Chaos Soliton. Fract. 106 330
    [48] Havlin S, Kiefer J E, Weiss G H 1987 Phys. Rev.A 36 1403 10.1103/PhysRevA.36.1403 doi: 10.1103/PhysRevA.36.1403
    [49] Yuste S B, Abad E, Baumgaertner A 2016 Phys. Rev.E 94 012118 10.1103/PhysRevE.94.012118 doi: 10.1103/PhysRevE.94.012118
    [50] Iomin A 2011 Phys. Rev.E 83 052106
    [51] Sandev T, Iomin A, Kantz H 2015 Phys. Rev.E 91 032108 10.1103/PhysRevE.91.032108 doi: 10.1103/PhysRevE.91.032108
    [52] Sandev T, Iomin A, Kantz H 2017 Phys. Rev.E 95 52107 10.1103/PhysRevE.95.052107 doi: 10.1103/PhysRevE.95.052107
    [53] Fan Y, Liu L, Zheng L, Li X 2019 Sci. Numer. Simul. 77 225
    [54] Wang Z, Zheng L 2020 Phys.A 549 123889
    [55] Iomin A 2021 Chaos Soliton. Fract. 142 110488
    [56] Suleiman K, Zheng L, Liu C, Zhang X, Wang E 2020 Int. Commun. Heat Mass Transf. 117 104733 10.1016/j.icheatmasstransfer.2020.104733 doi: 10.1016/j.icheatmasstransfer.2020.104733
    [57] Suleiman K, Song Q, Zhang X, Liu S, Zheng L 2022 Chaos Soliton. Fract. 155 111742
    [58] Saxton M J 2007 Biophys. J. 92 1178 10.1529/biophysj.106.092619 doi: 10.1529/biophysj.106.092619
    [59] Baskin E, Iomin A 2004 Phys. Rev. Lett. 93 120603 10.1103/PhysRevLett.93.120603 doi: 10.1103/PhysRevLett.93.120603
    [60] Sandev T, Iomin A, Kocarev L 2020 Phys. Rev.E 102 042109 10.1103/PhysRevE.102.042109 doi: 10.1103/PhysRevE.102.042109
    [61] Singh R K, Sandev T, Iomin A, Metzler R 2021 J. Phys.A 54 404006 10.1088/1751-8121/ac20ed doi: 10.1088/1751-8121/ac20ed
    [62] Wang Z, Lin P, Wang E 2021 Chaos Soliton. Fract. 148 111009
    [63] Sandev T, Domazetoski V, Iomin A, Kocarev L 2021 Mathematics 9 1 10.3390/math9010001 doi: 10.3390/math9010001
    [64] Dagdug L, Berezhkovskii A M, Makhnovskii Y A, Zitserman V Y 2007 J. Chem. Phys. 127 224712 10.1063/1.2805068 doi: 10.1063/1.2805068
    [65] Berezhkovskii A M, Dagdug L 2011 J. Chem. Phys. 134 124109 10.1063/1.3567187 doi: 10.1063/1.3567187
    [66] Ogawa N 2010 Phys. Rev.E 81 61113 10.1103/PhysRevE.81.061113 doi: 10.1103/PhysRevE.81.061113
    [67] Ogawa N 2013 Phys. Lett.A 377 2465 10.1016/j.physleta.2013.07.054 doi: 10.1016/j.physleta.2013.07.054
    [68] Chávez Y, Chacón-Acosta G, Dagdug L 2018 J. Chem. Phys. 148 214106 10.1063/1.5030892 doi: 10.1063/1.5030892
    [69] Castro-Villarreal P, Villada-Balbuena A, Méndez-Alcaraz J M, Castañeda-Priego R, Estrada-Jiménez S 2014 J. Chem. Phys. 140 214115 10.1063/1.4881060 doi: 10.1063/1.4881060
    [70] Assenza S, Mezzenga R 2018 J. Chem. Phys. 148 54902 10.1063/1.5019224 doi: 10.1063/1.5019224
    [71] Dobrowolski T, Jarmoliński A 2018 Results Phys. 8 48 10.1016/j.rinp.2017.11.020 doi: 10.1016/j.rinp.2017.11.020
    [72] Burada P S, Hnggi P, Marchesoni F, Schmid G, Talkner P 2010 Chem. Phys. Chem. 10
    [73] Bénichou O, Illien P, Oshanin G, Sarracino A, Voituriez R 2018 J. Phys. Condens. Matter 30 443001 10.1088/1361-648X/aae13a doi: 10.1088/1361-648X/aae13a
    [74] Hrabĕtová S, Hrabe J, Nicholson C 2003 J. Neurosci. 23 8351 10.1523/JNEUROSCI.23-23-08351.2003 doi: 10.1523/JNEUROSCI.23-23-08351.2003
    [75] Nicholson C, Hrabĕtová S 2017 Biophys. J. 113 2133 10.1016/j.bpj.2017.06.052 doi: 10.1016/j.bpj.2017.06.052
    [76] Tao L, Nicholson C 2004 J. Theor. Biol. 229 59 10.1016/j.jtbi.2004.03.003 doi: 10.1016/j.jtbi.2004.03.003
    [77] Sykova E, Nicholson C 2008 Physiol. Rev. 88 1277 10.1152/physrev.00027.2007 doi: 10.1152/physrev.00027.2007
    [78] Bueno-Orovio A, Kay D, Grau V, Rodriguez B, Burrage K 2014 J. R. Soc. Interface 11 20140352 10.1098/rsif.2014.0352 doi: 10.1098/rsif.2014.0352
    [79] Wolak D J, Thorne R G 2013 Mol. Pharm. 10 1492 10.1021/mp300495e doi: 10.1021/mp300495e
    [80] Vendel E, Rottschfer V, Lange E 2019 Fluids Barriers CNS 16 12 10.1186/s12987-019-0133-x doi: 10.1186/s12987-019-0133-x
    [81] Caratelli D, Natalini P, Ricci P E, Yarovoy A 2010 Appl. Math. Comput. 216 556
    [82] Domazetoski V, Masó-Puigdellosas A, Sandev T, Méndez V, Iomin A, Kocarev L 2020 Phys. Rev. Res. 2 033027 10.1103/PhysRevResearch.2.033027 doi: 10.1103/PhysRevResearch.2.033027
    [83] Ghosh S K, Cherstvy A G, Metzler R 2014 Phys. Chem. Chem. Phys. 17 1847 10.1039/C4CP03599B doi: 10.1039/C4CP03599B
  • 加载中
图(12)
计量
  • 文章访问数:  49
  • HTML全文浏览量:  50
  • PDF下载量:  0
  • 被引次数: 0
出版历程
  • 收稿日期:  2022-01-16
  • 网络出版日期:  2023-05-16
  • 刊出日期:  2023-01-01

目录

/

返回文章
返回