Anomalous reaction-diffusion as a model of nonexponential DNA escape kinetics

We show that data from recent experiments carried out on the kinetics of DNA escape from alpha-hemolysin nanopores [M. Wiggin, C. Tropini, C. T. Cossa, N. N. Jetha, and A. Marziali, Biophys. J. 95, 5317 (2008)] may be rationalized by a model of chain dynamics based on the anomalous diffusion of a particle moving in a harmonic well in the presence of a delta function sink. The experiments of Wiggin found, among other things, that the occasional occurrence of unusually long escape times in the distribution of chain trapping events led to nonexponential decays in the survival probability, S(t), of the DNA molecules within the nanopore. Wiggin ascribed this nonexponentiality to the existence of a distribution of trapping potentials, which they suggested was theresult of stochastic interactions between the bases of the DNA and the amino acids located on the surface of the nanopore. Based on this idea, they showed that the experimentally determined S(t) could be well fit in both the short and long time regimes by a function of the form (1+t/tau)(-alpha) (the so called Becquerel function). In our model, S(t) is found to be given by a Mittag-Leffler function at short times and by a generalized Mittag-Leffler function at long times. By suitable choice of certain parameter values, these functions are found to fit the experimental S(t) even better than the Becquerel function. Anomalous diffusion of DNA within the trap prior to escape over a barrier of fixed height may therefore provide a second, plausible explanation of the data, and may offer fresh perspectives on similar trapping and escape problems.

Dynamics of zeolite cage and its effect on the diffusion properties of sorbate: persistence of diffusion anomaly in NaA zeolite

Recent computer simulations on zeolites Y and A have found that the diffusion coefficient and the rate of intercage diffusion exhibit, apart from a linear dependence on the reciprocal of the square of the sorbate diameter, an anomalous peak as sorbate diameter approaches the window diameter. Here we report molecular dynamics simulations of zeolite NaA incorporating framework flexibility as a function of sorbate diameter in order to verify the existence of anomalous diffusion. Results suggest persistence of anomalous diffusion or ring effect. This suggests that the anomalous behavior is a general effect characteristic of zeolites Y and A. The barrier for diffusion across the eight-ring window is seen to be negative and is found to decrease with sorbate size. The effect of sorbate on the cage motion has also been investigated. Results suggest that the window expands during intercage migration only if the sorbate size is comparable to the window diameter. Flexible cage simulations yield a higher value for the diffusion coefficient and also the rate of intercage diffusion. This increase has been shown to be due to an increase in the intercage diffusions via the centralized diffusion mode rather than the surface-mediated mode. It is shown that this increase arises from an increase in the single particle density distribution in the region near the cage center.

Sintering of beta titanium: Role of dislocation pipe diffusion

Sintering of titanium in its high temperature beta phase was studied by isothermal dilatometry. The sintering shrinkage y did not follow the normal time exponent type of behaviour, instead being described by the equation y = Kt(m)/[1-(A+Bt)(2)], where m = 1.93 +/- 0.07, with an activation energy of 62-90 kJ mol(-1). A detailed analysis of these results, based on the 'anomalous' diffusion behaviour reported for beta titanium, is carried out. It is shown that the generation of a high density of dislocations during the alpha --> beta phase transformation, coupled with sluggish recovery at the sintering necks, enables sintering mass transport by pipe diffusion through dislocation cores from sources of matter within the particles to become dominant.

First passage distributions for long memory processes

We study the distribution of first passage time for Levy type anomalous diffusion. A fractional Fokker-Planck equation framework is introduced.For the zero drift case, using fractional calculus an explicit analytic solution for the first passage time density function in terms of Fox or H-functions is given. The asymptotic behaviour of the density function is discussed. For the nonzero drift case, we obtain an expression for the Laplace transform of the first passage time density function, from which the mean first passage time and variance are derived.

Handling large variations in mechanics: Some applications

There is a need to use probability distributions with power-law decaying tails to describe the large variations exhibited by some of the physical phenomena. The Weierstrass Random Walk (WRW) shows promise for modeling such phenomena. The theory of anomalous diffusion is now well established. It has found number of applications in Physics, Chemistry and Biology. However, its applications are limited in structural mechanics in general, and structural engineering in particular. The aim of this paper is to present some mathematical preliminaries related to WRW that would help in possible applications. In the limiting case, it represents a diffusion process whose evolution is governed by a fractional partial differential equation. Three applications of superdiffusion processes in mechanics, illustrating their effectiveness in handling large variations, are presented.

A Model of Anomalous Chain Translocation Dynamics

A model of polymer translocation based on the stochastic dynamics of the number of monomers on one side of a pore-containing surface is formulated in terms of a one-dimensional generalized Langevin equation, in which the random force is assumed to be characterized by long-ranged temporal correlations. The model is introduced to rationalize anomalies in measured and simulated values of the average time of passage through the pore, which in general cannot be satisfactorily accounted for by simple Brownian diffusion mechanisms. Calculations are presented of the mean first passage time for barrier crossing and of the mean square displacement of a monomeric segment, in the limits of strong and weak diffusive bias. The calculations produce estimates of the exponents in various scaling relations that are in satisfactory agreement with available data.

