Convergence to equilibrium distribution. The Klein-Gordon equation coupled to a particle

The Hamiltonian system formed by a Klein-Gordon vector field and a particle in ℝ³ is considered. The initial data of the system are given by a random function, with finite mean energy density, which also satisfies a Rosenblatt- or Ibragimov-type mixing condition. Moreover, initial correlation functions are assumed to be translation invariant. The distribution μ _t of the solution at time t ∈ ℝ is studied. The main result is the convergence of μ _t to a Gaussian measure as t → ∞, where μ_∞ is translation invariant.     

On Convergence to Equilibrium Distribution,I.¶The Klein–Gordon Equation with Mixing

Consider the Klein–Gordon equation (KGE) in ℝ ⁿ , n≥ 2, with constant or variable coefficients. We study the distribution μ _t of the random solution at time t∈ℝ. We assume that the initial probability measure μ₀ has zero mean, a translation-invariant covariance, and a finite mean energy density. We also assume that μ₀ satisfies a Rosenblatt- or Ibragimov–Linnik-type mixing condition. The main result is the convergence of μ _t to a Gaussian probability measure as t→∞ which gives a Central Limit Theorem for the KGE. The proof for the case of constant coefficients is based on an analysis of long time asymptotics of the solution in the Fourier representation and Bernstein's “room-corridor” argument. The case of variable coefficients is treated by using an “averaged” version ofthe scattering theory for infinite energy solutions, based on Vainberg's results on local energy decay. Received: 4 January 2001 / Accepted: 2 July 2001     

Distribution of a Particle’s Position in the ASEP with the Alternating Initial Condition

In this paper we give the distribution of the position of a particle in the asymmetric simple exclusion process (ASEP) with the alternating initial condition. That is, we find ℙ(X _m (t)≤x) where X _m (t) is the position of the particle at time t which was at m=2k−1, k∈ℤ at t=0. As in the ASEP with step initial condition, there arises a new combinatorial identity for the alternating initial condition, and this identity relates the integrand of the integral formula for ℙ(X _m (t)≤x) to a determinantal form together with an extra product.     

Renormalization group analysis of polymer cyclization with non-equilibrium initial conditions

We develop a renormalization group approach for cyclizing polymers for the case when chain ends are initially close together (ring initial conditions). We analyze the behavior at times much shorter than the longest polymer relaxation time. In agreement with our previous work (Europhys. Lett. 73, 621 (2006)) we find that the leading time dependence of the reaction rate k(t) for ring initial conditions and equilibrium initial conditions are related, namely k _ring(t) ∝ t ^-δ and k _eq(t) ∝ t ^1-δ for times less than the longest polymer relaxation time. Here δ is an effective exponent which approaches δ = 5/4 for very long Rouse chains. Our present analysis also suggests a “sub-leading” term proportional to (ln t)/t which should be particularly significant for smaller values of the renormalized reaction rate and early times. For Zimm dynamics, our RG analysis indicates that the leading time dependence for the reaction rate is k(t) ∼ 1/t for very long chains. The leading term is again consistent with the expected relation between ring and equilibrium initial conditions. We also find a logarithmic correction term which we “exponentiate” to a logarithmic form with a Landau pole. The presence of the logarithm is particularly important for smaller chains and, in the Zimm case, large values of the reaction rate.     

Non-linear scission/recombination kinetics of living polymerization

Living polymers are formed by reversible association of primary units (unimers). Generally the chain statistical weight involves a factor σ < 1 suppressing short chains in comparison with free unimers. Living polymerization is a sharp thermodynamic transition for σ ≪ 1 which is typically the case. We show that this sharpness has an important effect on the kinetics of living polymerization (one-dimensional association). The kinetic model involves i) the unimer activation step (a transition to an assembly-competent state); ii) the scission/recombination processes providing growth of polymer chains and relaxation of their length distribution. Analyzing the polymerization with no chains but unimers at t = 0 , with initial concentration of unimers M ≳ M ^* (M^* is the critical polymerization concentration), we determine the time evolution of the chain length distribution and find that: 1) for M ^* ≪ M ≪ M ^*/σ the kinetics is characterized by 5 distinct time stages demarcated by 4 characteristic times t₁, t₂, t₃ and t^*; 2) there are transient regimes (t ₁ ≲ t ≲ t ₃) when the molecular-weight distribution is strongly non-exponential; 3) the chain scissions are negligible at times shorter than t₂. The chain growth is auto-accelerated for t ₁ ≲ t ≲ t ₂ : the cut-off chain length (= polymerization degree 〈n〉_w N ₁ ∝ t ² in this regime. 4) For t ₂ < t < t ₃ the length distribution is characterized by essentially 2 non-linear modes; the shorter cut-off length N₁ is decreasing with time in this regime, while the length scale N₂ of the second mode is increasing. (5) The terminal relaxation time of the polymer length distribution, t^*, shows a sharp maximum in the vicinity of M^*; the effective exponent is as high as ∼ σ^-1/3 just above M^*.     

