Enhancedt −3/2 long-time tail for the stress-stress time correlation function

Nonequilibrium molecular dynamics is used to calculate the spectrum of shear viscosity for a Lennard-Jones fluid. The calculated zero-frequency shear viscosity agrees well with experimental argon results for the two state points considered. The low-frequency behavior of shear viscosity is dominated by an ^1/2 cusp. Analysis of the form of this cusp reveals that the stress-stress time correlation function exhibits at ^–3/2 long-time tail. It is shown that for the state points studied, the amplitude of this long-time tail is between 12 and 150 times larger than what has been predicted theoretically. If the low-frequency results are truly asymptotic, they imply that the cross and potential contributions to the Kubo-Green integrand for shear viscosity exhibit at ^–3/2 long-time tail. This result contradicts the established theory of such processes.     

The velocity correlation function for the Lorentz gas

The results of variational solutions of the repeated ring and self-consistent repeated ring equations for the two-and three-dimensional overlapping Lorentz gas (LG), as formulated in a previous report, are presented. Calculations of the full velocity correlation function (VCF) for the 2D LG, including long-time tails, are compared with those from molecular dynamics. The trial functions chosen lead to predictions for the long-time tails that improve as the density of the scatterers is increased. At a value of 0.24 for* (= ², where is the density and the radius of scatterers), the self-consistent amplitudes of the long-time tail are within 40% of the molecular dynamics. A limited number of 3D results for the short-time behavior of the repeated ring VCF are presented. The 3D solutions agree with the molecular dynamics to within 10%.     

Long-time behavior of the angular velocity autocorrelation function for a fluid of rough spheres

According to hydrodynamical and mode-coupling theories, the angular velocity autocorrelation function decays at long times as ₀(t/10^–14 sec)^–5/2. For rough spheres under the conditions reported here, the quantity ₀ is predicted to be 262. The molecular dynamics studies presented here yield a long-time tail of the form 230(t/10^–14 sec)^–2.38. The disagreement between theory and computer results probably arises from statistical error intrinsic to the computations.The authors are indebted to the National Science Foundation and the Computer Center of the University of Minnesota for financial support of the research reported here.     

Exact time correlation function for a nonlinearly coupled vibrational system

The time correlation function (t)=Re<[c(t), c (0)]>, which is related to the dipole spectrum and is the main focus of quantum molecular time scale generalized Langevin equation theory, is calculated for the Hamiltonian system in which a single oscillator is coupled by a nonlinear Davydov term to a chain of oscillators comprising a phonon heat bath. An exact expression for (t) is obtained. At long times we find that the time correlation function decays as a small power law atT=0K, but switches to exponential decay at higher temperature. This is a new result and bears on the long-standing issue of the existence of long-time tails.     

The goodness of ergodic adiabatic invariants

For a slowly time-dependent Hamiltonian system exhibiting chaotic motion that ergodically covers the energy surface, the phase space volume enclosed inside this surface is an adiabatic invariant. In this paper we examine, both numerically and theoretically, how the error in this ergodic adiabatic invariant scales with the slowness of the time variation of the Hamiltonian. It is found that under certain circumstances, the error is diffusive and scales likeT ^–1/2, whereT is the characteristic time over which the Hamiltonian changes. On the other hand, for other cases (where motion in the Hamiltonian has a long-time 1/t tail in a certain correlation function), the error scales like [T ^–1 ln(T)]^1/2. Both of these scalings are verified by numerical experiments. In the situation where invariant tori exist amid chaos, the motion may not be fully ergodic on the entire energy surface. The ergodic adiabatic invariant may still be useful in this case and the circumstances under which this is so are investigated numerically (in particular, the islands have to be small enough).     

Long-time behavior of Navier-Stokes flow on a two-dimensional torus excited by an external sinusoidal force

In this paper we study the Navier-Stokes flow on the two-dimensional torusS ¹ ×S ¹ excited by the external force (k ² sinky, 0) and find the long-time behavior for the flow starting from some states, whereS ¹=[0,2](mod 2). Especially for the casek=2, it follows from an analysis and computation that the Navier-Stokes flow with the initial state cos(mx+ny) or sin(mx+ny) will likely evolve through at most one step bifurcation to either a steady-state solution or a time-dependent periodic solution for any Reynolds number and integersm andn.     

Long-Time tails in a random diffusion model

