Computer simulation of shock waves in the completely asymmetric simple exclusion process

We study the evolution of the completely asymmetric simple exclusion process in one dimension, with particles moving only to the right, for initial configurations corresponding to average density _– ( ₊) left (right) of the origin, _{– +}. The microscopic shock position is identified by introducing a second-class particle. Results indicate that the shock profile is stable, and that the distribution as seen from the shock positionN(t) tends, as time increases, to a limiting distribution, which is locally close to an equilibrium distribution far from the shock. Moreover , withV=1– _–– ₊, as predicted, and the dispersion ofN(t), ²(t), behaves linearly, for not too small values of ₊– _–, i.e., , whereS is equal, up to a scaling factor, to the valueS _WA predicted in the weakly asymmetric case. For ₊= _– we find agreement with the conjecture .Dedicated to the memory of Paola Calderoni.     

Ensemble and trajectory statistics in a nonlinear partial differential equation

We have examined the influence of parametric noise on the solution behavioru(t, x) of a nonlinear initial value() problem arising in cell kinetics. In terms of ensemble statistics, the eventual limiting solution mean and variance are well-characterized functions of the noise statistics, and and depend on . When noise is continuously present along the trajectory, and are independent of the noise statistics and . However, in their evolution toward and , both _u (t, x) and _u ² (t, x) depend on the noise and.     

Generalized two-level quantum dynamics. II. Non-Hamiltonian state evolution

A theorem is derived that enables a systematic enumeration of all the linear superoperators (associated with a two-level quantum system) that generate, via the law of motion = , mappings (0) (t) restricted to the domain of statistical operators. Such dynamical evolutions include the usual Hamiltonian motion as a special case, but they also encompass more general motions, which are noncyclic and feature a destination state (t ) that is in some cases independent of (0).     

Existence of Pair Potential Corresponding to Specified Density and Pair Correlation

Given a potential of pair interaction and a value of activity, one can consider the Gibbs distribution in a finite domain . It is well known that for small values of activity there exist the infinite volume limiting Gibbs distribution and the infinite volume correlation functions. In this paper we consider the converse problem – we show that given ₁ and ₂(x), where ₁ is a constant and ₂(x) is a function on , which are sufficiently small, there exist a pair potential and a value of activity, for which ₁ is the density and ₂(x) is the pair correlation function.Partially supported by NSF Research Grant     

Lyapunov Exponent and Density of States of a One-Dimensional Non-Hermitian Schrödinger Equation

We calculate, using numerical methods, the Lyapunov exponent (E) and the density of states (E) at energy E of a one-dimensional non-Hermitian Schrödinger equation with off-diagonal disorder. For the particular case we consider, both (E) and (E) depend only on the modulus of E. We find a pronounced maximum of (|E|) at energy E=2/ , which seems to be linked to the fixed point structure of an associated random map. We show how the density of states (E) can be expanded in powers of E. We find (|E|)=(1/ ²)+(4/3 ³) |E|²+. This expansion, which seems to be asymptotic, can be carried out to an arbitrarily high order.     

Diffractive production of ρ-mesons and of ρπ-systems by neutrinos and antineutrinos on protons

Evidence is presented for diffractive production of -mesons and of -systems invp and chargedcurrent interactions. In the (anti-)neutrino energy range 10 GeVE _v <60 gev=" the=" cross=" sections=" for=" diffractive="> and diffractive production are found to be (0.64±0.14 (stat.)±0.08 (syst.))% and (0.28±0.08 (stat.)±0.04 (syst.))% of the charged-current cross section. The diffractive signal is consistent with being entirely due to diffractivea ₁ production. However, the data cannot distinguish between diffractivea ₁ and diffractive nonresonant production. The experimental distributions ofW, Q ²,x _Bj andy _Bj for diffractive and events are consistent with model predictions.     

Large Deviation of the Density Profile in the Steady State of the Open Symmetric Simple Exclusion Process

We consider an open one dimensional lattice gas on sites i=1,..., N, with particles jumping independently with rate 1 to neighboring interior empty sites, the simple symmetric exclusion process. The particle fluxes at the left and right boundaries, corresponding to exchanges with reservoirs at different chemical potentials, create a stationary nonequilibrium state (SNS) with a steady flux of particles through the system. The mean density profile in this state, which is linear, describes the typical behavior of a macroscopic system, i.e., this profile occurs with probability 1 when N. The probability of microscopic configurations corresponding to some other profile (x), x=i/N, has the asymptotic form exp[–N ({})]; is the large deviation functional. In contrast to equilibrium systems, for which _eq({}) is just the integral of the appropriately normalized local free energy density, the we find here for the nonequilibrium system is a nonlocal function of . This gives rise to the long range correlations in the SNS predicted by fluctuating hydrodynamics and suggests similar non-local behavior of in general SNS, where the long range correlations have been observed experimentally.     