Reacting polymers with highly correlated initial conditions

We propose and theoretically study an experiment designed to measure short-time polymer reaction kinetics in melts or dilute solutions. The photolysis of groups centrally located along chain backbones, one group per chain, creates pairs of spatially highly correlated macroradicals. We calculate time-dependent rate coefficients κ(t) governing their first-order recombination kinetics, which are novel on account of the far-from-equilibrium initial conditions. In dilute solutions (good solvents) reaction kinetics are intrinsically weak, despite the highly reactive radical groups involved. This leads to a generalised mean-field kinetics in which the rate of radical density decay - ∼S(t), where S(t) ∼t ^{- (1 + g/3)} is the equilibrium return probability for 2 reactive groups, given initial contact. Here g≈ 0.27 is the correlation hole exponent for self-avoiding chain ends. For times beyond the longest coil relaxation time τ, - ∼S(t) remains true, but center of gravity coil diffusion takes over with rms displacement of reactive groups x(t) ∼t ^1/2 and S(t) ∼ 1/x ³(t). At the shortest times ( t 10^-6s), recombination is inhibited due to spin selection rules and we find ∼tS(t). In melts, kinetics are intrinsically diffusion-controlled, leading to entirely different rate laws. During the regime limited by spin selection rules, the density of radicals decays linearly, n(0) - n(t) ∼t. At longer times the standard result - ∼d ³(t)/d (for randomly distributed ends) is replaced by ∼d2x ³(t)/d ² for these correlated initial conditions. The long-time behavior, t > τ, has the same scaling form in time as for dilute solutions. Received 18 May 2000     

Three-dimensional array structures associated with Richtmyer-Meshkov and Rayleigh-Taylor instability

A boundary separating adjacent gas or liquid media is frequently unstable. Richtmyer-Meshkov and Rayleigh-Taylor instability cause the growth of intricate structures on such boundaries. All the lattice symmetries [rectangular (pmm2), square (p4mm), hexagonal (p6mm), and triangular (p3m1) lattices] which are of interest in connection with the instability of the surface of a fluid are studied for the first time. They are obtained from initial disturbances consisting of one (planar case, two-dimensional flow), two (rectangular cells), or three (hexagons and triangles) harmonic waves. It is shown that the dynamic system undergoes a transition during development from an initial, weakly disturbed state to a limiting or asymptotic stationary state (stationary point). The stability of these points (stationary states) is investigated. It is shown that the stationary states are stable toward large-scale disturbances both in the case of Richtmyer-Meshkov instability and in the case of Rayleigh-Taylor instability. It is discovered that the symmetry increases as the system evolves in certain cases. In one example the initial Richtmyer-Meshkov or Rayleigh-Taylor disturbance is a sum of two waves perpendicular to one another with equal wave numbers, but unequal amplitudes: a ₁(t=0)≠a ₂(t=0). Then, during evolution, the flow has p2 symmetry (rotation relative to the vertical axis by 180°), which goes over to p4 symmetry (rotation by 90°) at t→∞, since the amplitudes equalize in the stationary state: a ₁(t=∞)=a ₂(t=∞). It is shown that the hexagonal and triangular arrays are complementary. Upon time inversion (t→−t), “rephasing” occurs, and the bubbles of a hexagonal array transform into jets of a triangular array and vice versa. Zh. éksp. Teor. Fiz. 116, 908–939 (September 1999)     

Possible Loss and Recovery of Gibbsianness¶During the Stochastic Evolution of Gibbs Measures

We consider Ising-spin systems starting from an initial Gibbs measure ν and evolving under a spin-flip dynamics towards a reversible Gibbs measure μ≠ν. Both ν and μ are assumed to have a translation-invariant finite-range interaction. We study the Gibbsian character of the measure νS(t) at time t and show the following: (1) For all ν and μ, νS(t) is Gibbs for small t. (2) If both ν and μ have a high or infinite temperature, then νS(t) is Gibbs for all t > 0. (3) If ν has a low non-zero temperature and a zero magnetic field and μ has a high or infinite temperature, then νS(t) is Gibbs for small t and non-Gibbs for large t. (4) If ν has a low non-zero temperature and a non-zero magnetic field and μ has a high or infinite temperature, then νS(t) is Gibbs for small t, non-Gibbs for intermediate t, and Gibbs for large t. The regime where μ has a low or zero temperature and t is not small remains open. This regime presumably allows for many different scenarios. Received: 26 April 2001 / Accepted: 10 October 2001     

