View source for Omnitruncated simplectic honeycomb

The term '''file dynamics''' is the motion of many particles in a narrow channel.

In science: in [[chemistry]], [[physics]], [[mathematics]] and related fields, '''file dynamics''' (sometimes called, '''single file dynamics''') is the diffusion of ''N'' (''N'' → ∞) identical [[Brownian motion|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 famous statistical properties of this process is that the [[mean squared displacement]] (MSD) of a particle in the file follows,  <math> \mathrm{MSD}\approx t^{\frac{1}{2}} </math>, and its [[probability density function]] (''PDF'') is [[Gaussian function|Gaussian]] in position with a variance MSD.<ref>Harris T. E. (1965) "Diffusion with 'Collisions' between Particles", ''Journal of Applied Probability'', 2 (2), 323-338 {{jstor|3212197}}</ref><ref>Jepsen D. W., J. Math. Phys. (N.Y.), 6 (1965) 405.</ref><ref>Lebowitz J. L., Percus J. K., [http://prola.aps.org/abstract/PR/v155/i1/p122_1 Phys. Rev., 155 (1967) 122].</ref>

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, <math>MSD\approx t^{\frac{1+a}{2}}</math>, with a Gaussian ''PDF''.<ref name=OF_2008EPL />
* When, in addition, the particles' diffusion coefficients are distributed like a power law with exponent γ (around the origin), the ''MSD'' follows, <math>MSD\approx t^{\frac{1-\gamma}{2/(1+a)-\gamma}}</math>, with a Gaussian ''PDF''.<ref name=OF_2010PRE />
* 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&nbsp;−&nbsp;''α'', the MSD scales like the MSD of the corresponding normal file, in the power of α.<ref name=OF_2010PLA />
* In anomalous files of independent particles, the ''MSD'' is very slow and scales like, <math>MSD\approx log^2(t)</math>. 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, <math> \xi \approx \sqrt{1-\alpha^3} </math>.<ref name=OF_2011EPL />
* Other generalizations include: when the particles can bypass each other with a constant probability upon encounter, an enhanced diffusion is seen.<ref name=mon/> When the particles interact with the channel, a slower diffusion is observed.<ref name=taloni/> Files in embedded in two-dimensions show similar characteristics of files in one dimension.<ref name=OF_2011EPL/>

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:<ref>Kärger J. and Ruthven D. M. (1992) ''Diffusion in Zeolites and Other Microscopic Solids'' (Wiley, NY).</ref><ref>Wei Q. H., Bechinger C. and Leiderer P., [http://www.sciencemag.org/content/287/5453/625.short ''Science'', 287 (2000) 625].</ref><ref>de Gennes J. P., J. Chem. Phys., 55 (1971) 572.</ref><ref>Richards P. M., [http://prb.aps.org/abstract/PRB/v16/i4/p1393_1 Phys. Rev. B, 16 (1977) 1393].</ref><ref>Maxfield F. R., Curr. Opin. Cell Biol., 14 (2002) 483.</ref><ref>Biological Membrane Ion Channels: Dynamics, Structure, And Applications, Chung S-h., Anderson O. S. and Krishnamurthy V. V., editors (Springer-verlag) 2006.</ref><ref>Howard J., Mechanics of Motor Proteins and the Cytoskeleton (Sinauer associates Inc. Sunderland, MA) 2001.</ref> 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.

== Mathematical formulation ==

===Simple files===
In simple Brownian files, <math> P(\mathbf{x},t\mid \mathbf{x_0}) </math>, the joint [[probability density function]] (PDF) for all the particles in file, obeys a normal diffusion equation:

{{NumBlk|:|<math>
\partial_t P(\mathbf{x},t\mid \mathbf{x_0})=D \Sigma_{j=-M}^M \partial^2_{x_j} P(\mathbf{x},t\mid \mathbf{x_0}).</math>|{{EquationRef|1}}}}
                     	  	       
In <math> P(\mathbf{x},t\mid \mathbf{x_0}) </math>, <math> \mathbf{x} = \{x_{-M}, x_{-M+1},\ldots,x_M\} </math> is the set of particles' positions at time <math>t</math> and <math>\mathbf{x_0}</math> is the set of the particles' initial positions at the initial time <math>t_0</math> (set to zero). Equation (1) is solved with the appropriate boundary conditions, which reflect the hard-sphere nature of the file:

{{NumBlk|:|<math>
\big( D \partial_{x_j} P(\mathbf{x},t\mid \mathbf{x_0}) \big)_{x_j=x_{j+1}} = \big( D\partial_{x_{j+1}} P(\mathbf{x},t\mid \mathbf{x_0})\big)_{x_{j+1}=x_j}; \qquad   j=-M,\ldots,M-1, 
</math>|{{EquationRef|2}}}}

and with the appropriate initial condition:

{{NumBlk|:|<math>
P(\mathbf{x},t\rightarrow\infty\mid x_0) = \Pi_{j=-M}^M \delta(x_j-x_{0,j}).     
</math>|{{EquationRef|3}}}}

In a simple file, the initial density is fixed, namely,<math> x_{0,j}=j\Delta</math>, where <math>\Delta</math> is a parameter that represents a microscopic length. The PDFs' coordinates must obey the order: <math>x_{-M}\le x_{-M+1}\le\cdots\le x_M</math>.

===Heterogeneous files===
In such files, the equation of motion follows,

{{NumBlk|:|<math>
\partial_t P(\mathbf{x},t\mid \mathbf{x_0})= \Sigma_{j=-M}^M D_j\partial^2_{x_j} P(\mathbf{x},t\mid \mathbf{x_0}).</math>|{{EquationRef|4}}}}	
	       
with the boundary conditions:

{{NumBlk|:|<math>
\big( D_j \partial_{x_j} P(\mathbf{x},t\mid \mathbf{x_0}) \big)_{x_j=x_{j+1}} = \big( D_{j+1}\partial_{x_{j+1}} P(\mathbf{x},t\mid \mathbf{x_0})\big)_{x_{j+1}=x_j};\qquad   j=-M,\ldots,M-1, 
</math>|{{EquationRef|5}}}}

and with the initial condition, Eq. ({{EquationNote|3}}), where the particles’ initial positions obey:

{{NumBlk|:|<math>
x0,j=sign(j)\Delta\mid j\mid ^{1/(1-a)};\qquad 0\le \mathit{a}\le1.		
</math>|{{EquationRef|6}}}}

The file diffusion coefficients are taken independently from the PDF,

{{NumBlk|:|<math>
W(D)= \frac{1-\gamma}{\Lambda} (D/\Lambda)^{-\gamma},\qquad    0\le\gamma\le1,		
</math>|{{EquationRef|7}}}}

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 waiting time probability density function (WT-PDF; see [[Continuous-time Markov process]] for more information) of the form: <math>\psi_{\alpha}(t)\sim k (kt)^{-1-\alpha}, 0<\alpha\le 1</math>, 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 <math>k_{\alpha}(t)</math>:

{{NumBlk|:|<math>
\partial_t P(\mathbf{x},t\mid \mathbf{x_0})=\Sigma_{j=-M}^M D_j \partial_{x_j}^2 \int_0^t k_{\alpha}(t-u) P(\mathbf{x},u\mid \mathbf{x_0}) \, du.                 
</math>|{{EquationRef|8}}}}

Here, the kernel <math>k_{\alpha}(t)</math>and the WT-PDF <math>\psi_{\alpha}(t)</math> are related in Laplace space, <math>\bar{k}_{\alpha} (s)=\frac{s\bar{\psi}_{\alpha}(s)}{1-\bar{\psi}_{\alpha}(s)}</math>. (The Laplace transform of a function <math>f(t)</math> reads, <math>\bar{f}(s)=\int_0^\infty f(t) e^{-st} \, dt</math>.) The reflecting boundary conditions accompanied Eq. ({{EquationNote|8}}) are obtained when convoluting the boundary conditions of a Brownian file with the kernel <math>k_{\alpha}(t)</math>, 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 <math>\psi_{\alpha}(t)</math> (<math>\psi_{\alpha}(t)</math> 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, <math>t_i</math> for particle ''i'', attempts a jump. Then, the waiting times for all the other particles are adjusted: we subtract <math>t_i</math> 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:

{{NumBlk|:|<math>
\partial_{t_i} P(\mathbf{x},\mathbf{t}\mid \mathbf{x_0}) = D_i \partial_{x_i}^2 \int_0^{t_i} k_{\alpha}(t_i-u_i) P(\mathbf{x},\mathbf{t}^{'(i)},u_i\mid \mathbf{x_0}) \, du_i;\qquad -M\le i \le M.           
</math>|{{EquationRef|9}}}}

Note that the time argument in the PDF <math>P(\mathbf{x},\mathbf{t}\mid \mathbf{x_0})</math> is a vector of times: <math>\mathbf{t}=\{t_i\}_{i=-M}^M</math>, and <math> \mathbf{t}^{'(i)}=\{t_c \}_{c=-M, c\ne i}^M</math>. 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. ({{EquationNote|9}}) is very complicated, and averaging further complicates things.

== Mathematical analysis ==

===Simple files===
The solution of Eqs. ({{EquationNote|1}})-({{EquationNote|2}}) is a complete set of permutations of all initial coordinates appearing in the Gaussians,<ref name=OF_2008EPL/>

{{NumBlk|:|<math>
P(\mathbf{x},\mid\mathbf{x_0}) = \frac{1}{c_N}  \Sigma_{\{p\}} e^{-1/4Dt \Sigma_{j=-M}^M (x_j - x_{0,j}(p))^2}. 	
</math>|{{EquationRef|10}}}}

Here, the index <math>p</math> goes on all the permutations of the initial coordinates, and contains <math>N!</math> permutations. From Eq.  ({{EquationNote|10}}), the PDF of a tagged particle in the file, <math>P(r,t\mid r_0)</math>, is calculated <ref name= OF_2008EPL/>

{{NumBlk|:|<math>
P(r,t\mid r_0 ) \sim \frac{1}{c_N}  e^{-R_d^2/\sqrt{2\tau}},
</math>|{{EquationRef|11}}}}

In Eq. ({{EquationNote|11}}), <math> R_d = r_d\Delta </math>, <math> r_d = r - r_0</math> (<math>r_0</math> is the initial condition of the tagged particle), and <math>\tau = \Delta^{-2} Dt</math>. The MSD for the tagged particle is obtained directly from Eq. ({{EquationNote|11}}):

{{NumBlk|:|<math>
\langle R_d^2\rangle\sim \sqrt{\tau}.
</math>|{{EquationRef|12}}}}

===Heterogeneous files===
The solution of Eqs. ({{EquationNote|4}})-({{EquationNote|7}}) is approximated with the expression,<ref name=OF_2010PRE/>

{{NumBlk|:|<math>
P(x,t\mid x_0) \approx \frac{1}{c_N}  \Sigma_{\{p\}} e^{-\Sigma_{j=-M}^M (x_j-x_{0,j}(p))^2/4tD_j}.
</math>|{{EquationRef|13}}}}

Starting from Eq. ({{EquationNote|13}}), the PDF of the tagged particle in the heterogeneous file follows,<ref name=OF_2010PRE/>

{{NumBlk|:|<math>
P(r,t\mid r_0 ) \sim \frac{1}{c_N}  e^{-R_d^2/4\tau \tau^{(1-a))/(2-\gamma(1+a))} }; \qquad \tau=\Delta^{-2}\Lambda t.
</math>|{{EquationRef|14}}}}

The MSD of a tagged particle in a heterogeneous file is taken from Eq. ({{EquationNote|14}}):

{{NumBlk|:|<math>
\langle R_d^2\rangle =2\tau^{(1-\gamma)/(2c-\gamma)},  \qquad  c=1/((1+a)).	     
</math>|{{EquationRef|15}}}}

===Renewal anomalous heterogeneous files===
The results of renewal-anomalous files are simply derived from the results of Brownian files. Firstly, the ''PDF'' in Eq. ({{EquationNote|8}}) 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:

{{NumBlk|:|<math>
\bar{P}(x,s\mid x_0) = \frac{1}{\bar{k}_{\alpha}(s)} \bar{P}_\text{nrml}(x,s/\bar{k}_\alpha(s)\mid x_0).
</math>|{{EquationRef|16}}}}

(The subscript ''nrml'' stands for normal dynamics.) From Eq. ({{EquationNote|16}}), it is straightforward relating the ''MSD'' of Brownian heterogeneous files and renewal-anomalous heterogeneous files,<ref name=OF_2010PLA/>

{{NumBlk|:|<math>
\langle\bar{r}^2(s)\rangle = \frac{1}{\bar{k}_\alpha(s)}\langle\bar{r}^2(s/\bar{k}_\alpha(s))\rangle_\text{nrml}.
</math>|{{EquationRef|17}}}}

From Eq. ({{EquationNote|18}}), one finds that the ''MSD'' of a file with normal dynamics in the power of <math>\alpha</math> is the ''MSD'' of the corresponding renewal-anomalous file,<ref name=OF_2010PLA/>

{{NumBlk|:|<math>
\langle r^2 (t)\rangle \sim \langle r^2 (t) \rangle_\text{nrml}^\alpha. 						
</math>|{{EquationRef|19}}}}

===Anomalous files with independent particles===
The equation of motion for anomalous files with independent particles, ({{EquationNote|9}}), 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:<ref name=OF_2008EPL/><ref name=OF_2010PRE/><ref name=OF_2011EPL/>

{{NumBlk|:|<math>
\langle\mid r \mid\rangle \sim \langle\mid r \mid\rangle_\text{free}/n.		
</math>|{{EquationRef|20}}}}

Here, <math>n</math> is the number of particles in the covered-length <math>\langle\mid r\mid \rangle </math>, and <math> \langle\mid r\mid \rangle_\text{free}</math> is the ''MAD'' of a free anomalous particle, <math> \langle\mid r\mid \rangle_\text{free}\sim t^{\alpha/2})</math>. In Eq. ({{EquationNote|20}}), <math>n</math> enters the calculations since all the particles within the distance <math>\langle\mid r\mid\rangle</math> from the tagged one must move in the same direction in order that the tagged particle will reach a distance <math>\langle\mid r\mid\rangle</math> from its initial position. Based on Eq. ({{EquationNote|20}}), we write a generalized scaling law for anomalous files of independent particles:

{{NumBlk|:|<math>
\langle |r|\rangle \sim \frac{\langle\mid r\mid\rangle_\text{free}}{n} f(n);\qquad 0<f(n)<1.
</math>|{{EquationRef|21}}}}

The first term on the right hand side of Eq. ({{EquationNote|21}}) 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, (<math>\langle\mid r\mid \rangle_\text{free}/n)</math>), 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,<math> 0<f(n)<1</math>, where ''f''(''n'')&nbsp;=&nbsp;1 for renewal files since the particles jump together, yet also in files of independent particles with  <math>\alpha>1</math>, 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, <math>(1/n!)\le f(n)</math>. Yet, although not optimal, propagation is also possible in many other configurations; when m is the number of particles that move, then,

{{NumBlk|:|<math>
f(n)\sim\dbinom{n}{m} (n-m)! 1/n!,
</math>|{{EquationRef|22}}}}

where <math>\dbinom{n}{m}(n-m)!</math> 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, <math> f(n)\sim e^{-n/2}</math>. Using in Eq. ({{EquationNote|21}}), <math>f(n)\sim e^{-n/n_0}</math> (<math>n_0</math> a small number larger than 1), we see,

{{NumBlk|:|<math>
\mathrm{MSD} \sim \left(\frac{\alpha}{n_0}\right)^2 \ln^2(t).				
</math>|{{EquationRef|23}}}}

(In Eq. ({{EquationNote|23}}), we use, <math>MSD\sim \mid MAD\mid^2</math>.) Equation ({{EquationNote|23}})  shows that asymptotically the particles are extremely slow in anomalous files of independent particles.

====Numerical studies of anomalous files of independent particles====
[[File:File indep par.pdf|thumb|right|'''Figure 1''' Trajectories from simulations of 501 anomalous, independent, particles, with <math>\alpha=0.9, 0.1</math> (recommended: open the file in a new window)]]

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, <math>\xi(\alpha)</math>, follows,

{{NumBlk|:|<math>
\xi(\alpha) \approx \sqrt{1 - \alpha^3}						
</math>|{{EquationRef|24}}}}

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 <math>\alpha=0.9</math> and the lower panels show trajectories for <math>\alpha=0.1</math>. For each value of <math>\alpha</math> 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.

== See also ==
*[[Brownian motion]]
*[[Langevin dynamics]]
*[[System dynamics]]

== References ==
{{reflist|refs=

<ref name=mon>Mon K. K. and Percus J. K., [http://jcp.aip.org/resource/1/jcpsa6/v117/i5/p2289_s1?isAuthorized=no J. Chem. Phys., 117 (2002) 2289].</ref> 
<ref name=OF_2008EPL >Flomenbom O. and Taloni A., [http://iopscience.iop.org/0295-5075/83/2/20004 Europhys. Lett., 83 (2008) 20004].</ref>
<ref name=OF_2010PRE>Flomenbom O., [http://pre.aps.org/abstract/PRE/v82/i3/e031126 Phys. Rev. E, 82 (2010) 031126].</ref>
<ref name=taloni>Taloni A. and Marchesoni F., [http://prl.aps.org/abstract/PRL/v96/i2/e020601 Phys. Rev. Lett., 96 (2006) 020601].</ref>
<ref name=OF_2010PLA >Flomenbom O., [http://www.sciencedirect.com/science?_ob=ArticleURL&_udi=B6TVM-50T41FC-1&_user=10&_coverDate=09%2F20%2F2010&_rdoc=12&_fmt=high&_orig=browse&_origin=browse&_zone=rslt_list_item&_srch=doc-info(%23toc%235538%232010%23996259957%232378741%23FLA%23display%23Volume)&_cdi=5538&_sort=d&_docanchor=&_ct=20&_acct=C000050221&_version=1&_urlVersion=0&_userid=10&md5=3cb1301417736a056fe0e0da19d86269&searchtype=a Phys. Lett. A, 374 (2010) 4331].</ref>
<ref name=OF_2011EPL >Flomenbom O., [http://iopscience.iop.org/0295-5075/94/5/58001/?rss=2.0 EPL 94, 58001 (2011)].</ref>

}}

[[Category:Diffusion]]
[[Category:Statistical mechanics]]
[[Category:Stochastic processes]]