Power-law falloff in the kosterlitz-Thouless phase of a two-dimensional lattice Coulomb gas

We give a rigorous proof of power-law falloff in the Kosterlitz-Thouless phase of a two-dimensional Coulomb gas in the sense that there exists a critical inverse temperaturegb and a constant >0 such that for all> and all external charges R we have , whereG (x) is the two-point external charges correlation function,=dist(, Z), and for 0$$ " align="middle" border="0"> . In the case of a hard-core or standard Coulomb gas with activityz, we may choose=(z) such that(z)24 asz0.     

Metastable fluid flow described via a discrete-velocity coagulation-fragmentation model

A discrete-velocity Boltzmann model is introduced. It is based on two principles: (i) clusters of particles move in ³ with seven fixed momenta; (ii) clusters may gain or lose particles according to the rules of Becker-Döring cluster equations. The model provides a kinetic representation of evaporation and condensation. The model is used to obtain macroscopic fluid equations which are valid into the metastable fluid regime, , where is any positive number, is the inelastic Knudsen number, and _s is the saturation density.     

Bounded Fluctuations and Translation Symmetry Breaking: A Solvable Model

The variance of the particle number (equivalently the total charge) in a domain of length of a one-component plasma (OCP) on a cylinder of circumference W at the reciprocal temperature =2, is shown to remain bounded as . This exactly solvable system with average density has a density profile which is periodic with period (W)^–1 along the axis of the infinitely long cylinder. This illustrates the connection between bounded variance and periodicity in (quasi) one-dimensional systems.⁽¹⁾ When W the system approaches the two-dimensional OCP and the variance in a domain grows like its perimeter ||. In this limit, the system is translation invariant with rapid decay of correlations.     

Electron-electron interaction and magnon scattering in the d-electron dominated resistivity of Ni80−xFe x Si8B12 glasses

Our investigation of the electrical resistivity (T) betweenT=0.1 K and 295 K revealed for all Fe-concentrations, 2.4x16, a minimum of betweenT _min=7 K and 15 K. Above the minimum, scattering from structural disorder prevails in the paramagnetic alloy (x=2.4), while in the ferromagnetic glasses (x5.6) coherent magnon scattering appears to dominate. With one exception, (TT _min) increases as down to the lowest temperatures for all investigated Fe-contents (x10.4) in zero magnetic field as well as in applied 6 T. The value of the slopes, 6.3(6) ( cm)^–1 K^–1/2, agrees with the universal scaling behaviour reported recently by Cochrane and Strom-Olsen for systems with strong electron-electron interaction. A large, linear rise of (T<T _min,B=0) observed in the most diluted sample (x=2.4) is tentatively attributed to quantum localization in the non-reentrant spin glass phase. The magnitude of the temperature effects on (T) above as well as belowT _min indicates significantd-band contributions to the conductivity.     

The decays of theρ′1 andω′1 mesons

New measurements of the reactione ⁺ e ^– confirm previous evidence for the decay₁ and provide a strong constraint on the mass of the₁. The results from a detailed analysis of these data are compared with those from a corresponding analysis of new data on the reactione ⁺ e ^–. It is shown that in accord with expectation. This identification of a strong₁ decay supports a previous claim that the₁ and theC(1480), observed in the channel, cannot be different manifestations of the same particle.     

Surface (110) of elementary semiconductors with diamond-type structure during stratified growth from an atomic beam in a vacuum

The stratified growth of an ideal crystal with diamond structure from an atomic beam in a vacuum is investigated by using the kinetic method of the theory of crystal growth. Equations are obtained to determine the mean concentrations _n of molecules consisting of n atoms on the (110) surface of the crystal. In the stationary case, an approximate analytic solution is obtained with a nonlinear dependence of ₁, ₂, and on the atomic beam density.Translated from Izvestiya Vysshikh Uchebnykh Zavedenii, Fizika, No. 12, pp. 92–96, Decmeber, 1981.     

Model for γγ→VV reactions near threshold

We propose a Veneziano-like model for VV reactions near threshold involving exchanges of 0^– trajectories int- andu-channels. The model is applied to , , , and processes. The model in its simplest form cannot account for the disparities in the cross sections for andK ^*+ K ^*– and for ⁰⁰ and ^+–. Further refinements of the model are suggested.     