A Model of Anomalous Enzyme-Catalyzed Gel Degradation Kinetics

We show that a model of target location involving n noninteracting particles moving subdiffusively along a line segment (a generalization of a model introduced by Sokolov et al. [Biophys. J. 2005, 89, 895.]) provides a basis for understanding recent experiments by Pelta et al. [Phys. Rev. Lett. 2007, 98, 228302.] on the kinetics of diffusion-limited gel degradation. These experiments find that the time t(c) taken by the enzyme thermolysin to completely hydrolyze a gel varies inversely as roughly the 3/2 power of the initial enzyme concentration [E]. In general, however, this time would be expected to vary either as [E](-1) or as [E](-2), depending on whether the Brownian diffusion of the enzyme to the site of cleavage took place along the network chains (1-d diffusion) or through the pore spaces (3-d diffusion). In our model, the unusual dependence of t(c) on [E] is explained in terms of a reaction-diffusion equation that is formulated in terms of fractional rather than ordinary time derivatives.

Anomalous behavior of hydrodynamic modes in the two dimensional shear flow of a granular material

The growth rates of the hydrodynamic modes in the homogeneous sheared state of a granular material are determined by solving the Boltzmann equation. The steady velocity distribution is considered to be the product of the Maxwell Boltzmann distribution and a Hermite polynomial expansion in the velocity components; this form is inserted into them Boltzmann equation and solved to obtain the coeificients of the terms in the expansion. The solution is obtained using an expansion in the parameter epsilon =(1 - e)(1/2), and terms correct to epsilon(4) are retained to obtain an approximate solution; the error due to the neglect of higher terms is estimated at about 5% for e = 0.7. A small perturbation is placed on the distribution function in the form of a Hermite polynomial expansion for the velocity variations and a Fourier expansion in the spatial coordinates: this is inserted into the Boltzmann equation and the growth rate of the Fourier modes is determined. It is found that in the hydrodynamic limit, the growth rates of the hydrodynamic modes in the flow direction have unusual characteristics. The growth rate of the momentum diffusion mode is positive, indicating that density variations are unstable in the limit k--> 0, and the growth rate increases proportional to kslash} k kslash}(2/3) in the limit k --> 0 (in contrast to the k(2) increase in elastic systems), where k is the wave vector in the flow direction. The real and imaginary parts of the growth rate corresponding to the propagating also increase proportional to kslash k kslash(2/3) (in contrast to the k(2) and k increase in elastic systems). The energy mode is damped due to inelastic collisions between particles. The scaling of the growth rates of the hydrodynamic modes with the wave vector I in the gradient direction is similar to that in elastic systems. (C) 2000 Elsevier Science B.V. All rights reserved.

Large-amplitude chirped coherent phonons in tellurium mediated by ultrafast photoexcited carrier diffusion

