Observation of anomalous diffusion and fractional self-similarity in one dimension

We experimentally study anomalous diffusion of ultracold atoms in a one dimensional polarization optical lattice. The atomic spatial distribution is recorded at different times and its dynamics and shape are analyzed. We find that the width of the cloud exhibits a power-law time dependence with an exponent that depends on the lattice depth. Moreover, the distribution exhibits fractional self-similarity with the same characteristic exponent. The self-similar shape of the distribution is found to be well fitted by a Lévy distribution, but with a characteristic exponent that differs from the temporal one. Numerical simulations suggest that this is due to long trapping times in the lattice and correlations between the atom's velocity and flight duration.     

Strong and Weak Chaos in Weakly Nonintegrable Many-Body Hamiltonian Systems

We study properties of chaos in generic one-dimensional nonlinear Hamiltonian lattices comprised of weakly coupled nonlinear oscillators by numerical simulations of continuous-time systems and symplectic maps. For small coupling, the measure of chaos is found to be proportional to the coupling strength and lattice length, with the typical maximal Lyapunov exponent being proportional to the square root of coupling. This strong chaos appears as a result of triplet resonances between nearby modes. In addition to strong chaos we observe a weakly chaotic component having much smaller Lyapunov exponent, the measure of which drops approximately as a square of the coupling strength down to smallest couplings we were able to reach. We argue that this weak chaos is linked to the regime of fast Arnold diffusion discussed by Chirikov and Vecheslavov. In disordered lattices of large size we find a subdiffusive spreading of initially localized wave packets over larger and larger number of modes. The relations between the exponent of this spreading and the exponent in the dependence of the fast Arnold diffusion on coupling strength are analyzed. We also trace parallels between the slow spreading of chaos and deterministic rheology.     

Conductivity exponent in three-dimensional percolation by diffusion based on molecular trajectory algorithm and blind-ant rules

Diffusion on random systems above and at their percolation threshold in three dimensions is carried out by a molecular trajectory method and a simple lattice random walk method, respectively. The classical regimes of diffusion on percolation near the threshold are observed in our simulations by both methods. Our Monte Carlo simulations by the simple lattice random walk method give the conductivity exponent μ/ν=2.32±0.02 for diffusion on the incipient infinite clusters and μ/ν=2.21±0.03 for diffusion on a percolating lattice above the threshold. However, while diffusion is performed by the molecular trajectory algorithm either on the incipient infinite clusters or on a percolating lattice above the threshold, the result is found to be μ/ν=2.26±0.02. In addition, it takes less time step for diffusion based on the molecular trajectory algorithm to reach the asymptotic limit comparing with the simple lattice random walk.     

Dynamical behavior of the multiplicative diffusion coupled map lattices

We report a dynamical study of multiplicative diffusion coupled map lattices with the coupling between the elements only through the bifurcation parameter of the mapping function. We discuss the diffusive process of the lattice from an initially random distribution state to a homogeneous one as well as the stable range of the diffusive homogeneous attractor. For various coupling strengths we find that there are several types of spatiotemporal structures. In addition, the evolution of the lattice into chaos is studied. A largest Lyapunov exponent and a spatial correlation function have been used to characterize the dynamical behavior. (c) 1996 American Institute of Physics.     

On the non-universality of a critical exponent for self-avoiding walks

We have extended the enumeration of self-avoiding walks on the Manhattan lattice from 28 to 53 steps and for self-avoiding polygons from 48 to 84 steps. Analysis of this data suggests that the walk generating function exponent γ = 1.3385 ± 0.003, which is different from the corresponding exponent on the square, triangular and honeycomb lattices. This provides numerical support for an argument recently advanced by Cardy, to the effect that excluding walks with parallel nearest-neighbour steps should cause a change in the exponent γ. The lattice topology of the Manhattan lattice precludes such parallel steps.     

Electronic shot noise in fractal conductors