Radius,perimeter, and density profile for percolation clusters and lattice animals

Monte Carlo simulation and series expansion shows the radius of gyration of large clusters withs sites each to vary ass with0.56 in two and0.47 in three dimensions at the percolation threshold, and with(d=2)0.65 and(d=3)0.53 for random lattice animals (zero concentration). Clusters up tos=100 were used. The perimeter of random animals approaches 2.8s for larges on the simple cubic lattice. Monte Carlo simulation of the Eden process (growing animals) up tos=5,000 indicates a systematic variation of about ±0.05 for the effective exponent=(s) and thus suggests that the true asymptotic exponents may be compatible with the predictions of hyper-scaling.     

Trapping,percolation, and anomalous diffusion of particles in a two-dimensional random field

We analyze from first principles the advection of particles in a velocity field with HamiltonianH(x, y)=¯ V ₁ y–¯ V ₂ x+W ₁ (y)-W ₂ (x), whereW _i , i=1, 2, are random functions with stationary, independent increments. In the absence of molecular diffusion, the particle dynamics are very sensitive to the streamline topology, which depends on the mean-to-fluctuations ratio=max(|¯V₁¦/; ¦¯V₂|/), with =|W ₁ |^21/2=rms fluctuations. Remarkably, the model is exactly solvable for 1 and well suited for Monte Carlo simulations for all , providing a nice setting for studying seminumerically the influence of streamline topology on large-scale transport. First, we consider the statistics of streamlines for=0, deriving power laws for p_nc(L) and (L), which are, respectively, the escape probability and the length of escaping trajectories for a box of sizeL, L » 1. We also obtain a characterization of the statistical topography of the HamiltonianH. Second, we study the large-scale transport of advected particles with > 0. For 0 < < 1, a fraction of particles is trapped in closed field lines and another fraction undergoes unbounded motions; while for 1 all particles evolve in open streamlines. The fluctuations of the free particle positions about their mean is studied in terms of the normalized variablest ^– v/2[x(t)–x(t)] andt ^–v/2 [y(t)-(t)]. The large-scale motions are shown to be either Fickian (=1), or superdiffusive (=3/2) with a non-Gaussian coarse-grained probability, according to the direction of the mean velocity relative to the underlying lattice. These results are obtained analytically for 1 and extended to the regime 0<<1 by Monte Carlo simulations. Moreover, we show that the effective diffusivity blows up for resonant values of ) for which stagnation regions in the flow exist. We compare the results with existing predictions on the topology of streamlines based on percolation theory, as well as with mean-field calculations of effective diffusivities. The simulations are carried out with a CM 200 massively parallel computer with 8192 SIMD processors.     

Dynamical Localization II with an Application to the Almost Mathieu Operator

Several recent works have established dynamical localization for Schrödinger operators, starting from control on the localization length of their eigenfunctions, in terms of their centers of localization. We provide an alternative way to obtain dynamical localization, without resorting to such a strong condition on the exponential decay of the eigenfunctions. Furthermore, we illustrate our purpose with the almost Mathieu operator, H _{, ,} =–+ cos(2(+x)), 15 and with good Diophantine properties. More precisely, for almost all , for all q>0, and for all functions ²( ) of compact support, we show that The proof applies equally well to discrete and continuous random Hamiltonians. In all cases, it uses as input a repulsion principle of singular boxes, supplied in the random case by the multi-scale analysis.     

Examples of Twisted Cyclic Cocycles from Covariant Differential Calculi

For two covariant differential *-calculi, the twisted cyclic cocycle associated with the volume form is represented in terms of commutators for some self-adjoint operator and some *-representation of the underlying *-algebra.     

Absence of symmetry breaking for systems of rotors with random interactions

We prove that Gibbs states for the Hamiltonian , with thes _x varying on theN-dimensional unit sphere, obtained with nonrandom boundary conditions (in a suitable sense), are almost surely rotationally invariant if withJ _xy i.i.d. bounded random variables with zero average, 1 in one dimension, and 2 in two dimensions.     

Microscopic selection principle for a diffusion-reaction equation

We consider a model of stochastically interacting particles on , where each site is assumed to be empty or occupied by at most one particle. Particles jump to each empty neighboring site with rate/2 and also create new particles with rate 1/2 at these sites. We show that as seen from the rightmost particle, this process has precisely one invariant distribution. The average velocity of this particle V() then satisfies ^–1/2V() as. This limit corresponds to that of the macroscopic density obtained by rescaling lengths by a factor ^1/2 and letting. This density solves the reaction-diffusion equation , and under Heaviside initial data converges to a traveling wave moving at the same rate .