We report femtosecond time-resolved reflectivity measurements of coherent phonons in tellurium performed over a wide range of temperatures (3-296 K) and pump-laser intensities. A totally symmetric A(1) coherent phonon at 3.6 THz responsible for the oscillations in the reflectivity data is observed to be strongly positively chirped (i.e., phonon time period decreases at longer pump-probe delay times) with increasing photoexcited carrier density, more so at lower temperatures. We show that the temperature dependence of the coherent phonon frequency is anomalous (i.e, increasing with increasing temperature) at high photoexcited carrier density due to electron-phonon interaction. At the highest photoexcited carrier density of (1.4 x 10(21) cm(-3) and the sample temperature of 3 K, the lattice displacement of the coherent phonon mode is estimated to be as high as similar to 0.24 angstrom. Numerical simulations based on coupled effects of optical absorption and carrier diffusion reveal that the diffusion of carriers dominates the nonoscillatory electronic part of the time-resolved reflectivity. Finally, using the pump-probe experiments at low carrier density of 6 x 10(18) cm(-3), we separate the phonon anharmonicity to obtain the electron-phonon coupling contribution to the phonon frequency and linewidth.

Diffusion of sorbates in zeolites Y and A novel dependence on sorbate size and strength of sor ba t e-Zeoli t e interaction

A systematic investigation of monatomic spherical sorbates in the supercages of zeolites Y and A by molecular dynamics technique is presented. Rates of intercage diffusion, rates of cage visits, and the diffusion coefficients have been calculated as a function of the sorbate-zeolite interaction strength. These properties exhibit markedly different dependences on interaction strength for the two zeolites. The observed behavior is shown to be a consequence of the two principal mechanisms of intercage diffusion and the energetic barrier associated with them. The diffusion coefficient and other properties associated with intercage diffusion are found to be directly proportional to the reciprocal of the square of the sorbate diameter when the sorbate diameter is significantly smaller than the window diameter. As the sorbate diameter increases, a peak is observed in all the transport properties investigated including the diffusion coefficient. We call this surprising effect as the ring or levitation effect and it explains several anomalous results reported in the literature and suggests a breakdown of the geometrical criterion for diffusion of sorbates. It shows that under certain conditions nongeometrical factors play a major role and geometrical factors become secondary in the determination of the molecular sieve property. A generalized parameter has been proposed which suggests conditions under which one can expect the ring or levitation effect in any porous medium. Inverse size selectivity becomes operative under these conditions.

Diffusion Anomaly in Silicalite and VPI-5 from Molecular Dynamics Simulations

Molecular dynamics (MD) simulations on rigid and flexible framework models of silicalite and a rigid framework model of the aluminophosphate VPI-5 for different sorbate diameters are reported. The sorbate-host interactions are modeled in terms of simple atom-atom Lennard-Jones interactions. The results suggest that the diffusion coefficient exhibits an anomaly as gamma approaches unity. The MD results confirm the existence of a linear regime for sorbate diameters significantly smaller than the channel diameter and an anomalous regime observed for sorbate diameters comparable to the channel diameter. The power spectra obtained by Fourier transformation of the velocity autocorrelation function indicate that there is an increase in the intensity of the low-frequency component for the velocity component parallel to the direction of motion for the sorbate diameter in the anomalous regime. The present results suggest that the diffusion anomaly is observed irrespective of (1) the geometry and topology of the pore structure and (2) the nature of the host material. The results are compared with the work of Derouane and co-workers, who have suggested the existence of ''floating molecules'' on the basis of earlier theoretical and computational approaches.

Memory-induced anomalous dynamics in a minimal random walk model

A discrete-time dynamics of a non-Markovian random walker is analyzed using a minimal model where memory of the past drives the present dynamics. In recent work N. Kumar et al., Phys. Rev. E 82, 021101 (2010)] we proposed a model that exhibits asymptotic superdiffusion, normal diffusion, and subdiffusion with the sweep of a single parameter. Here we propose an even simpler model, with minimal options for the walker: either move forward or stay at rest. We show that this model can also give rise to diffusive, subdiffusive, and superdiffusive dynamics at long times as a single parameter is varied. We show that in order to have subdiffusive dynamics, the memory of the rest states must be perfectly correlated with the present dynamics. We show explicitly that if this condition is not satisfied in a unidirectional walk, the dynamics is only either diffusive or superdiffusive (but not subdiffusive) at long times.

Anomalous transport at weak coupling

We evaluate the contribution of chiral fermions in d = 2, 4, 6, chiral bosons, a chiral gravitino like theory in d = 2 and chiral gravitinos in d = 6 to all the leading parity odd transport coefficients at one loop. This is done by using finite temperature field theory to evaluate the relevant Kubo formulae. For chiral fermions and chiral bosons the relation between the parity odd transport coefficient and the microscopic anomalies including gravitational anomalies agree with that found by using the general methods of hydrodynamics and the argument involving the consistency of the Euclidean vacuum. For the gravitino like theory in d = 2 and chiral gravitinos in d = 6, we show that relation between the pure gravitational anomaly and parity odd transport breaks down. From the perturbative calculation we clearly identify the terms that contribute to the anomaly polynomial, but not to the transport coefficient for gravitinos. We also develop a simple method for evaluating the angular integrals in the one loop diagrams involved in the Kubo formulae. Finally we show that charge diffusion mode of an ideal 2 dimensional Weyl gas in the presence of a finite chemical potential acquires a speed, which is equal to half the speed of light.

Anomalous lifetime statistics of ligand-receptor binding in a tethered polymer system

In this paper, motivated by observations of non-exponential decay times in the stochastic binding and release of ligand-receptor systems, exemplified by the work of Rogers et al on optically trapped DNA-coated colloids (Rogers et al 2013 Soft Matter 9 6412), we explore the general problem of polymer-mediated surface adhesion using a simplified model of the phenomenon in which a single polymer molecule, fixed at one end, binds through a ligand at its opposite end to a flat surface a fixed distance L away and uniformly covered with receptor sites. Working within the Wilemski-Fixman approximation to diffusion-controlled reactions, we show that for a flexible Gaussian chain, the predicted distribution of times f(t) for which the ligand and receptor are bound is given, for times much shorter than the longest relaxation time of the polymer, by a power law of the form t(-1/4). We also show when the effects of chain stiffness are incorporated into this model (approximately), the structure of f(t) is altered to t(-1/2). These results broadly mirror the experimental trends in the work cited above.

Determination of diffusion parameters and activation energy of diffusion in V3Si phase with A15 crystal structure

Diffusion such is the integrated diffusion coefficient of the phase, the tracer diffusion coefficient of species at different temperatures and the activation energy for diffusion, are determined in V3Si phase with A15 crystal structure. The tracer diffusion coefficient of Si Was found to be negligible compared to the tracer diffusion coefficient of V. The calculated diffusion parameters will help to validate the theoretical analysis of defect structure of the phase, which plays an important role in the superconductivity.