Exact solution of a Lévy walk model for anomalous heat transport

Transition from ballistic to diffusive heat transfer in a chain with breaks

The transition from a ballistic to a diffusive regime of heat transfer is studied using two models. The first model is a one-dimensional chain with bonds, capable of dissociation. Interparticle forces in the chain are harmonic for bond deformations below a critical value, corresponding to the dissociation, and zero above this value. A kinetic description of heat transfer in the chain is proposed using the second model, namely, a gas of noninteracting quasiparticles, reflecting from randomly occurring barriers. The motion of quasiparticles mimicks heat (energy) transfer in the chain, while the barriers mimic dissociated bonds. For the gas, a kinetic equation is derived and solved analytically. The solution demonstrates the transition from the ballistic regime at small times to the diffusive regime at large times. In the diffusive limit, the distance traveled by a heat obeys square-root asymptotics as in the case of classical diffusion. However, the shape of the fundamental solution for temperature differs from the Gaussian function and therefore the Fourier law is not satisfied. Two examples are considered to demonstrate that the presented kinetic model is in good qualitative agreement with the results of the numerical solution of the chain dynamics. The presented results show that bond dissociation is an important mechanism underlying the transition from ballistic to diffusive heat transfer in one-dimensional chains.

Nonequilibrium spin transport in integrable and nonintegrable classical spin chains

Anomalous transport in low dimensional spin chains is an intriguing topic that can offer key insights into the interplay of integrability and symmetry in many-body dynamics. Recent studies have shown that spin-spin correlations in spin chains, where integrability is either perfectly preserved or broken by symmetry-preserving interactions, fall in the Kardar-Parisi-Zhang (KPZ) universality class. Similarly, energy transport can show ballistic or diffusive-like behavior. Although such behavior has been studied under equilibrium conditions, no results on nonequilibrium spin transport in classical spin chains has been reported so far. In this work, we investigate both spin and energy transport in classical spin chains (integrable and nonintegrable) when coupled to two reservoirs at two different temperatures/magnetizations. In both the integrable case and the broken-integrability (but spin-symmetry preserving) case, we report anomalous scaling of spin current with system size (J^{s}∝L^{-μ}), with an exponent value that, within error bars, is close to the KPZ universality class value μ≈2/3. On the other hand, it is noteworthy that the energy current remains ballistic (J^{e}∝L^{-η} with η≈0) in the purely integrable case while there is departure from ballistic behavior (η>0) when integrability is broken regardless of spin-symmetry. We also present results on other interesting observables in the nonequilibrium steady state such as the spatial profiles of magnetization and energy, and spin-spin correlations.

One-dimensional Lévy quasicrystal

Space-fractional quantum mechanics (SFQM) is a generalization of the standard quantum mechanics when the Brownian trajectories in Feynman path integrals are replaced by Lévy flights. We introduce Lévy quasicrystal by discretizing the space-fractional Schrödinger equation using the Grünwald-Letnikov derivatives and adding on-site quasiperiodic potential. The discretized version of the usual Schrödinger equation maps to the Aubry-André (AA) Hamiltonian, which supports localization-delocalization transition even in one dimension. We find the similarities between Lévy quasicrystal and the AA model with power-law hopping, and show that the Lévy quasicrystal supports a delocalization-localization transition as one tunes the quasiperiodic potential strength and shows the coexistence of localized and delocalized states separated by mobility edge. Hence, a possible realization of SFQM in optical experiments should be a new experimental platform to test the predictions of AA models in the presence of power-law hopping.

Universal Kardar-Parisi-Zhang transient diffusion in nonequilibrium anharmonic chains