By solving a master equation in the Sierpiński lattice and in a planar random-resistor network, we determine the scaling with size L of the shot noise power P due to elastic scattering in a fractal conductor. We find a power-law scaling P proportional, variantL;{d_{f}-2-alpha}, with an exponent depending on the fractal dimension d_{f} and the anomalous diffusion exponent alpha. This is the same scaling as the time-averaged current I[over ], which implies that the Fano factor F=P/2eI[over ] is scale-independent. We obtain a value of F=1/3 for anomalous diffusion that is the same as for normal diffusion, even if there is no smallest length scale below which the normal diffusion equation holds. The fact that F remains fixed at 1/3 as one crosses the percolation threshold in a random-resistor network may explain recent measurements of a doping-independent Fano factor in a graphene flake.     

Size Dependent Heat Conduction in One-Dimensional Diatomic Lattices

We study the size dependency of heat conduction in one-dimensional diatomic FPU-β lattices and establish that for low dimensional material,contribution from optical phonons is found more effective to the thermal conductivity and enhance heat transport in the thermodynamic limit N →∞.For the finite size,thermal conductivity of 1D diatomic lattice is found to be lower than 1D monoatomic chain of the same size made up of the constituent particle of the diatomic chain.For the present 1D diatomic chain,obtained value of power divergent exponent of thermal conductivity0.428±0.001 and diffusion exponent 1.2723 lead to the conclusions that increase in the system size,increases the thermal conductivity and existence of anomalous energy diffusion.Existing numerical data supports our findings.     

Dynamical chaos in the Lorentz lattice gas

This paper provides an introduction to the applications of dynamical systems theory to nonequilibrium statistical mechanics, in particular to a study of nonequilibrium phenomena in Lorentz lattice gases with stochastic collision rules. Using simple arguments, based upon discussions in the mathematical literature, we show that such lattice gases belong to the category of dynamical systems with positive Lyapunov exponents. This is accomplished by showing how such systems can be expressed in terms of continuous phase space variables. Expressions for the Lyapunov exponent of a one-dimensional Lorentz lattice gas with periodic boundaries are derived. Other quantities of interest for the theory of irreversible processes are discussed.     

Exact solutions of triple-order time-fractional differential equations for anomalous relaxation and diffusion I: The accelerating case

In recent years the interest around the study of anomalous relaxation and diffusion processes is increased due to their importance in several natural phenomena. Moreover, a further generalization has been developed by introducing time-fractional differentiation of distributed order which ranges between 0 and 1. We refer to accelerating processes when the driving power law has a changing-in-time exponent whose modulus tends from less than 1 to 1, and to decelerating processes when such an exponent modulus decreases in time moving away from the linear behaviour. Accelerating processes are modelled by a time-fractional derivative in the Riemann-Liouville sense, while decelerating processes by a time-fractional derivative in the Caputo sense. Here the focus is on the accelerating case while the decelerating one is considered in the companion paper. After a short reminder about the derivation of the fundamental solution for a general distribution of time-derivative orders, we consider in detail the triple-order case for both accelerating relaxation and accelerating diffusion processes and the exact results are derived in terms of an infinite series of H-functions. The method adopted is new and it makes use of certain properties of the generalized Mittag-Leffler function and the H-function, moreover it provides an elegant generalization of the method introduced by Langlands (2006) [T.A.M. Langlands, Physica A 367 (2006) 136] to study the double-order case of accelerating diffusion processes.     

Behavior of general one-dimensional diffusion processes

We develop simple rigorous techniques to estimate the behavior of general one-dimensional diffusion processes. Any one-dimensional diffusion process (with drift) can be mapped onto a symmetric diffusion through an explicit change of variable. For such processes we can estimate explicitly the diffusion exponent, the recurrence properties, and the large fluctuations. In a second part, we apply these results to different models (including the Sinaï random walk: diffusion in a random drift) and we show how the main features of the diffusion can be readily handled.     

Universality in diffusion front growth dynamics