On ergodic properties for harmonic crystals

We consider the dynamics of a harmonic crystal in n dimensions with d components, where d and n are arbitrary, d, n ⩾ 1. The initial data are given by a random function with finite mean energy density which also satisfies a Rosenblatt-or Ibragimov-type mixing condition. The random function is close to diverse space-homogeneous processes as x _n → ±∞, with the distributions μ±. We prove that the phase flow is mixing with respect to the limit measure of statistical solutions. Partially supported by RFBR under grant no. 06-01-00096.     

Formulas for ASEP with Two-Sided Bernoulli Initial Condition

For the asymmetric simple exclusion process η _t on the integer lattice with two-sided Bernoulli initial condition, we derive exact formulas for the following quantities: (1) ℙ(η _t (x)=1), the probability that site x is occupied at time t; (2) the correlation function ℙ(η _t (x)=1,η ₀(0)=1); (3) the distribution function for Q _t , the total flux across 0 at time t, and its exponential generating function.     

Smoothening of a random surface as it results from the edwards-wilkinson kinetic equation

The time evolution of a random surfacez=h(r, t) (r=x, y) formed by a deposition process of the Edwards-Wilkinson type is discussed. The discussion is based on the author’s former derivation of the autocorrelation functionA _h(|r − r′|,t, t′)=〈h(r,t)h(r′,t′)〉 of the height functionh(r,t) under the assumption of a stochastic initial condition [V. Bezák: Acta Physica Univ. Comenianae39 (1998) 135]. Under the assumption of a steady (non-zero) deposition rate, the varianceσ _h ² (t)=〈[h(r,t)]²〉 increases logarithmically in time whilst the correlation lengthl _h(t) of the height functionh(r,t) increases as ∼t ^1/2. Therefore, the ratioσ _h(t)/l _h (t) tends to zero and the surfacez=h(r,t) does always tend towards a smoothened appearance. This work has been supported by the Slovak Grant Agency VEGA under contract No. 1/4319/97.     

Dynamical Localization of Quantum Walks in Random Environments

We study shock statistics in the scalar conservation law ∂ _t u+∂ _x f(u)=0, x∈ℝ, t>0, with a convex flux f and spatially random initial data. We show that the Markov property (in x) is preserved for a large class of random initial data (Markov processes with downward jumps and derivatives of Lévy processes with downward jumps). The kinetics of shock clustering is then described completely by an evolution equation for the generator of the Markov process u(x,t), x∈ℝ. We present four distinct derivations for this evolution equation, and show that it takes the form of a Lax pair. The Lax equation admits a spectral parameter as in Manakov (Funct. Anal. Appl. 10:328–329, 1976), and has remarkable exact solutions for Burgers equation (f(u)=u ²/2). This suggests the kinetic equations of shock clustering are completely integrable.     

Temporal correlations and persistence in the kinetic Ising model: the role of temperature

We study the statistical properties of the sum S _t = dt'σ _t', that is the difference of time spent positive or negative by the spin σ _t, located at a given site of a D-dimensional Ising model evolving under Glauber dynamics from a random initial configuration. We investigate the distribution of S_t and the first-passage statistics (persistence) of this quantity. We discuss successively the three regimes of high temperature ( T > T _c), criticality ( T = T _c), and low temperature ( T < T _c). We discuss in particular the question of the temperature dependence of the persistence exponent , as well as that of the spectrum of exponents (x), in the low temperature phase. The probability that the temporal mean S _t/t was always larger than the equilibrium magnetization is found to decay as t ^{- - ?}. This yields a numerical determination of the persistence exponent in the whole low temperature phase, in two dimensions, and above the roughening transition, in the low-temperature phase of the three-dimensional Ising model. Received 4 December 2000     

Effect of time reversal in the magnetization of an atom by a resonant light pulse

It is shown that the even dependence of the light-induced magnetic moment on the detuning ω-ω _ba from resonance in the case of a circularly polarized pulse and an isotropic initial state of the atom and the odd dependence on ω-ω _ba in the case of a linearly polarized pulse and an anisotropic initial state in the form of alignment of the atom are consequences of the symmetry under time reversal t→−t and of the initial conditions at time t=0. In a number of cases, this fundamental law makes it possible to determine the vector properties of a light-induced magnetic moment and its dependence on the time t and ω-ω _ba without solving the equation for the density matrix in detail and without calculating the sum over the projections of the angular momenta in the formula for the magnetization of an atom by light. Pis’ma Zh. éksp. Teor. Fiz. 65, No. 3, 231–236 (10 February 1997)     