The well known nonlinear fluctuating hydrodynamics theory has grouped diffusions in anharmonic chains into two universality classes: one is the Kardar-Parisi-Zhang (KPZ) class for chains with either asymmetric potential or nonzero static pressure and the other is the Gaussian class for chains with symmetric potential at zero static pressure, such as Fermi-Pasta-Ulam-Tsingou (FPUT)-β chains. However, little is known of the nonequilibrium transient diffusion in anharmonic chains. Here, we reveal that the KPZ class is the only universality class for nonequilibrium transient diffusion, manifested as the KPZ scaling of the side peaks of momentum correlation (corresponding to the sound modes correlation), which was completely unexpected in equilibrium FPUT-β chains. The underlying mechanism is that the nonequilibrium soliton dynamics cause nonzero transient pressure so that the sound modes satisfy approximately the noisy Burgers equation, in which the collisions of solitons was proved to yield the KPZ dynamic exponent of the soliton dispersion. Therefore, the unexpected KPZ universality class is obtained in the nonequilibrium transient diffusion in FPUT-β chains and the corresponding carriers of nonequilibrium transient diffusion are attributed to solitons.

Infinite density and relaxation for Lévy walks in an external potential: Hermite polynomial approach

Lévy walks are continuous-time random-walk processes with a spatiotemporal coupling of jump lengths and waiting times. We here apply the Hermite polynomial method to study the behavior of LWs with power-law walking time density for four different cases. First we show that the known result for the infinite density of an unconfined, unbiased LW is consistently recovered. We then derive the asymptotic behavior of the probability density function (PDF) for LWs in a constant force field, and we obtain the corresponding qth-order moments. In a harmonic external potential we derive the relaxation dynamic of the LW. For the case of a Poissonian walking time an exponential relaxation behavior is shown to emerge. Conversely, a power-law decay is obtained when the mean walking time diverges. Finally, we consider the case of an unconfined, unbiased LW with decaying speed v(τ)=v_{0}/sqrt[τ]. When the mean walking time is finite, a universal Gaussian law for the position-PDF of the walker is obtained explicitly.

Numerical Solution to Anomalous Diffusion Equations for Levy Walks

The process of Levy random walks is considered in view of the constant velocity of a particle. A kinetic equation is obtained that describes the process of walks, and fractional differential equations are obtained that describe the asymptotic behavior of the process. It is shown that, in the case of finite and infinite mathematical expectation of paths, these equations have a completely different form. To solve the obtained equations, the method of local estimation of the Monte Carlo method is described. The solution algorithm is described and the advantages and disadvantages of the considered method are indicated.

Too Close to Integrable: Crossover from Normal to Anomalous Heat Diffusion

Energy transport in one-dimensional chains of particles with three conservation laws is generically anomalous and belongs to the Kardar-Parisi-Zhang dynamical universality class. Surprisingly, some examples where an apparent normal heat diffusion is found over a large range of length scales were reported. We propose a novel physical explanation of these intriguing observations. We develop a scaling analysis that explains how this may happen in the vicinity of an integrable limit, such as, but not only, the famous Toda model. In this limit, heat transport is mostly supplied by quasiparticles with a very large mean free path ℓ. Upon increasing the system size L, three different regimes can be observed: a ballistic one, an intermediate diffusive range, and, eventually, the crossover to the anomalous (hydrodynamic) regime. Our theoretical considerations are supported by numerical simulations of a gas of diatomic hard-point particles for almost equal masses and of a weakly perturbed Toda chain. Finally, we discuss the case of the perturbed harmonic chain, which exhibits a yet different scenario.

Lévy walk dynamics in an external harmonic potential

Lévy walks (LWs) are spatiotemporally coupled random-walk processes describing superdiffusive heat conduction in solids, propagation of light in disordered optical materials, motion of molecular motors in living cells, or motion of animals, humans, robots, and viruses. We here investigate a key feature of LWs-their response to an external harmonic potential. In this generic setting for confined motion we demonstrate that LWs equilibrate exponentially and may assume a bimodal stationary distribution. We also show that the stationary distribution has a horizontal slope next to a reflecting boundary placed at the origin, in contrast to correlated superdiffusive processes. Our results generalize LWs to confining forces and settle some longstanding puzzles around LWs.

Universality in the Onset of Superdiffusion in Lévy Walks

