1
|
Kosztołowicz T. Subdiffusion equation with fractional Caputo time derivative with respect to another function in modeling transition from ordinary subdiffusion to superdiffusion. Phys Rev E 2023; 107:064103. [PMID: 37464604 DOI: 10.1103/physreve.107.064103] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/25/2022] [Accepted: 05/11/2023] [Indexed: 07/20/2023]
Abstract
We use a subdiffusion equation with fractional Caputo time derivative with respect to another function g (g-subdiffusion equation) to describe a smooth transition from ordinary subdiffusion to superdiffusion. Ordinary subdiffusion is described by the equation with the "ordinary" fractional Caputo time derivative, superdiffusion is described by the equation with a fractional Riesz-type spatial derivative. We find the function g for which the solution (Green's function, GF) to the g-subdiffusion equation takes the form of GF for ordinary subdiffusion in the limit of small time and GF for superdiffusion in the limit of long time. To solve the g-subdiffusion equation we use the g-Laplace transform method. It is shown that the scaling properties of the GF for g-subdiffusion and the GF for superdiffusion are the same in the long time limit. We conclude that for a sufficiently long time the g-subdiffusion equation describes superdiffusion well, despite a different stochastic interpretation of the processes. Then, paradoxically, a subdiffusion equation with a fractional time derivative describes superdiffusion. The superdiffusive effect is achieved here not by making anomalously long jumps by a diffusing particle, but by greatly increasing the particle jump frequency which is derived by means of the g-continuous-time random walk model. The g-subdiffusion equation is shown to be quite general, it can be used in modeling of processes in which a kind of diffusion change continuously over time. In addition, some methods used in modeling of ordinary subdiffusion processes, such as the derivation of local boundary conditions at a thin partially permeable membrane, can be used to model g-subdiffusion processes, even if this process is interpreted as superdiffusion.
Collapse
Affiliation(s)
- Tadeusz Kosztołowicz
- Institute of Physics, Jan Kochanowski University, Uniwersytecka 7, 25-406 Kielce, Poland
| |
Collapse
|
2
|
Angelani L, Garra R. g-fractional diffusion models in bounded domains. Phys Rev E 2023; 107:014127. [PMID: 36797940 DOI: 10.1103/physreve.107.014127] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2022] [Accepted: 01/04/2023] [Indexed: 01/19/2023]
Abstract
In the recent literature, the g-subdiffusion equation involving Caputo fractional derivatives with respect to another function has been studied in relation to anomalous diffusions with a continuous transition between different subdiffusive regimes. In this paper we study the problem of g-fractional diffusion in a bounded domain with absorbing boundaries. We find the explicit solution for the initial boundary value problem, and we study the first-passage time distribution and the mean first-passage time (MFPT). The main outcome is the proof that with a particular choice of the function g it is possible to obtain a finite MFPT, differently from the anomalous diffusion described by a fractional heat equation involving the classical Caputo derivative.
Collapse
Affiliation(s)
- L Angelani
- ISC-CNR, Institute for Complex Systems, Piazzale A. Moro 2, 00185 Rome, Italy.,Dipartimento di Fisica, Sapienza Università di Roma, Piazzale A. Moro 2, 00185 Rome, Italy
| | - R Garra
- Institute of Marine Sciences, National Research Council (CNR), Via del Fosso del Cavaliere, I-00133 Rome, Italy
| |
Collapse
|
3
|
Kosztołowicz T, Dutkiewicz A. Composite subdiffusion equation that describes transient subdiffusion. Phys Rev E 2022; 106:044119. [PMID: 36397481 DOI: 10.1103/physreve.106.044119] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/25/2022] [Accepted: 09/27/2022] [Indexed: 06/16/2023]
Abstract
A composite subdiffusion equation with fractional Caputo time derivative with respect to another function g is used to describe a process of a continuous transition from subdiffusion with parameters α and D_{α} to subdiffusion with parameters β and D_{β}. The parameters are defined by the time evolution of the mean square displacement of diffusing particle σ^{2}(t)=2D_{i}t^{i}/Γ(1+i), i=α,β. The function g controls the process at intermediate times. The composite subdiffusion equation is more general than the ordinary fractional subdiffusion equation with constant parameters; it has potentially wide application in modeling diffusion processes with changing parameters.
Collapse
Affiliation(s)
- Tadeusz Kosztołowicz
- Institute of Physics, Jan Kochanowski University, Uniwersytecka 7, 25-406 Kielce, Poland
| | - Aldona Dutkiewicz
- Faculty of Mathematics and Computer Science, Adam Mickiewicz University, Uniwersytetu Poznańskiego 4, 61-614 Poznań, Poland
| |
Collapse
|
4
|
Kosztołowicz T, Dutkiewicz A, Lewandowska KD, Wąsik S, Arabski M. Subdiffusion equation with Caputo fractional derivative with respect to another function in modeling diffusion in a complex system consisting of a matrix and channels. Phys Rev E 2022; 106:044138. [PMID: 36397549 DOI: 10.1103/physreve.106.044138] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/26/2021] [Accepted: 10/07/2022] [Indexed: 06/16/2023]
Abstract
Anomalous diffusion of an antibiotic (colistin) in a system consisting of packed gel (alginate) beads immersed in water is studied experimentally and theoretically. The experimental studies are performed using the interferometric method of measuring concentration profiles of a diffusing substance. We use the g-subdiffusion equation with the fractional Caputo time derivative with respect to another function g to describe the process. The function g and relevant parameters define anomalous diffusion. We show that experimentally measured time evolution of the amount of antibiotic released from the system determines the function g. The process can be interpreted as subdiffusion in which the subdiffusion parameter (exponent) α decreases over time. The g-subdiffusion equation, which is more general than the "ordinary" fractional subdiffusion equation, can be widely used in various fields of science to model diffusion in a system in which parameters, and even a type of diffusion, evolve over time. We postulate that diffusion in a system composed of channels and a matrix can be described by the g-subdiffusion equation, just like diffusion in a system of packed gel beads placed in water.
Collapse
Affiliation(s)
- Tadeusz Kosztołowicz
- Institute of Physics, Jan Kochanowski University, Uniwersytecka 7, 25-406 Kielce, Poland
| | - Aldona Dutkiewicz
- Faculty of Mathematics and Computer Science, Adam Mickiewicz University, Uniwersytetu Poznańskiego 4, 61-614 Poznań, Poland
| | - Katarzyna D Lewandowska
- Department of Radiological Informatics and Statistics, Medical University of Gdańsk, Tuwima 15, 80-210 Gdańsk, Poland
| | - Sławomir Wąsik
- Department of Medical Physics and Biophysics, Institute of Physics, Jan Kochanowski University, Uniwersytecka 7, 25-406 Kielce, Poland
| | - Michał Arabski
- Institute of Biology, Jan Kochanowski University, Uniwersytecka 7, 25-406 Kielce, Poland
| |
Collapse
|
5
|
Wang X, Chen Y. Random diffusivity processes in an external force field. Phys Rev E 2022; 106:024112. [PMID: 36109990 DOI: 10.1103/physreve.106.024112] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/03/2022] [Accepted: 07/21/2022] [Indexed: 06/15/2023]
Abstract
Brownian yet non-Gaussian processes have recently been observed in numerous biological systems, and corresponding theories have been constructed based on random diffusivity models. Considering the particularity of random diffusivity, this paper studies the effect of an external force acting on two kinds of random diffusivity models whose difference is embodied in whether the fluctuation-dissipation theorem is valid. Based on the two random diffusivity models, we derive the Fokker-Planck equations with an arbitrary external force, and we analyze various observables in the case with a constant force, including the Einstein relation, the moments, the kurtosis, and the asymptotic behaviors of the probability density function of particle displacement at different timescales. Both the theoretical results and numerical simulations of these observables show a significant difference between the two kinds of random diffusivity models, which implies the important role of the fluctuation-dissipation theorem in random diffusivity systems.
Collapse
Affiliation(s)
- Xudong Wang
- School of Mathematics and Statistics, Nanjing University of Science and Technology, Nanjing 210094, People's Republic of China
| | - Yao Chen
- College of Sciences, Nanjing Agricultural University, Nanjing 210095, People's Republic of China
| |
Collapse
|
6
|
Kosztołowicz T. First-passage time for the g-subdiffusion process of vanishing particles. Phys Rev E 2022; 106:L022104. [PMID: 36110021 DOI: 10.1103/physreve.106.l022104] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2022] [Accepted: 08/05/2022] [Indexed: 06/15/2023]
Abstract
Subdiffusion equation and molecule survival equation, both with Caputo fractional time derivatives with respect to other functions g_{1} and g_{2}, respectively, are used to describe diffusion of a molecule that can disappear at any time with a constant probability. The process can be interpreted as an "ordinary" subdiffusion and "ordinary" molecule survival process in which timescales are changed by the functions g_{1} and g_{2}. We derive the first-passage time distribution for the process. The mutual influence of subdiffusion and molecule-vanishing processes can be included in the model when the functions g_{1} and g_{2} are related to each other. As an example, we consider the processes in which subdiffusion and molecule survival are highly related, which corresponds to the case of g_{1}≡g_{2}.
Collapse
Affiliation(s)
- Tadeusz Kosztołowicz
- Institute of Physics, Jan Kochanowski University, Uniwersytecka 7, 25-406 Kielce, Poland
| |
Collapse
|
7
|
Abstract
The internal energy (U-energy) conversion to free energy (F-energy) and energy dissipation (S-energy) is a basic process that enables the continuity of life on Earth. Here, we present a novel method of evaluating F-energy in a membrane system containing ternary solutions of non-electrolytes based on the Kr version of the Kedem–Katchalsky–Peusner (K–K–P) formalism for concentration polarization conditions. The use of this formalism allows the determination of F-energy based on the production of S-energy and coefficient of the energy conversion efficiency. The K–K–P formalism requires the calculation of the Peusner coefficients Kijr and Kdetr (i, j ∈ {1, 2, 3}, r = A, B), which are necessary to calculate S-energy, the degree of coupling and coefficients of energy conversion efficiency. In turn, the equations for S-energy and coefficients of energy conversion efficiency are used in the F-energy calculations. The Kr form of the Kedem–Katchalsky–Peusner model equations, containing the Peusner coefficients Kijr and Kdetr, enables the analysis of energy conversion in membrane systems and is a useful tool for studying the transport properties of membranes. We showed that osmotic pressure dependences of indicated Peusner coefficients, energy conversion efficiency coefficient, entropy and energy production are nonlinear. These nonlinearities were caused by pseudophase transitions from non-convective to convective states or vice versa. The method presented in the paper can be used to assess F-energy resources. The results can be adapted to various membrane systems used in chemical engineering, environmental engineering or medical applications. It can be used in designing new technologies as a part of process management.
Collapse
|
8
|
Kosztołowicz T, Dutkiewicz A. Stochastic interpretation of g-subdiffusion process. Phys Rev E 2021; 104:L042101. [PMID: 34781426 DOI: 10.1103/physreve.104.l042101] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/23/2021] [Accepted: 09/23/2021] [Indexed: 06/13/2023]
Abstract
Recently, we considered the g-subdiffusion equation with a fractional Caputo time derivative with respect to another function g, T. Kosztołowicz et al. [Phys. Rev. E 104, 014118 (2021)2470-004510.1103/PhysRevE.104.014118]. This equation offers different possibilities for modeling diffusion such as a process in which a type of diffusion evolves continuously over time. However, the equation has not been derived from a stochastic model and the stochastic interpretation of g subdiffusion is still unknown. In this Letter, we show the stochastic foundations of this process. We derive the equation by means of a modified continuous time random walk model. An interpretation of the g-subdiffusion process is also discussed.
Collapse
Affiliation(s)
- Tadeusz Kosztołowicz
- Institute of Physics, Jan Kochanowski University, Uniwersytecka 7, 25-406 Kielce, Poland
| | - Aldona Dutkiewicz
- Faculty of Mathematics and Computer Science, Adam Mickiewicz University, Uniwersytetu Poznańskiego 4, 61-614 Poznań, Poland
| |
Collapse
|