Maximum of a Fractional Brownian Motion: Probabilities of Small Values

Let b _γ (t), b _γ(0)= 0 be a fractional Brownian motion, i.e., a Gaussian process with the structure function , 0 < γ < 2. We study the logarithmic asymptotics of P _T = P{b _γ (t) < 1,□t∈TΔ} as T→∞, where Δ is either the interval (0,1) or a bounded region that contains a vicinity of 0 for the case of multidimensional time. It is shown that ln P _T = - D ln T(1 + o(1)), where D is the dimension of zeroes of b _γ (t) in the former case and the dimension of time in the latter. Received: 28 September 1998 / Accepted: 19 February 1999     

Current Fluctuations of the One Dimensional Symmetric Simple Exclusion Process with Step Initial Condition

For the symmetric simple exclusion process on an infinite line, we calculate exactly the fluctuations of the integrated current Q _t during time t through the origin when, in the initial condition, the sites are occupied with density ρ _a on the negative axis and with density ρ _b on the positive axis. All the cumulants of Q _t grow like . In the range where , the decay exp [−Q _t ³/t] of the distribution of Q _t is non-Gaussian. Our results are obtained using the Bethe ansatz and several identities derived recently by Tracy and Widom for exclusion processes on the infinite line. We acknowledge the support of the French Ministry of Education through the ANR BLAN07-2184264 grant.     

The corrections to scaling within Mazenko’s theory in the limit of low and high dimensions

We consider corrections to scaling within an approximate theory developed by Mazenko for nonconserved order parameter in the limit of low (d → 1) and high (d → ∞) dimensions. The corrections to scaling considered here follows from the departures of the initial condition from the scaling morphology. Including corrections to scaling, the equal time correlation function has the form: C(r, t) = f ₀(r/L)+L ^−ω f ₁(r/L)+…, where L is a characteristic length scale (i.e. domain size). The correction-to-scaling exponent ω and the correction-to-scaling functions f ₁(x) are calculated for both low and high dimensions. In both dimensions the value of ω is found to be ω = 4 similar to 1D Glauber model and OJK theory (the theory developed by Ohta, Jasnow and Kawasaki).     

Discrete-time quantum walks on one-dimensional lattices

In this paper, we study discrete-time quantum walks on one-dimensional lattices. We find that the coherent dynamics depends on the initial states and coin parameters. For infinite size of lattices, we derive an explicit expression for the return probability, which shows scaling behavior P(0, t) ~ t ^-1 and does not depends on the initial states of the walk. In the long-time limit, the probability distribution shows various patterns, depending on the initial states, coin parameters and the lattice size. The time-averaged probability mixes to the limiting probability distribution in linear time, i.e., the mixing time M _ε is a linear function of N (size of the lattices) for large values of thresholds ϵ. Finally, we introduce another kind of quantum walk on infinite or even-numbered size of lattices, and show that by the method of mathematical induction, the walk is equivalent to the traditional quantum walk with symmetrical initial state and coin parameter.     

Propagation of electromagnetic waves generated by moving sources in dispersive media

The propagation of electromagnetic waves issued by modulated moving sources of the form j( t,x ) = a( t )e ^{- iw₀ t} [(x)\dot]₀ ( t )d( x - x₀ ( t ) )j\left( {t,x} \right) = a\left( t \right)e^{ - i\omega _0 t} \dot x_0 \left( t \right)\delta \left( {x - x_0 \left( t \right)} \right) is considered, where j(t, x) stands for the current density vector, x = (x ₁, x ₂, x ₃) ∈ ℝ³ for the space variables, t ∈ ℝ for time, t → x ₀(t) ∈ ℝ³ for the vector function defining the motion of the source, ω ₀ for the eigenfrequency of the source, a(t) for a narrow-band amplitude, and δ for the standard δ function. Suppose that the media under consideration are dispersive. This means that the electric and magnetic permittivity ɛ(ω), μ(ω) depends on the frequency ω. We obtain a representation of electromagnetic fields in the form of time-frequency oscillating integrals whose phase contains a large parameter λ > 0 characterizing the slowness of the change of the amplitude a(t) and the velocity [(x)\dot]₀ ( t )\dot x_0 \left( t \right) and a large distance between positions of the source and the receiver. Applying the two-dimensional stationary phase method to the integrals, we obtain explicit formulas for the electromagnetic field and for the Doppler effects. As an application of our approach, we consider the propagation of electromagnetic waves produced by moving source in a cold nonmagnetized plasma and the Cherenkov radiation in dispersive media.