Letw = {w(x)xZ^d} be a positive random field with i.i.d. distribution. Given its realization, letX _t be the position at timet of a particle starting at the origin and performing a simple random walk with jump rate w^–1(X_t). The processX={X _t:t0} combined withw on a common probability space is an example of random walk in random environment. We consider the quantities _t =(d/dt) E (X _t ² –M ^–1 t and _t(w) = (d/dt)E_w(X _t ² – M^– 1t). Here E_w. is expectation overX at fixedw and E = E_w (dw) is the expectation over bothX andw. We prove the following long-time tail results: (1) lim_t t^d/2_t= V²M^d/2–3(d/2)^d/2 and (2) lim_t t^d/4 _st(w)= Z_s weakly in path space, with {Z_s:s>0} the Gaussian process with EZ_s=0 and EZ_rZ_s= V²M^d/2–4(d)^d/2 (r + s)^–d/2. HereM and V² are the mean and variance of w(0) under . The main surprise is that fixingw changes the power of the long-time tail fromd/2 tod/4. Since , with ₀ the stationary measure for the environment process, our result (1) exhibits a long-time tail in an equilibrium autocorrelation function.     

Effect of disorder on relaxation in the one-dimensional Glauber model

We study the long-time relaxation of magnetization in a disordered linear chain of Ising spins from an initially aligned state. The coupling constants are ferromagnetic and nearest-neighbor only, taking valuesJ ₀ andJ ₁ with probabilitiesp and 1–p, respectively. The time evolution of the system is governed by the Glauber master equation. It is shown that for large timest, the magnetizationM(t) varies as [exp(–₀ t](t), where ₀ is a function of the stronger bond strengthJ ₀ only, and (t) decreases slower than an exponential. For very long times, we find that ln (t) varies as –t ^1/3. For low enough temperatures, there is an intermediate time regime when ln (t) varies as –t ^1/2. The results can be extended to more general probability distributions of ferromagnetic coupling constants, assuming thatM(t) can only increase if any bond in the chain is strengthened. If the coupling constants have a continuous distribution in which the probability density varies as a power law near some maximum valueJ ₀, we find that ln (t) varies as –t ^1/3(lnt)^2/3 for large times.     

Lévy walks for turbulence: A numerical study

We present a numerical study of enhanced diffusion, for which the mean-squared displacement follows asymptotically r ²(t) t , > 1. We simulate continuous time random walks with waiting-time distributions which couple the spatial and temporal parameters; this gives rise to Lévy-walks. Our results confirm the theoretically predicted long-time behavior and demonstrate its temporal regime of validity. Furthermore, the simulations document the appearance of (parameter-dependent) transitions between regular and enhanced diffusion regimes.     

Asymptotic time behavior of correlation functions. II. Kinetic and potential terms

On the basis of the mode-coupling theory we obtain the long-time behavior t ^–d/2 for the kinetic, potential, and cross-terms in the Green-Kubo integrands, expressed completely in terms of transport coefficients and thermodynamic quantities. All two-mode amplitudes are explicitly evaluated in terms of measurable quantities such as specific heats, thermal expansion coefficients, etc.     

Nonequilibrium transitions driven by external dichotomous noise

The stationary probability densityP _s for a class of nonlinear one-dimensional models driven by a dichotomous Markovian process (DMP)I _t , can be calculated explicitly. For the specific case of the Stratonovich model, x=ax –-x ³ +I _t x, the qualitative shape ofP _s and its support is discussed in the whole parameter region. The location of the maxima ofP _s shows a behavior similar to order parameters in continuous phase transitions. The possibility of a noiseinduced change from continuous to a discontinuous transition in an extended model, in which the DMP couples also to the cubic term, is discussed. The time-dependent moments x_t ⁿ can be represented as an infinite series of terms, which are determined by a recursion formula. For negative even moments the series terminates and the long-time behavior can be obtained analytically. As a function of the physical parameters, qualitative changes of this behavior may occur which can be partially related to the behavior ofP _s . All results reproduce those for Gaussian white noise in the corresponding limit. The influence of the finite correlation time and the discreteness of the space of states of the DMP are discussed. An extensive list of references is contained in U. Behn, K. Schiele, and A. Teubel,Wiss. Z. Karl-Marx-Univ. Leipzig, Mathem.-Naturwiss. R. 34:602 (1985).Contribution to the symposium on the Statistical Mechanics of Phase Transitions —Mathematical and Physical Aspects, Tebo, Czechoslovakia, September 1–6, 1986.     

Long-time behavior of the nonlocal shear viscosity of a one-component plasma: A microscopic approach

The nonlocal shear viscosity(t) of a classical one-component plasma is shown to have anoscillatory long-time tail. This result is obtained on the basis of amicroscopic theory which does not rely on expansions in a small parameter such as the plasma expansion parameter. Our major approximation is the restriction to the coupling oftwo hydrodynamic propagators in the computation of the long-time behavior of the transport matrix. The Coulomb divergence is correctly accounted for, while the nonanalyticities of both the plasma parameter and gradient expansions are discussed at the level of the kinetic as well as the hydrodynamic equations.     

Generalized q-Modified Laguerre Functions

In this paper we define a new q-special function A _n (x, b, c; q). The new function is a generalization of the q-Laguerre function and the Stieltjes–Wigert function. We deduced all the properties of the function A _n (x, b, c; q). Finally, lim _q1 A _n ((1 – q)x, –, 1;q) gives L _n ^(,)(x,q), which is a -modification of the ordinary Laguerre function.     

Bond percolation in two and three dimensions: Numerical evaluation of time-dependent transport properties

For random walks on two- and three-dimensional cubic lattices, numerical results are obtained for the static,D(), and time-dependent diffusion coefficientD(t), as well as for the velocity autocorrelation function (VACF). The results cover all times and include linear and quadratic terms in the density expansions. Within the context of kinetic theory this is the only model in two and three dimensions for which the time-dependent transport properties have been calculated explicitly, including the long-time tails.     

A long-time tail for random walk in random scenery

Consider a simple random walk on ^d whose sites are colored black or white independently with probabilityq, resp. 1–q. Walk and coloring are independent. Letn _k be the number of steps by the walk between itskth and (k+1) th visits to a black site (i.e., the length of itskth white run), and let _k =E(n _k )–q ^–1. Our main result is a proof that (*) lim _k k ^d/2 _k = (1 –q)q ^{d/2 – 2}(d/2) ^d/2. Since it is known thatq ^{– 1} _k =E(n ₁ n _{k + 1} B) –E(n ₁ B)E(n _{k + 1} B), withB the event that the origin is black, (*) exhibits a long-time tail in the run length autocorrelation function. Numerical calculations of _k (1k100) ind=1, 2, and 3 show that there is an oscillatory behavior of _k for smallk. This damps exponentially fast, following which the power law sets in fairly rapidly. We prove that if the coloring is not independent, but is convex in the sense of FKG, then the decay of _k cannot be faster than (*).     

Contractive Metrics for a Boltzmann Equation for Granular Gases: Diffusive Equilibria

We quantify the long-time behavior of a system of (partially) inelastic particles in a stochastic thermostat by means of the contractivity of a suitable metric in the set of probability measures. Existence, uniqueness, boundedness of moments and regularity of a steady state are derived from this basic property. The solutions of the kinetic model are proved to converge exponentially as t to this diffusive equilibrium in this distance metrizing the weak convergence of measures. Then, we prove a uniform bound in time on Sobolev norms of the solution, provided the initial datum has a finite norm in the corresponding Sobolev space. These results are then combined, using interpolation inequalities, to obtain exponential convergence to the diffusive equilibrium in the strong L¹-norm, as well as various Sobolev norms.     

Time-dependent correlation functions for random walks on bond disordered cubic lattices

For hopping models on cubic lattices with a fractionc of impurity bonds, time-dependent transport properties and correlation functions (long-time tails) are calculated through a systematicc-expansion (in the percolation literature referred to as high-density expansion), using a method developed in an earlier paper. The time-dependent diffusion coefficient, velocity autocorrelation function (VACF), and Burnett functions are calculated exact toO(c) for allt, and exact toO(c ₂ ) for long times only. A comparison is made with the results of the effective medium approximation, and numerical results are given for the square lattice.     

New universal short-time scaling behaviour of critical relaxation processes

We study the critical relaxation properties of Model A (purely dissipative relaxation) starting from a macroscopically prepared initial state characterised by non-equilibrium values for order parameter and correlations. Using a renormalisation group approach we observe that even (macroscopically)early stages of the relaxation process display universal behaviour governed by a new, independent initial slip exponent. For large times, the system crosses over to the well-known long-time relaxation behaviour.The new exponent is calculated toO(²) in =4–d, whered is the spatial dimension of the system. The initial slip scaling form of general correlation and response functions as well as the order parameter is derived, exploiting a short-time operator expansion. The leading scaling behaviour is determined by initial states with sharp values of the order parameter. Non-vanishing correlations generate corrections to scaling.     

Measurement of the diffractive structure function in deep inelastic scattering at HERA

This paper presents an analysis of the inclusive properties of diffractive deep inelastic scattering events produced inep interactions at HERA. The events are characterised by a rapidity gap between the outgoing proton system and the remaining hadronic system. Inclusive distributions are presented and compared with Monte Carlo models for diffractive processes. The data are consistent with models where the pomeron structure function has a hard and a soft contribution. The diffractive structure function is measured as a function ofx , the momentum fraction lost by the proton, of , the momentum fraction of the struck quark with respect tox , and ofQ ² in the range 6.3·10^–4x <>^–2, 0.1<0.8 and=">Q ²<100>². The dependence is consistent with the formx wherea=1.30±0.08(stat) _–0.14 ^+0.08 (sys) in all bins of andQ ². In the measuredQ ² range, the diffractive structure function approximately scales withQ ² at fixed . In an Ingelman-Schlein type model, where commonly used pomeron flux factor normalisations are assumed, it is found that the quarks within the pomeron do not saturate the momentum sum rule.supported by Worldlab, Lausanne, Switzerland     

Variational spiked oscillator

A variational analysis of the spiked harmonic oscillator Hamiltonian –d ²/dr² +r ² +/r ^5/2,>0, is reported. A trial function automatically satisfying both the Dirichlet boundary condition at the origin and the boundary condition at infinity is introduced. The results are excellent for a very large range of values of the coupling parameter, suggesting that the present variational function is appropriate for the treatment of the spiked oscillator in all its regimes (strong, moderate, and weak interactions).