Anomalous dynamics in which local perturbations spread faster than diffusion are ubiquitously observed in the long-time behavior of a wide variety of systems. Here, the manner by which such systems evolve towards their asymptotic superdiffusive behavior is explored using the 1D Lévy walk of order 1<β<2. The approach towards superdiffusion, as captured by the leading correction to the asymptotic behavior, is shown to remarkably undergo a transition as β crosses the critical value β_{c}=3/2. Above β_{c}, this correction scales as |x|∼t^{1/2}, describing simple diffusion. However, below β_{c} it is instead found to remain superdiffusive, scaling as |x|∼t^{1/(2β-1)}. This transition is shown to be independent of the precise model details and is thus argued to be universal.

Lévy walks on finite intervals: A step beyond asymptotics

A Lévy walk of order β is studied on an interval of length L, driven out of equilibrium by different-density boundary baths. The anomalous current generated under these settings is nonlocally related to the density profile through an integral equation. While the asymptotic solution to this equation is known, its finite-L corrections remain unstudied despite their importance in the study of anomalous transport. Here a perturbative method for computing such corrections is presented and explicitly demonstrated for the leading correction to the asymptotic transport of a Lévy walk of order β=5/3, which represents a broad universal class of anomalous transport models. Surprisingly, many other physical problems are described by similar integral equations, to which the method introduced here can be directly applied.

Solitons as candidates for energy carriers in Fermi-Pasta-Ulam lattices

Currently, effective phonons (renormalized or interacting phonons) rather than solitary waves (for short, solitons) are regarded as the energy carriers in nonlinear lattices. In this work, by using the approximate soliton solutions of the corresponding equations of motion and adopting the Boltzmann distribution for these solitons, the average velocities of solitons are obtained and are compared with the sound velocities of energy transfer. Excellent agreements with the numerical results and the predictions of other existing theories are shown in both the symmetric Fermi-Pasta-Ulam-β lattices and the asymmetric Fermi-Pasta-Ulam-αβ lattices. These clearly indicate that solitons are suitable candidates for energy carriers in Fermi-Pasta-Ulam lattices. In addition, the root-mean-square velocity of solitons can be obtained from the effective phonons theory.

Energy Current Cumulants in One-Dimensional Systems in Equilibrium

A recent theory based on fluctuating hydrodynamics predicts that one-dimensional interacting systems with particle, momentum, and energy conservation exhibit anomalous transport that falls into two main universality classes. The classification is based on behavior of equilibrium dynamical correlations of the conserved quantities. One class is characterized by sound modes with Kardar-Parisi-Zhang scaling, while the second class has diffusive sound modes. The heat mode follows Lévy statistics, with different exponents for the two classes. Here we consider heat current fluctuations in two specific systems, which are expected to be in the above two universality classes, namely, a hard particle gas with Hamiltonian dynamics and a harmonic chain with momentum conserving stochastic dynamics. Numerical simulations show completely different system-size dependence of current cumulants in these two systems. We explain this numerical observation using a phenomenological model of Lévy walkers with inputs from fluctuating hydrodynamics. This consistently explains the system-size dependence of heat current fluctuations. For the latter system, we derive the cumulant-generating function from a more microscopic theory, which also gives the same system-size dependence of cumulants.

A violation of universality in anomalous Fourier's law

Since the discovery of long-time tails, it has been clear that Fourier's law in low dimensions is typically anomalous, with a size-dependent heat conductivity, though the nature of the anomaly remains puzzling. The conventional wisdom, supported by renormalization-group arguments and mode-coupling approximations within fluctuating hydrodynamics, is that the anomaly is universal in 1d momentum-conserving systems and belongs in the Lévy/Kardar-Parisi-Zhang universality class. Here we challenge this picture by using a novel scaling method to show unambiguously that universality breaks down in the paradigmatic 1d diatomic hard-point fluid. Hydrodynamic profiles for a broad set of gradients, densities and sizes all collapse onto an universal master curve, showing that (anomalous) Fourier's law holds even deep into the nonlinear regime. This allows to solve the macroscopic transport problem for this model, a solution which compares flawlessly with data and, interestingly, implies the existence of a bound on the heat current in terms of pressure. These results question the renormalization-group and mode-coupling universality predictions for anomalous Fourier's law in 1d, offering a new perspective on transport in low dimensions.