We have studied the scaling properties of diffusion fronts by numerical calculations based on the mean field approach in the context of a lattice gas model, performed in a triangular lattice. We find that the height-height correlation function scales with time t and length l as C(l, t) ≈l ^α f (t/l ^α/β) with α = 0.62±0.01 and β = 0.39±0.02. These exponent values are identical to those characterising the roughness of the diffusion fronts evolving through a square lattice [1,2], thus confirming their universality. Received 14 November 2001 / Received in final form 20 April 2002 Published online 31 July 2002     

Mobility of large clusters on a semiconductor surface:Kinetic Monte Carlo simulation results

The mobility of clusters on a semiconductor surface for various values of cluster size is studied as a function of temperature by kinetic Monte Carlo method. The cluster resides on the surface of a square grid. Kinetic processes such as the diffusion of single particles on the surface, their attachment and detachment to/from clusters, diffusion of particles along cluster edges are considered. The clusters considered in this study consist of 150–6000 atoms per cluster on average.A statistical probability of motion to each direction is assigned to each particle where a particle with four nearest neighbors is assumed to be immobile. The mobility of a cluster is found from the root mean square displacement of the center of mass of the cluster as a function of time. It is found that the diffusion coefficient of clusters goes as D = A(T)Nαwhere N is the average number of particles in the cluster, A(T) is a temperature-dependent constant and α is a parameter with a value of about-0.64 α -0.75. The value of α is found to be independent of cluster sizes and temperature values(170–220 K)considered in this study. As the diffusion along the perimeter of the cluster becomes prohibitive, the exponent approaches a value of-0.5. The diffusion coefficient is found to change by one order of magnitude as a function of cluster size.     

Concentration and temperature dependence of diluent dynamics studied by line shape experiments on a phosphate ester in polycarbonate

³¹P Hahn spin echo line shape and proton line shape experiments are reported on bisphenol A polycarbonate (BPAPC)-tris(2-ethylhexyl)phosphate (TOP) systems to study the concentration and temperature dependence of the local dynamics. In an earlier ³¹P line shape study a lattice model was presented as a framework to interpret the plasticization and antiplasticization behavior of the diluent based on a fractional population given by the type of nearest neighbor contacts in the mixed polymer-diluent glass. In this study, ³¹P spin echo line shapes of BPAPC, with 5%, 10% and 15% TOP, which monitor the diluent dynamics, at different temperatures and echo delay times are simulated in terms of fast- and slow-moving components, and the resulting fractional populations are compared with that predicted by the lattice model. Comparisons with the lattice model calculations are also made in the simulation of the ¹H line shapes on BPAPC with 5% and 10% TOP, which probes both the polymer and diluent dynamics, and on BPAPC with 5% and 10% perdeuterated trioctylphosphate (DTOP), which detects only the polymer motion. Fairly good line shape simulations and agreement between the lattice model and the fitting results at low diluent concentrations are obtained in all cases. Restricted cone motion best describes the slow-moving component in the ³¹P line shape fittings. For the fast component, rotational Brownian diffusion with a distribution of correlation times corresponding to a stretched exponential function is used. An activation energy E_a of 56 kJ/mol and an exponent of 0.7 for the fractional exponential correlation function are obtained and used to calculate the mechanical loss peak which was compared with the experimental loss data. The plateau character of the fractional population as a function of temperature can also be interpreted and understood in terms of the lattice model.     

Correlation decay in quantum chaotic billiards with bulk or surface disorder

We study the decay properties of correlation functions in quantum billiards with surface or bulk disorder. The quantum system is modeled by means of a tight-binding Hamiltonian with diagonal disorder, solved on LxL clusters of the square lattice. The correlation function is calculated by launching the system at t=0 into a wave function of the regular (clean) system and following its time evolution. The results show that the correlation function decays exponentially with a characteristic correlation time (inverse of the Lyapunov exponent lambda). For small enough disorder the Lyapunov exponent is approximately given by the imaginary part of the self-energy induced by disorder. On the other hand, if the scaling of the Lyapunov exponent with L is investigated by keeping constant l/L, where l is the mean free path, the results show that lambda is proportional to 1/L.     

Numerical simulation of the modulation transfer function in planar InGaAs dense arrays

