Algebraic Tail Research Articles

Testing the instanton approach to the large amplification limit of a diffraction–amplification problem

Abstract The validity of the instanton analysis approach is tested numerically in the case of the diffraction-amplification problem ∂zψ−(i/2m)∂2 x2 ψ=g│S│2 for ln U>>1, where U=│ψ(0,L)│2. Here, S(x,z) is a complex Gaussian random field, z and x respectively are the axial and transverse coordinates, with 0≤z≤L, and both m≠0 and g>0 are real parameters. We consider a class of S, called the `one-max class', for which we devise a specific biased sampling procedure. As an application, p(U), the probability distribution of U, is obtained down to values less than 10-2270 in the far right tail. We find that the agreement of our numerical results with the instanton analysis predictions in Mounaix (2023 J. Phys. A: Math. Theor. 56 305001) is remarkable. Both the predicted algebraic tail of p(U) and concentration of the realizations of S onto the leading instanton are clearly confirmed, which validates the instanton analysis numerically in the large ln U limit for S in the one-max class.

Precise rates of propagation in reaction–diffusion equations with logarithmic Allee effect

This paper focuses on propagation phenomena in reaction–diffusion equations with a weakly degenerate monostable nonlinearity. The kind of reaction term we consider can be seen as an intermediate between the classical logistic one (or Fisher–KPP) and more usual power laws that usually model Allee effects. We investigate the effect of the decay rate of the initial data on the propagation rate. When the tail of the initial data is sub-exponential, both finite speed propagation and acceleration may happen. We derive the exact separation between the two situations. When the initial data is sub-exponentially unbounded, acceleration unconditionally occurs. Estimates for the locations of the level sets are expressed in terms of the decay of the initial data. In addition, sharp exponents of acceleration for initial data with sub-exponential and algebraic tails are given. Numerical simulations are presented to illustrate the above findings.

Stochastic modeling of injection induced seismicity based on the continuous time random walk model

The spatiotemporal evolution of earthquakes induced by fluid injections into the subsurface can be erratic owing to the complexity of the physical process. To effectively mitigate the associated hazard and to draft appropriate regulatory strategies, a detailed understanding of how induced seismicity may evolve is needed. In this work, we build on the well-established continuous-time random walk (CTRW) theory to develop a purely stochastic framework that can delineate the essential characteristics of this process. We use data from the 2003 and 2012 hydraulic stimulations in the Cooper Basin geothermal field that induced thousands of microearthquakes to test and demonstrate the applicability of the model. Induced seismicity in the Cooper Basin shows all the characteristics of subdiffusion, as indicated by the fractional order power-law growth of the mean square displacement with time and broad waiting-time distributions with algebraic tails. We further use an appropriate master equation and the time-fractional diffusion equation to map the spatiotemporal evolution of seismicity. The results show good agreement between the model and the data regarding the peak earthquake concentration close to the two injection wells and the stretched exponential relaxation of seismicity with distance, suggesting that the CTRW model can be efficiently incorporated into induced seismicity forecasting.

Critical dynamics of long-range quantum disordered systems.

Long-range hoppings in quantum disordered systems are known to yield quantum multifractality, the features of which can go beyond the characteristic properties associated with an Anderson transition. Indeed, critical dynamics of long-range quantum systems can exhibit anomalous dynamical behaviors distinct from those at the Anderson transition in finite dimensions. In this paper, we propose a phenomenological model of wave packet expansion in long-range hopping systems. We consider both their multifractal properties and the algebraic fat tails induced by the long-range hoppings. Using this model, we analytically derive the dynamics of moments and inverse participation ratios of the time-evolving wave packets, in connection with the multifractal dimension of the system. To validate our predictions, we perform numerical simulations of a Floquet model that is analogous to the power law random banded matrix ensemble. Unlike the Anderson transition in finite dimensions, the dynamics of such systems cannot be adequately described by a single parameter scaling law that solely depends on time. Instead, it becomes crucial to establish scaling laws involving both the finite size and the time. Explicit scaling laws for the observables under consideration are presented. Our findings are of considerable interest towards applications in the fields of many-body localization and Anderson localization on random graphs, where long-range effects arise due to the inherent topology of the Hilbert space.

Anomalous Diffusion in the Long-Range Haken-Strobl-Reineker Model.

We analyze the propagation of excitons in a d-dimensional lattice with power-law hopping ∝1/r^{α} in the presence of dephasing, described by a generalized Haken-Strobl-Reineker model. We show that in the strong dephasing (quantum Zeno) regime the dynamics is described by a classical master equation for an exclusion process with long jumps. In this limit, we analytically compute the spatial distribution, whose shape changes at a critical value of the decay exponent α_{cr}=(d+2)/2. The exciton always diffuses anomalously: a superdiffusive motion is associated to a Lévy stable distribution with long-range algebraic tails for α≤α_{cr}, while for α>α_{cr} the distribution corresponds to a surprising mixed Gaussian profile with long-range algebraic tails, leading to the coexistence of short-range diffusion and long-range Lévy flights. In the many-exciton case, we demonstrate that, starting from a domain-wall exciton profile, algebraic tails appear in the distributions for any α, which affects thermalization: the longer the hopping range, the faster equilibrium is reached. Our results are directly relevant to experiments with cold trapped ions, Rydberg atoms, and supramolecular dye aggregates. They provide a way to realize an exclusion process with long jumps experimentally.

Observation of 1/k^{4}-Tails after Expansion of Bose-Einstein Condensates with Impurities.

We measure the momentum density in a Bose-Einstein condensate (BEC) with dilute spin impurities after an expansion in the presence of interactions. We observe tails decaying as 1/k^{4} at large momentum k in the condensate and in the impurity cloud. These algebraic tails originate from the impurity-BEC interaction, but their amplitudes greatly exceed those expected from two-body contact interactions at equilibrium in the trap. Furthermore, in the absence of impurities, such algebraic tails are not found in the BEC density measured after the interaction-driven expansion. These results highlight the key role played by impurities when present, a possibility that had not been considered in our previous work [Chang etal., Phys. Rev. Lett. 117, 235303 (2016)PRLTAO0031-900710.1103/PhysRevLett.117.235303]. Our measurements suggest that these unexpected algebraic tails originate from the nontrivial dynamics of the expansion in the presence of impurity-bath interactions.

Universality and Control of Fat Tails

Motivated by applications in hydrodynamics and networks of thermostatically-control loads in buildings we study control of linear dynamical systems driven by additive and also multiplicative noise of a general position. Utilizing mathematical theory of stochastic multiplicative processes we present a universal way to estimate fat, algebraic tails of the state vector probability distributions. This prompts us to introduce and analyze mean-q-power stability criterion, generalizing the mean-square stability criterion, and then juxtapose it to other tools in control.

Nonequilibrium relaxation of a trapped particle in a near-critical Gaussian field.

We study the nonequilibrium relaxational dynamics of a probe particle linearly coupled to a thermally fluctuating scalar field and subject to a harmonic potential, which provides a cartoon for an optically trapped colloid immersed in a fluid close to its bulk critical point. The average position of the particle initially displaced from the position of mechanical equilibrium is shown to feature long-time algebraic tails as the critical point of the field is approached, the universal exponents of which are determined in arbitrary spatial dimensions. As expected, this behavior cannot be captured by adiabatic approaches which assumes fast field relaxation. The predictions of the analytic, perturbative approach are qualitatively confirmed by numerical simulations.

The Boltzmann Equation for Uniform Shear Flow

The uniform shear flow for the rarefied gas is governed by the time-dependent spatially homogeneous Boltzmann equation with a linear shear force. The main feature of such flow is that the temperature may increase in time due to the shearing motion that induces viscous heat and the system becomes far from equilibrium. For Maxwell molecules, we establish the unique existence, regularity, shear-rate-dependent structure and non-negativity of self-similar profiles for any small shear rate. The non-negativity is justified through the large time asymptotic stability even in spatially inhomogeneous perturbation framework, and the exponential rates of convergence are also obtained with the size proportional to the second order shear rate. The analysis supports the numerical result that the self-similar profile admits an algebraic high-velocity tail that is the key difficulty to overcome in the proof.

Asymptotic spreading of a time periodic diffusion equation with degenerate monostable nonlinearity

This paper deals with the asymptotic spreading for a class of diffusion equations with degenerate monostable nonlinearity, where the initial value is asymptotically front-like. By the method of squeezing technique and super- and sub-solutions, we prove that the speed of asymptotic spreading may be finite or infinite, which depends on the degeneracy of nonlinearity as well as the decaying rate of the initial value. Different from the non-degenerate case, the traveling wave solution with the critical speed is globally asymptotically stable if the initial value has light tails, while an acceleration propagation may occur with the initial value having heavy tails in the degenerate case. Furthermore, if the initial value has heavy tails but lighter than algebraic, then the speed of the asymptotic spreading is finite, and the balance between finite or infinite speed is established if the initial value has algebraic tails. Our results reflect that the degeneracy may be harmful to the acceleration propagation.

Wigner function for noninteracting fermions in hard-wall potentials

The Wigner function $W_N({\bf x}, {\bf p})$ is a useful quantity to characterize the quantum fluctuations of an $N$-body system in its phase space. Here we study $W_N({\bf x}, {\bf p})$ for $N$ noninteracting spinless fermions in a $d$-dimensional spherical hard box of radius $R$ at temperature $T=0$. In the large $N$ limit, the local density approximation (LDA) predicts that $W_N({\bf x}, {\bf p}) \approx 1/(2 \pi \hbar)^d$ inside a finite region of the $({\bf x}, {\bf p})$ plane, namely for $|{\bf x}| < R$ and $|{\bf p}| < k_F$ where $k_F$ is the Fermi momentum, while $W_N({\bf x}, {\bf p})$ vanishes outside this region, or "droplet", on a scale determined by quantum fluctuations. In this paper we investigate systematically, in this quantum region, the structure of the Wigner function along the edge of this droplet, called the Fermi surf. In one dimension, we find that there are three distinct edge regions along the Fermi surf and we compute exactly the associated nontrivial scaling functions in each regime. We also study the momentum distribution $\hat \rho_N(p)$ and find a striking algebraic tail for very large momenta $\hat \rho_N(p) \propto 1/p^4$, well beyond $k_F$, reminiscent of a similar tail found in interacting quantum systems (discussed in the context of Tan's relation). We then generalize these results to higher $d$ and find, remarkably, that the scaling function close to the edge of the box is universal, i.e., independent of the dimension~$d$.

Instantons for rare events in heavy-tailed distributions

Large deviation theory and instanton calculus for stochastic systems are widely used to gain insight into the evolution and probability of rare events. At its core lies the fact that rare events are, under the right circumstances, dominated by their least unlikely realization. Their computation through a saddle-point approximation of the path integral for the corresponding stochastic field theory then reduces an inefficient stochastic sampling problem into a deterministic optimization problem: finding the path of smallest action, the instanton. In the presence of heavy tails, though, standard algorithms to compute the instanton critically fail to converge. The reason for this failure is the divergence of the scaled cumulant generating function (CGF) due to a non-convex large deviation rate function. We propose a solution to this problem by ‘convexifying’ the rate function through a nonlinear reparametrization of the observable, which allows us to compute instantons even in the presence of super-exponential or algebraic tail decay. The approach is generalizable to other situations where the existence of the CGF is required, such as exponential tilting in importance sampling for Monte-Carlo algorithms. We demonstrate the proposed formalism by applying it to rare events in several stochastic systems with heavy tails, including extreme power spikes in fiber optics induced by soliton formation.

Kink-antikink collisions and multi-bounce resonance windows in higher-order field theories

We study collisions of coherent structures in higher-order field-theoretic models, such as the ϕ8,ϕ10 and ϕ12 ones. The main distinguishing feature, of the example models considered herein, is that the collision arises due to the long-range interacting algebraic tails of these solitary waves. We extend the approach to suitably initialize the relevant kinks, in the additional presence of finite initial velocity, in order to minimize the dispersive wave radiation potentially created by their slow spatial decay. We find that, when suitably initialized, these models still feature the multi-bounce resonance windows earlier found in models in which the kinks bear exponential tails, such as the ϕ4 and ϕ6 field theories among others. Also present is the self-similar structure of the associated windows with three- and more-bounce windows at the edges of two- and lower-bounce ones. Moreover, phenomenological but highly accurate (and predictive), scaling relations are derived for the dependence of the time between consecutive collisions and, e.g., the difference in kinetic energy between the incoming one and the critical one for one-bounces. Such scalings are traced extensively over two-bounce collision windows throughout the three models, hinting at the possibility of an analytical theory in this direction.

Edge states of the long-range Kitaev chain: An analytical study

We analyze the properties of the edge states of the one-dimensional Kitaev model with long-range anisotropic pairing and tunneling. Tunneling and pairing are assumed to decay algebraically with exponents $\ensuremath{\alpha}$ and $\ensuremath{\beta}$, respectively, and $\ensuremath{\alpha},\ensuremath{\beta}&gt;1$. We determine analytically the decay of the edge modes. We show that the decay is exponential for $\ensuremath{\alpha}=\ensuremath{\beta}$ and when the coefficients scaling tunneling and pairing terms are equal. Otherwise, the decay is exponential at sufficiently short distances and then algebraic at the asymptotics. We show that the exponent of the algebraic tail is determined by the smallest exponent between $\ensuremath{\alpha}$ and $\ensuremath{\beta}$. Our predictions are in agreement with numerical results found by exact diagonalization and reported in the literature.

Dynamics of a noninteracting colloidal fluid in a quenched Gaussian random potential: a time-reversal-symmetry-preserving field-theoretic approach

We develop a field-theoretic perturbation method preserving the fluctuation–dissipation relation (FDR) for the dynamics of the density fluctuations of a noninteracting colloidal gas plunged in a quenched Gaussian random field. It is based on an expansion about the Brownian noninteracting gas and can be considered and justified as a low-disorder or high-temperature expansion. The first-order bare theory yields the same memory integral as the mode-coupling theory (MCT) developed for (ideal) fluids in random environments, apart from the bare nature of the correlation functions involved. It predicts an ergodic dynamical behavior for the relaxation of the density fluctuations, in which the memory kernels and correlation functions develop long-time algebraic tails. An FDR-consistent renormalized theory is also constructed from the bare theory. It is shown to display a dynamic ergodic–nonergodic transition similar to the one predicted by the MCT at the level of the density fluctuations, but, at variance with the MCT, the transition does not fully carry over to the self-diffusion, which always reaches normal diffusive behavior at long time, in agreement with known rigorous results.

Nonlocal emergent hydrodynamics in a long-range quantum spin system

Generic short-range interacting quantum systems with a conserved quantity exhibit universal diffusive transport at late times. We employ non-equilibrium quantum field theory and semi-classical phase-space simulations to show how this universality is replaced by a more general transport process in a long-range XY spin chain at infinite temperature with couplings decaying algebraically with distance as $r^{-\alpha}$. While diffusion is recovered for $\alpha>1.5$, longer-ranged couplings with $0.5<\alpha\leq 1.5 $ give rise to effective classical L\'evy flights; a random walk with step sizes drawn from a distribution with algebraic tails. We find that the space-time dependent spin density profiles are self-similar, with scaling functions given by the stable symmetric distributions. As a consequence, for $0.5<\alpha\leq1.5$ autocorrelations show hydrodynamic tails decaying in time as $t^{-1/(2\alpha-1)}$ and linear-response theory breaks down. Our findings can be readily verified with current trapped ion experiments.

Optimal in situ electromechanical sensing of molecular species

We investigate protocols for optimal molecular detection with electromechanical nanoscale sensors under ambient conditions. Our models are representative of suspended graphene nanoribbons, which due to their piezoelectric and electronic properties provide responsive and versatile sensors. In particular, we analytically account for the corrections in the electronic transmission function and signal-to-noise ratio originating in environmental perturbations, such as thermal fluctuations and solvation effects. We also investigate the role of the sampling time in the current statistics. As a result, we formulate a protocol for optimal sensing based on the modulation of the Fermi level at a fixed bias and provide approximate forms for the current, linear susceptibility, and current fluctuations. We show how the algebraic tails in the thermally broadened transmission function affect the behavior of the signal-to-noise ratio and optimal sensing. These results provide further insights into the operation of graphene deflectometers and other techniques for electromechanical sensing.

Microscopic Slip Boundary Conditions in Unsteady Fluid Flows.

An algebraic tail in the Green-Kubo integral for the solid-fluid friction coefficient hampers its use in the determination of the slip length. A simple theory for discrete nonlocal hydrodynamics near parallel solid walls with extended friction forces is given. We explain the origin of the algebraic tail and give a solution of the plateau problem in the Green-Kubo expressions. We derive the slip boundary condition with a microscopic expression for the slip length and the hydrodynamic wall position, and assess it through simulations of an unsteady plug flow.

Fractional Hydrodynamic Memory and Superdiffusion in Tilted Washboard Potentials.

Diffusion in tilted washboard potentials can paradoxically exceed free normal diffusion. The effect becomes much stronger in the underdamped case due to inertial effects. What happens upon inclusion of usually neglected fractional hydrodynamics memory effects (Basset-Boussinesq frictional force), which result in a heavy algebraic tail of the velocity autocorrelation function of the potential-free diffusion making it transiently superdiffusive? Will a giant enhancement of diffusion become even stronger, and the transient superdiffusion last even longer? These are the questions that we answer in this Letter based on an accurate numerical investigation. We show that a resonancelike enhancement of normal diffusion becomes indeed much stronger and sharper. Moreover, a long-lasting transient regime of superdiffusion, including Richardson-like diffusion, ⟨δx^{2}(t)⟩∝t^{3} and ballistic supertransport, ⟨δx(t)⟩∝t^{2}, is revealed.

Persistent Anti-Correlations in Brownian Dynamics Simulations of Dense Colloidal Suspensions Revealed by Noise Suppression.

Transport properties of a hard-sphere colloidal fluid are investigated by Brownian dynamics simulations. We implement a novel algorithm for the time-dependent velocity-autocorrelation function (VACF) essentially eliminating the noise of the bare random motion. The measured VACF reveals persistent anti-correlations manifested by a negative algebraic power-law tail t^{-5/2} at all densities. At small packing fractions the simulations fully agree with the analytic low-density prediction, yet the amplitude of the tail becomes dramatically suppressed as the packing fraction is increased. The mode-coupling theory of the glass transition provides a qualitative explanation for the strong variation in terms of the static compressibility as well as the slowing down of the structural relaxation.