1
|
Rajyaguru A, Metzler R, Cherstvy AG, Berkowitz B. Quantifying anomalous chemical diffusion through disordered porous rock materials. Phys Chem Chem Phys 2025. [PMID: 40223812 DOI: 10.1039/d5cp00654f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/15/2025]
Abstract
Fickian (normal) diffusion models show limitations in quantifying diffusion-controlled migration of solute species through porous rock structures, as observed in experiments. Anomalous diffusion prevails and can be interpreted using a Continuous Time Random Walk (CTRW) framework with a clear mechanistic underpinning. From the associated fractional diffusion equation we derive solutions over a broad range of anomalous diffusion behaviours, from highly anomalous to nearly Fickian, that yield temporal breakthrough curves and spatial concentration profiles of diffusing solutes. We illustrate that these solutions can be tailored to match realistic experimental conditions and resulting measurements that display anomalous diffusion. In particular, our analysis enables clear differentiation between early-time Fickian and anomalous diffusion, which becomes more pronounced over longer durations. It is shown that recent measurements of diffusion in natural rocks display distinct anomalous behaviour, with significant implications for critical assessment of solute migration in diverse geological and engineering applications.
Collapse
Affiliation(s)
- Ashish Rajyaguru
- Paul Scherrer Institut, 5232 Villigen, Switzerland.
- Department of Earth and Planetary Sciences, Weizmann Institute of Science, Rehovot 7610001, Israel
| | - Ralf Metzler
- Institute for Physics & Astronomy, University of Potsdam, 14476 Potsdam, Germany
- Asia Pacific Centre for Theoretical Physics, Pohang 37673, Republic of Korea
| | - Andrey G Cherstvy
- Institute for Physics & Astronomy, University of Potsdam, 14476 Potsdam, Germany
| | - Brian Berkowitz
- Department of Earth and Planetary Sciences, Weizmann Institute of Science, Rehovot 7610001, Israel
| |
Collapse
|
2
|
Adejumobi OP, Mantsevich VN, Palyulin VV. Diffusion of fast and slow excitons with an exchange in quasi-two-dimensional systems. Phys Rev E 2024; 110:054139. [PMID: 39690650 DOI: 10.1103/physreve.110.054139] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/10/2024] [Accepted: 10/24/2024] [Indexed: 12/19/2024]
Abstract
By means of analytical calculations and numerical simulations, we study the diffusion properties in quasi-two-dimensional structures with two exciton subsystems with an exchange between them. The experimental realization is possible in systems where fast and slow exciton subsystems appear. For substantially different diffusion coefficients of the species, the negative diffusion can be observed if one measures the transport properties of only a single subsystem, just as was obtained in experimental studies for quasi-two-dimensional semiconductor systems. The initial transition from a fast subsystem to a slow one results in a delayed release of fast excitons in the area close to the original excitation spot. Hence, the signal from the fast excitons alone includes the delayed replenishment from the transition of slow quasiparticles. This is seen as the narrowing of the exciton density profile or decrease of mean-squared displacement which is then interpreted as a negative diffusion. We show that the analytical theory matches the available experimental data for negative diffusion quite well. The average diffusion coefficients for the combined population are analytically expressed through the diffusion coefficients of fast and slow excitons. Simple analytical expressions for effective diffusion coefficients obtained in limiting cases are of interest both for theoretical and experimental analysis of not only the exciton transport, but also for a variety of systems, where fast and slow moving subsystems are present.
Collapse
|
3
|
Pfannstiel A, Imlau M, Bazzan M, Vittadello L. Defect related anomalous mobility of small polarons in dielectric oxides at the example of congruent lithium niobate. Sci Rep 2024; 14:22047. [PMID: 39333554 PMCID: PMC11436959 DOI: 10.1038/s41598-024-71858-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/17/2024] [Accepted: 09/01/2024] [Indexed: 09/29/2024] Open
Abstract
Polarons play a major role in the description of optical, electrical and dielectrical properties of several ferroelectric oxides. The motion of those particles occurs by elementary hops among the material lattice sites. In order to compute macroscopic transport parameters such as charge mobility, normal (i.e. Fickian) diffusion laws are generally assumed. In this paper we show that when defect states able to trap the polarons for long times are considered, significant deviations from the normal diffusion behaviour arise. As an example of this behavior, we consider here the case of lithium niobate (LN). This can be considered as a prototypical system, having a rich landscape of interacting polaron types and for which a significant wealth of information is available in literature. Our analysis considers the case of a stoichiometric, defect-free lithium niobate containing a certain concentration of small electron polarons hopping on regular Nb sites, and compares it to the material in congruent composition, which is generally found in real-life applications and which is characterized by a large concentration of antisite NbLi defects. While in the first case the charge carriers are free polarons hopping on a regular Nb sublattice, in the second case a fraction of polarons is trapped on antisite defects. Thus, in the congruent material, a range of different hopping possibilities arises, depending on the type of starting and destination sites. We develop a formalism encompassing all these microscopic processes in the framework of a switching diffusion model which can be well approximated by a mobile-immobile transport model providing explicit expressions for the polaron mobility. Finally, starting from the Marcus-Holstein's model for the polaron hopping frequency we verify by means of a Monte Carlo approach the diffusion/mobility of the different polarons species showing that, while free polarons obey the laws for normal diffusion as expected, bound polarons follow an anomalous diffusion behaviour and that in the case of the congruent crystal where mixed free and bound polaron transport is involved, our expressions indeed provide a satisfactory description.
Collapse
Affiliation(s)
- Anton Pfannstiel
- Department of Mathematics/Informatics/Physics, Institute of Physics, University of Osnabrück, Barbarastraße 7, 49076, Osnabrück, Germany
| | - Mirco Imlau
- Department of Mathematics/Informatics/Physics, Institute of Physics, University of Osnabrück, Barbarastraße 7, 49076, Osnabrück, Germany
- Center for Cellular Nanoanalytics, University of Osnabrück, Barbarastraße 11, 49076, Osnabrück, Germany
| | - Marco Bazzan
- Dipartimento di Fisica e Astronomia, Università di Padova, Via Marzolo 8, 35131, Padova, Italy
| | - Laura Vittadello
- Department of Mathematics/Informatics/Physics, Institute of Physics, University of Osnabrück, Barbarastraße 7, 49076, Osnabrück, Germany.
- Center for Cellular Nanoanalytics, University of Osnabrück, Barbarastraße 11, 49076, Osnabrück, Germany.
| |
Collapse
|
4
|
Federbush A, Moscovich A, Bar-Sinai Y. Hidden Markov modeling of single-particle diffusion with stochastic tethering. Phys Rev E 2024; 109:034129. [PMID: 38632757 DOI: 10.1103/physreve.109.034129] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/03/2023] [Accepted: 02/14/2024] [Indexed: 04/19/2024]
Abstract
The statistics of the diffusive motion of particles often serve as an experimental proxy for their interaction with the environment. However, inferring the physical properties from the observed trajectories is challenging. Inspired by a recent experiment, here we analyze the problem of particles undergoing two-dimensional Brownian motion with transient tethering to the surface. We model the problem as a hidden Markov model where the physical position is observed and the tethering state is hidden. We develop an alternating maximization algorithm to infer the hidden state of the particle and estimate the physical parameters of the system. The crux of our method is a saddle-point-like approximation, which involves finding the most likely sequence of hidden states and estimating the physical parameters from it. Extensive numerical tests demonstrate that our algorithm reliably finds the model parameters and is insensitive to the initial guess. We discuss the different regimes of physical parameters and the algorithm's performance in these regimes. We also provide a free software implementation of our algorithm.
Collapse
Affiliation(s)
- Amit Federbush
- Department of Condensed Matter Physics, Tel Aviv University, Tel Aviv 69978, Israel
- Center for Physics and Chemistry of Living Systems, Tel Aviv University, Tel Aviv 69978, Israel
| | - Amit Moscovich
- Department of Statistics and Operations Research, Tel Aviv University, Tel Aviv 69978, Israel
| | - Yohai Bar-Sinai
- Department of Condensed Matter Physics, Tel Aviv University, Tel Aviv 69978, Israel
- Center for Physics and Chemistry of Living Systems, Tel Aviv University, Tel Aviv 69978, Israel
| |
Collapse
|
5
|
Kurilovich AA, Mantsevich VN, Chechkin AV, Palyulin VV. Negative diffusion of excitons in quasi-two-dimensional systems. Phys Chem Chem Phys 2024; 26:922-935. [PMID: 38088027 DOI: 10.1039/d3cp03521b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2024]
Abstract
We show how two different mobile-immobile type models explain the observation of negative diffusion of excitons reported in experimental studies in quasi-two-dimensional semiconductor systems. The main reason for the effect is the initial trapping and a delayed release of free excitons in the area close to the original excitation spot. The density of trapped excitons is not registered experimentally. Hence, the signal from the free excitons alone includes the delayed release of not diffusing trapped particles. This is seen as the narrowing of the exciton density profile or decrease of mean-squared displacement which is then interpreted as a negative diffusion. The effect is enhanced with the increase of recombination intensity as well as the rate of the exciton-exciton binary interactions.
Collapse
Affiliation(s)
- Aleksandr A Kurilovich
- Center for Energy Science and Technology, Skolkovo Institute of Science and Technology, Bolshoy Boulevard 30, 121205, Moscow, Russia
| | - Vladimir N Mantsevich
- Chair of Semiconductors and Cryoelectronics, Physics department, Lomonosov Moscow State University, 119991, Moscow, Russia
| | - Aleksei V Chechkin
- Faculty of Pure and Applied Mathematics, Hugo Steinhaus Center, Wroclaw University of Science and Technology, Wyspianskiego 27, 50-370 Wroclaw, Poland
- Institute for Physics & Astronomy, University of Potsdam, D-14476 Potsdam-Golm, Germany
- Akhiezer Institute for Theoretical Physics National Science Center "Kharkov Institute of Physics and Technology", 61108, Kharkov, Ukraine
| | - Vladimir V Palyulin
- Applied AI centre, Skolkovo Institute of Science and Technology, Bolshoy Boulevard 30, Moscow, 121205, Russia.
| |
Collapse
|
6
|
Seckler H, Szwabiński J, Metzler R. Machine-Learning Solutions for the Analysis of Single-Particle Diffusion Trajectories. J Phys Chem Lett 2023; 14:7910-7923. [PMID: 37646323 DOI: 10.1021/acs.jpclett.3c01351] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/01/2023]
Abstract
Single-particle traces of the diffusive motion of molecules, cells, or animals are by now routinely measured, similar to stochastic records of stock prices or weather data. Deciphering the stochastic mechanism behind the recorded dynamics is vital in understanding the observed systems. Typically, the task is to decipher the exact type of diffusion and/or to determine the system parameters. The tools used in this endeavor are currently being revolutionized by modern machine-learning techniques. In this Perspective we provide an overview of recently introduced methods in machine-learning for diffusive time series, most notably, those successfully competing in the anomalous diffusion challenge. As such methods are often criticized for their lack of interpretability, we focus on means to include uncertainty estimates and feature-based approaches, both improving interpretability and providing concrete insight into the learning process of the machine. We expand the discussion by examining predictions on different out-of-distribution data. We also comment on expected future developments.
Collapse
Affiliation(s)
- Henrik Seckler
- Institute of Physics & Astronomy, University of Potsdam, 14476 Potsdam-Golm, Germany
| | - Janusz Szwabiński
- Hugo Steinhaus Center, Faculty of Pure and Applied Mathematics, Wrocław University of Science and Technology, Wybrzeże Wyspiańskiego 27, 50-370 Wrocław, Poland
| | - Ralf Metzler
- Institute of Physics & Astronomy, University of Potsdam, 14476 Potsdam-Golm, Germany
- Asia Pacific Center for Theoretical Physics, Pohang 37673, Republic of Korea
| |
Collapse
|
7
|
Liang Y, Wang W, Metzler R. Anomalous diffusion, non-Gaussianity, and nonergodicity for subordinated fractional Brownian motion with a drift. Phys Rev E 2023; 108:024143. [PMID: 37723819 DOI: 10.1103/physreve.108.024143] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/20/2023] [Accepted: 08/11/2023] [Indexed: 09/20/2023]
Abstract
The stochastic motion of a particle with long-range correlated increments (the moving phase) which is intermittently interrupted by immobilizations (the trapping phase) in a disordered medium is considered in the presence of an external drift. In particular, we consider trapping events whose times follow a scale-free distribution with diverging mean trapping time. We construct this process in terms of fractional Brownian motion with constant forcing in which the trapping effect is introduced by the subordination technique, connecting "operational time" with observable "real time." We derive the statistical properties of this process such as non-Gaussianity and nonergodicity, for both ensemble and single-trajectory (time) averages. We demonstrate nice agreement with extensive simulations for the probability density function, skewness, kurtosis, as well as ensemble and time-averaged mean-squared displacements. We place a specific emphasis on the comparisons between the cases with and without drift.
Collapse
Affiliation(s)
- Yingjie Liang
- College of Mechanics and Materials, Hohai University, 211100 Nanjing, China
- University of Potsdam, Institute of Physics and Astronomy, 14476 Potsdam-Golm, Germany
| | - Wei Wang
- University of Potsdam, Institute of Physics and Astronomy, 14476 Potsdam-Golm, Germany
| | - Ralf Metzler
- University of Potsdam, Institute of Physics and Astronomy, 14476 Potsdam-Golm, Germany
- Asia Pacific Centre for Theoretical Physics, Pohang 37673, Republic of Korea
| |
Collapse
|