Single integrodifferential wave equation for a Lévy walk

We derive the single integrodifferential wave equation for the probability density function of the position of a classical one-dimensional Lévy walk with continuous sample paths. This equation involves a classical wave operator together with memory integrals describing the spatiotemporal coupling of the Lévy walk. It is valid at all times, not only in the long time limit, and it does not involve any large-scale approximations. It generalizes the well-known telegraph or Cattaneo equation for the persistent random walk with the exponential switching time distribution. Several non-Markovian cases are considered when the particle's velocity alternates at the gamma and power-law distributed random times. In the strong anomalous case we obtain the asymptotic solution to the integrodifferential wave equation. We implement the nonlinear reaction term of Kolmogorov-Petrovsky-Piskounov type into our equation and develop the theory of wave propagation in reaction-transport systems involving Lévy diffusion.

Numerical test of hydrodynamic fluctuation theory in the Fermi-Pasta-Ulam chain

Recent work has developed a nonlinear hydrodynamic fluctuation theory for a chain of coupled anharmonic oscillators governing the conserved fields, namely, stretch, momentum, and energy. The linear theory yields two propagating sound modes and one diffusing heat mode, all three with diffusive broadening. In contrast, the nonlinear theory predicts that, at long times, the sound mode correlations satisfy Kardar-Parisi-Zhang scaling, while the heat mode correlations have Lévy-walk scaling. In the present contribution we report on molecular dynamics simulations of Fermi-Pasta-Ulam chains to compute various spatiotemporal correlation functions and compare them with the predictions of the theory. We obtain very good agreement in many cases, but also some deviations.

Anomalous heat diffusion

Consider anomalous energy spread in solid phases, i.e., <Δx(2)(t)>E≡∫(x-<x>E)(2)ρE(x,t)dx∝t(β), as induced by a small initial excess energy perturbation distribution ρE(x,t=0) away from equilibrium. The second derivative of this variance of the nonequilibrium excess energy distribution is shown to rigorously obey the intriguing relation d(2)<Δx(2)(t)>E/dt2=2CJJ(t)/(kBT(2)c), where CJJ(t) equals the thermal equilibrium total heat flux autocorrelation function and c is the specific volumetric heat capacity. Its integral assumes a time-local Helfand-like relation. Given that the averaged nonequilibrium heat flux is governed by an anomalous heat conductivity, the energy diffusion scaling determines a corresponding anomalous thermal conductivity scaling behavior.

Space-time velocity correlation function for random walks

Space-time correlation functions constitute a useful instrument from the research toolkit of continuous-media and many-body physics. Here we adopt this concept for single-particle random walks and demonstrate that the corresponding space-time velocity autocorrelation functions reveal correlations which extend in time much longer than estimated with the commonly employed temporal correlation functions. A generic feature of considered random-walk processes is an effect of velocity echo identified by the existence of time-dependent regions where most of the walkers are moving in the direction opposite to their initial motion. We discuss the relevance of the space-time velocity correlation functions for the experimental studies of cold atom dynamics in an optical potential and charge transport on micro- and nanoscales.

Anomalous fluctuations of currents in Sinai-type random chains with strongly correlated disorder

We study properties of a random walk in a generalized Sinai model, in which a quenched random potential is a trajectory of a fractional Brownian motion with arbitrary Hurst parameter H, 0<H<1, so that the random force field displays strong spatial correlations. In this case, the disorder-average mean-square displacement grows in proportion to log(2/H)(n), n being time. We prove that moments of arbitrary order k of the steady-state current J(L) through a finite segment of length L of such a chain decay as L(-(1-H)), independently of k, which suggests that despite a logarithmic confinement the average current is much higher than its Fickian counterpart in homogeneous systems. Our results reveal a paradoxical behavior such that, for fixed n and L, the mean-square displacement decreases when one varies H from 0 to 1, while the average current increases. This counterintuitive behavior is explained via an analysis of representative realizations of disorder.