Three-dimensional simulation methodology has been used to evaluate the performance of lattice matched InGaAs/InP double layer planar heterointerface detector arrays. The device characteristics under optical illumination and dark conditions have been computed. The modulation transfer function (MTF) profiles have been calculated with varying device geometries and carrier dynamics. It is found that the p well diffusion radius and minority carrier recombination play important roles in the MTF behaviors of dense arrays. Moderate p well diffusion dimension should be used to balance the device performances between the dark current and MTF profile. Moreover, better MTF characteristic under low light condition can be achieved with higher quality material which has longer recombination lifetime. The influences of underlying mechanisms including photon generated carriers diffusion and carrier recombination processes have been discussed. These simulation methods and results should provide a useful tool for the evaluation and improvement of imaging power of InGaAs focal plane arrays.     

Ordering kinetics in bcc alloys with allowance for diffusion processes

The kinetics of inhomogeneous ordering in binary bcc substitution alloys of arbitrary stoichiometry is considered with the account for for diffusion processes. A system of kinetic equations describing the joint evolution of the occupancies of two sublattices of a bcc lattice is derived using a phenomenological approach. It is shown that, even within the mean-field approximation, this approach allows one to describe simultaneously the processes of establishment of the long-range order and diffusion of the alloy components. Numerical analysis of the obtained system of equations showed that quick onset of long-range order first occurs, which is followed by slow diffusion of the alloy components.     

Grain boundary sliding during ambient-temperature creep in hexagonal close-packed metals

Even at ambient temperature or less, below their 0.2% proof stresses all hexagonal close-packed metals and alloys show creep behaviour because they have dislocation arrays lying on a single slip system with no tangled dislocation inside each grain. In this case, lattice dislocations move without obstacles and pile-up in front of a grain boundary. Then these dislocations must be accommodated at the grain boundary to continue creep deformation. Atomic force microscopy revealed the occurrence of grain boundary sliding (GBS) in the ambient-temperature creep region. Lattice rotation of 5° was observed near grain boundaries by electron backscatter diffraction pattern analyses. Because of an extra low apparent activation energy of 20 kJ/mol, conventional diffusion processes are not activated. To accommodate these piled-up dislocations without diffusion processes, lattice dislocations must be absorbed by grain boundaries through a slip-induced GBS mechanism.     

Bounded fitness landscapes and the evolution of the linguistic diversity

We have recently introduced a simple spatial computer simulation model to study the evolution of the linguistic diversity. The model considers processes of selective geographic colonization, linguistic anomalous diffusion and mutation. In the approach, we ascribe to each language a fitness function which depends on the number of people that speak that language. Here, we extend the aforementioned model to examine the role of saturation of the fitness on the language dynamics. We found that the dependence of the linguistic diversity on the area after colonization displays a power law regime with a nontrivial exponent in very good agreement with the measured exponent associated with the actual distribution of languages on the Earth.     

多孔介质中流体流动及扩散的耦合格子Boltzmann模型

多孔介质中高Péclet数和大黏性比下混溶流体的流动和扩散广泛存在于二氧化碳驱油、化工生产等工业过程中.用数值方法对该问题进行研究时,关键在于如何正确描述高Péclet数和大黏性比下多孔介质内流体的行为.为此,提出了一种基于多松弛模型和格子动理模型的耦合格子Boltzmann模型.通过Chapman-Enskog分析,证明该模型能有效求解不可压Navier-Stokes方程和对流扩散方程.数值结果表明,该模型不仅具有二阶精度和良好的稳健性,而且对于高Péclet数和大黏性比的问题具有良好的数值稳定性,为模拟此类问题提供了有效工具.     

Numerical simulation of diffusion in a three-dimensional disordered lattice

The conductivity in the Anderson's model of a disordered diamond lattice is studied numerically via a direct simulation of the particle diffusion. In the vicinity of mobility edges the results reveal the continuous variation of diffusivity and the divergence of the amplitude fluctuations. The critical exponent for diffusivity t_w = 1.45 ± 0.1, obtained at fixed energy E = 0, seems to support the analogy with the classical percolation.