File dynamics
Encyclopedia
In chemistry
, physics
, mathematics
and related fields, file dynamics (sometimes called, single file dynamics) is the diffusion of N (N → ∞) identical Brownian hard spheres
in a quasi-one-dimensional channel of length L (L → ∞), such that the spheres do not jump one on top of the other, and the average particle's density is approximately fixed. The most well known statistical properties of this process is that the mean square displacement (MSD
) of a particle in the file follows,
, and its probability density function
(PDF) is Gaussian in poistion with a variance MSD.
Results in files that generalize the basic file include:
Generalizations of the basic file are important since these models represent reality much more accurately than the basic file. Indeed, file dynamics are used in modeling numerous microscopic processes: the diffusion within biological and synthetic pores and porous material, the diffusion along 1D objects, such as in biological roads, the dynamics of a monomer in a polymer, etc.
, the joint probability density function
(PDF) for all the particles in file, obeys a normal diffusion equation:
In
, 
is the set of particles' positions at time 
and 
is the set of the particles' initial positions at the initial time 
(set to zero). Equation (1) is solved with the appropriate boundary conditions, which reflect the hard-sphere nature of the file:
and with the appropriate initial condition:
In a simple file, the initial density is fixed, namely,
, where 
is a parameter that represents a microscopic length. The PDFs' coordinates must obey the order: 
.
with the boundary conditions:
and with the initial condition, Eq. , where the particles’ initial positions obey:
The file diffusion coefficients are taken independently from the PDF,
where Λ has a finite value that represents the fastest diffusion coefficient in the file.
, where k is a parameter. Then, all the particles in the file stand still for this random period, where afterwards, all the particles attempt jumping in accordance with the rules of the file. This procedure is carried on over and over again. The equation of motion for the particles’ PDF in a renewal-anomalous file is obtained when convoluting the equation of motion for a Brownian file with a kernel 
:
Here, the kernel
and the WT-PDF 
are related in Laplace space, 
. (The Laplace transform of a function 
reads, 
.) The reflecting boundary conditions accompanied Eq. are obtained when convoluting the boundary conditions of a Brownian file with the kernel 
, where here and in a Brownian file the initial conditions are identical.
(
is the same for all the particles), the anomalous file is not a renewal file. The basic dynamical cycle in such a file consists of the following steps: a particle with the fastest jumping time in the file, say, 
for particle i, attempts a jump. Then, the waiting times for all the other particles are adjusted: we subtract 
from each of them. Finally, a new waiting time is drawn for particle i. The most crucial difference among renewal anomalous files and anomalous files that are not renewal is that when each particle has its own clock, the particles are in fact connected also in the time domain, and the outcome is further slowness in the system (proved in the main text). The equation of motion for the PDF in anomalous files of independent particles reads:
Note that the time argument in the PDF
is a vector of times: 
, and 
. Adding all the coordinates and performing the integration in the order of faster times first (the order is determined randomly from a uniform distribution in the space of configurations) gives the full equation of motion in anomalous files of independent particles (averaging of the equation over all configurations is therefore further required). Indeed, even Eq. is very complicated, and averaging further complicates things.
Here, the index
goes on all the permutations of the initial coordinates, and contains 
permutations. From Eq. , the PDF of a tagged particle in the file, 
, is calculated
In Eq. ,
, 
(
is the initial condition of the tagged particle), and 
. The MSD for the tagged particle is obtained directly from Eq. :
Starting from Eq. , the PDF of the tagged particle in the heterogeneous file follows,
The MSD of a tagged particle in a heterogeneous file is taken from Eq. :
(The subscript nrml stands for normal dynamics.) From Eq. , it is straightforward relating the MSD of Brownian heterogeneous files and renewal-anomalous heterogeneous files,
From Eq. , one finds that the MSD of a file with normal dynamics in the power of
is the MSD of the corresponding renewal-anomalous file,
Here,
is the number of particles in the covered-length 
, and 
is the MAD of a free anomalous particle, 
. In Eq. , 
enters the calculations since all the particles within the distance 
from the tagged one must move in the same direction in order that the tagged particle will reach a distance 
from its initial position. Based on Eq. , we write a generalized scaling law for anomalous files of independent particles:
The first term on the right hand side of Eq. appears also in renewal files; yet, the term f(n) is unique. f(n) is the probability that accounts for the fact that for moving n anomalous independent particles in the same direction, when these particles indeed try jumping in the same direction (expressed with the term, (
), the particles in the periphery must move first so that the particles in the middle of the file will have the free space for moving, demanding faster jumping times for those in the periphery. f(n) appears since there is not a typical timescale for a jump in anomalous files, and the particles are independent, and so a particular particle can stand still for a very long time, substantially limiting the options of progress for the particles around him, during this time. Clearly,
, where f(n) = 1 for renewal files since the particles jump together, yet also in files of independent particles with 
, since in such files there is a typical timescale for a jump, considered the time for a synchronized jump. We calculate f(n) from the number of configurations in which the order of the particles’ jumping times enables motion; that is, an order where the faster particles are always located towards the periphery. For n particles, there are n! different configurations, where one configuration is the optimal one; so, 
. Yet, although not optimal, propagation is also possible in many other configurations; when m is the number of particles that move, then,
where
counts the number of configurations in which those m particles around the tagged one have the optimal jumping order. Now, even when m~n/2, 
. Using in Eq. , 
(
a small number larger than 1), we see,
(In Eq. , we use,
.) Equation shows that asymptotically the particles are extremely slow in anomalous files of independent particles.
, follows,
In Figure 1 we show trajectories from 9 particles in a file of 501 particles. (It is recommended opening the file in a new window). The upper panels show trajetcories for
and the lower panels show trajectories for 
. For each value of 
shown are trajectories in the early stages of the simulations (left) and in all stages of the simulation (right). The panels exhibit the phenomenon of the clustering, where the trajectories attract each other and then move pretty much together.
Chemistry
Chemistry is the science of matter, especially its chemical reactions, but also its composition, structure and properties. Chemistry is concerned with atoms and their interactions with other atoms, and particularly with the properties of chemical bonds....
, physics
Physics
Physics is a natural science that involves the study of matter and its motion through spacetime, along with related concepts such as energy and force. More broadly, it is the general analysis of nature, conducted in order to understand how the universe behaves.Physics is one of the oldest academic...
, mathematics
Mathematics
Mathematics is the study of quantity, space, structure, and change. Mathematicians seek out patterns and formulate new conjectures. Mathematicians resolve the truth or falsity of conjectures by mathematical proofs, which are arguments sufficient to convince other mathematicians of their validity...
and related fields, file dynamics (sometimes called, single file dynamics) is the diffusion of N (N → ∞) identical Brownian hard spheres
Brownian motion
Brownian motion or pedesis is the presumably random drifting of particles suspended in a fluid or the mathematical model used to describe such random movements, which is often called a particle theory.The mathematical model of Brownian motion has several real-world applications...
in a quasi-one-dimensional channel of length L (L → ∞), such that the spheres do not jump one on top of the other, and the average particle's density is approximately fixed. The most well known statistical properties of this process is that the mean square displacement (MSD
MSD
MSD may refer to:*Macromolecular Structure Database, one of the services provided by the European Bioinformatics Institute*Mahendra Singh Dhoni, Indian cricketer, wicketkeeper batsman*Marine Sanitation Device*Maryland School for the Deaf...
) of a particle in the file follows,
, and its probability density functionProbability density function
In probability theory, a probability density function , or density of a continuous random variable is a function that describes the relative likelihood for this random variable to occur at a given point. The probability for the random variable to fall within a particular region is given by the...
(PDF) is Gaussian in poistion with a variance MSD.
Results in files that generalize the basic file include:
- In files with a density law that is not fixed, but decays as a power law with an exponent a with the distance from the origin, the particle in the origin has a MSD that scales like,
, with a Gaussian PDF. - When, in addition, the particles' diffusion coefficients are distributed like a power law with exponent γ (around the origin), the MSD follows,
, with a Gaussian PDF. - In anomalous files that are renewal, namely, when all particles attempt a jump together, yet, with jumping times taken from a distribution that decays as a power law with an exponent, −1 − α, the MSD scales like the MSD of the corresponding normal file, in the power of α.
- In anomalous files of independent particles, the MSD is very slow and scales like,
. Even more exciting, the particles form clusters in such files, defining a dynamical phase transition. This depends on the anomaly power α: the percentage of particles in clusters ξ follows, 
. - Other generalizations include: when the particles can bypass each other with a constant probability upon encounter, an enhanced diffusion is seen. When the particles interact with the channel, a slower diffusion is observed. Files in embedded in two-dimensions show similar characteristics of files in one dimension.
Generalizations of the basic file are important since these models represent reality much more accurately than the basic file. Indeed, file dynamics are used in modeling numerous microscopic processes: the diffusion within biological and synthetic pores and porous material, the diffusion along 1D objects, such as in biological roads, the dynamics of a monomer in a polymer, etc.
Simple files
In simple Brownian files,, the joint probability density functionProbability density function
In probability theory, a probability density function , or density of a continuous random variable is a function that describes the relative likelihood for this random variable to occur at a given point. The probability for the random variable to fall within a particular region is given by the...
(PDF) for all the particles in file, obeys a normal diffusion equation:
In
, is the set of particles' positions at time and is the set of the particles' initial positions at the initial time (set to zero). Equation (1) is solved with the appropriate boundary conditions, which reflect the hard-sphere nature of the file:and with the appropriate initial condition:
In a simple file, the initial density is fixed, namely,
, where is a parameter that represents a microscopic length. The PDFs' coordinates must obey the order: .Heterogeneous files
In such files, the equation of motion follows,with the boundary conditions:
and with the initial condition, Eq. , where the particles’ initial positions obey:
The file diffusion coefficients are taken independently from the PDF,
where Λ has a finite value that represents the fastest diffusion coefficient in the file.
Renewal, anomalous, heterogeneous files
In renewal-anomalous files, a random period is taken independently from a WT-PDF of the form:, where k is a parameter. Then, all the particles in the file stand still for this random period, where afterwards, all the particles attempt jumping in accordance with the rules of the file. This procedure is carried on over and over again. The equation of motion for the particles’ PDF in a renewal-anomalous file is obtained when convoluting the equation of motion for a Brownian file with a kernel :Here, the kernel
and the WT-PDF are related in Laplace space, . (The Laplace transform of a function reads, .) The reflecting boundary conditions accompanied Eq. are obtained when convoluting the boundary conditions of a Brownian file with the kernel , where here and in a Brownian file the initial conditions are identical.Anomalous files with independent particles
When each particle in the anomalous file is assigned with its own jumping time drawn form ( is the same for all the particles), the anomalous file is not a renewal file. The basic dynamical cycle in such a file consists of the following steps: a particle with the fastest jumping time in the file, say, for particle i, attempts a jump. Then, the waiting times for all the other particles are adjusted: we subtract from each of them. Finally, a new waiting time is drawn for particle i. The most crucial difference among renewal anomalous files and anomalous files that are not renewal is that when each particle has its own clock, the particles are in fact connected also in the time domain, and the outcome is further slowness in the system (proved in the main text). The equation of motion for the PDF in anomalous files of independent particles reads:Note that the time argument in the PDF
is a vector of times: , and . Adding all the coordinates and performing the integration in the order of faster times first (the order is determined randomly from a uniform distribution in the space of configurations) gives the full equation of motion in anomalous files of independent particles (averaging of the equation over all configurations is therefore further required). Indeed, even Eq. is very complicated, and averaging further complicates things.Simple files
The solution of Eqs. - is a complete set of permutations of all initial coordinates appearing in the Gaussians,Here, the index
goes on all the permutations of the initial coordinates, and contains permutations. From Eq. , the PDF of a tagged particle in the file, , is calculated In Eq. ,
, ( is the initial condition of the tagged particle), and . The MSD for the tagged particle is obtained directly from Eq. :Heterogeneous files
The solution of Eqs. - is approximated with the expression,Starting from Eq. , the PDF of the tagged particle in the heterogeneous file follows,
The MSD of a tagged particle in a heterogeneous file is taken from Eq. :
Renewal anomalous heterogeneous files
The results of renewal-anomalous files are simply derived from the results of Brownian files. Firstly, the PDF in Eq. is written in terms of the PDF that solves the un-convoluted equation, that is, the Brownian file equation; this relation is made in Laplace space:(The subscript nrml stands for normal dynamics.) From Eq. , it is straightforward relating the MSD of Brownian heterogeneous files and renewal-anomalous heterogeneous files,
From Eq. , one finds that the MSD of a file with normal dynamics in the power of
is the MSD of the corresponding renewal-anomalous file,Anomalous files with independent particles
The equation of motion for anomalous files with independent particles, , is very complicated. Solutions for such files are reached while deriving scaling laws and with numerical simulations.Scaling laws for anomalous files of independent particles
Firstly, we write down the scaling law for the mean absolute displacement (MAD) in a renewal file with a constant density :Here,
is the number of particles in the covered-length , and is the MAD of a free anomalous particle, . In Eq. , enters the calculations since all the particles within the distance from the tagged one must move in the same direction in order that the tagged particle will reach a distance from its initial position. Based on Eq. , we write a generalized scaling law for anomalous files of independent particles:The first term on the right hand side of Eq. appears also in renewal files; yet, the term f(n) is unique. f(n) is the probability that accounts for the fact that for moving n anomalous independent particles in the same direction, when these particles indeed try jumping in the same direction (expressed with the term, (
), the particles in the periphery must move first so that the particles in the middle of the file will have the free space for moving, demanding faster jumping times for those in the periphery. f(n) appears since there is not a typical timescale for a jump in anomalous files, and the particles are independent, and so a particular particle can stand still for a very long time, substantially limiting the options of progress for the particles around him, during this time. Clearly,, where f(n) = 1 for renewal files since the particles jump together, yet also in files of independent particles with , since in such files there is a typical timescale for a jump, considered the time for a synchronized jump. We calculate f(n) from the number of configurations in which the order of the particles’ jumping times enables motion; that is, an order where the faster particles are always located towards the periphery. For n particles, there are n! different configurations, where one configuration is the optimal one; so, . Yet, although not optimal, propagation is also possible in many other configurations; when m is the number of particles that move, then,where
counts the number of configurations in which those m particles around the tagged one have the optimal jumping order. Now, even when m~n/2, . Using in Eq. , ( a small number larger than 1), we see,(In Eq. , we use,
.) Equation shows that asymptotically the particles are extremely slow in anomalous files of independent particles.Numerical studies of anomalous files of independent particles
With numerical studies, one sees that anomalous files of independent particles form clusters. This phenomenon defines a dynamical phase transition. At steady state, the percentage of particles in cluster,, follows,In Figure 1 we show trajectories from 9 particles in a file of 501 particles. (It is recommended opening the file in a new window). The upper panels show trajetcories for
and the lower panels show trajectories for . For each value of shown are trajectories in the early stages of the simulations (left) and in all stages of the simulation (right). The panels exhibit the phenomenon of the clustering, where the trajectories attract each other and then move pretty much together.