Spectral Degeneracies in the Totally Asymmetric Exclusion Process

We study the spectrum of the Markov matrix of the totally asymmetric exclusion process (TASEP) on a one-dimensional periodic lattice at ARBITRARY filling. Although the system does not possess obvious symmetries except translation invariance, the spectrum presents many multiplets with degeneracies of high order. This behaviour is explained by a hidden symmetry property of the Bethe Ansatz. Combinatorial formulae for the orders of degeneracy and the corresponding number of multiplets are derived and compared with numerical results obtained from exact diagonalisation of small size systems. This unexpected structure of the TASEP spectrum suggests the existence of an underlying large invariance group. Keywords: ASEP, Markov matrix, Bethe Ansatz, Symmetries.Comment: 19 pages, 1 figur

Interacting Random Walkers and Non-Equilibrium Fluctuations

We introduce a model of interacting Random Walk, whose hopping amplitude depends on the number of walkers/particles on the link. The mesoscopic counterpart of such a microscopic dynamics is a diffusing system whose diffusivity depends on the particle density. A non-equilibrium stationary flux can be induced by suitable boundary conditions, and we show indeed that it is mesoscopically described by a Fourier equation with a density dependent diffusivity. A simple mean-field description predicts a critical diffusivity if the hopping amplitude vanishes for a certain walker density. Actually, we evidence that, even if the density equals this pseudo-critical value, the system does not present any criticality but only a dynamical slowing down. This property is confirmed by the fact that, in spite of interaction, the particle distribution at equilibrium is simply described in terms of a product of Poissonians. For mesoscopic systems with a stationary flux, a very effect of interaction among particles consists in the amplification of fluctuations, which is especially relevant close to the pseudo-critical density. This agrees with analogous results obtained for Ising models, clarifying that larger fluctuations are induced by the dynamical slowing down and not by a genuine criticality. The consistency of this amplification effect with altered coloured noise in time series is also proved.Comment: 8 pages, 7 figure

The triangular Ising antiferromagnet in a staggered field

We study the equilibrium properties of the nearest-neighbor Ising antiferromagnet on a triangular lattice in the presence of a staggered field conjugate to one of the degenerate ground states. Using a mapping of the ground states of the model without the staggered field to dimer coverings on the dual lattice, we classify the ground states into sectors specified by the number of ``strings''. We show that the effect of the staggered field is to generate long-range interactions between strings. In the limiting case of the antiferromagnetic coupling constant J becoming infinitely large, we prove the existence of a phase transition in this system and obtain a finite lower bound for the transition temperature. For finite J, we study the equilibrium properties of the system using Monte Carlo simulations with three different dynamics. We find that in all the three cases, equilibration times for low field values increase rapidly with system size at low temperatures. Due to this difficulty in equilibrating sufficiently large systems at low temperatures, our finite-size scaling analysis of the numerical results does not permit a definite conclusion about the existence of a phase transition for finite values of J. A surprising feature in the system is the fact that unlike usual glassy systems, a zero-temperature quench almost always leads to the ground state, while a slow cooling does not.Comment: 12 pages, 18 figures: To appear in Phys. Rev.

Differential Equations for Definition and Evaluation of Feynman Integrals

It is shown that every Feynman integral can be interpreted as Green function of some linear differential operator with constant coefficients. This definition is equivalent to usual one but needs no regularization and application of $R$-operation. It is argued that presented formalism is convenient for practical calculations of Feynman integrals.Comment: pages, LaTEX, MSU-PHYS-HEP-Lu2/9

On the Fluctuation Relation for Nose-Hoover Boundary Thermostated Systems

We discuss the transient and steady state fluctuation relation for a mechanical system in contact with two deterministic thermostats at different temperatures. The system is a modified Lorentz gas in which the fixed scatterers exchange energy with the gas of particles, and the thermostats are modelled by two Nos\'e-Hoover thermostats applied at the boundaries of the system. The transient fluctuation relation, which holds only for a precise choice of the initial ensemble, is verified at all times, as expected. Times longer than the mesoscopic scale, needed for local equilibrium to be settled, are required if a different initial ensemble is considered. This shows how the transient fluctuation relation asymptotically leads to the steady state relation when, as explicitly checked in our systems, the condition found in [D.J. Searles, {\em et al.}, J. Stat. Phys. 128, 1337 (2007)], for the validity of the steady state fluctuation relation, is verified. For the steady state fluctuations of the phase space contraction rate \zL and of the dissipation function \zW, a similar relaxation regime at shorter averaging times is found. The quantity \zW satisfies with good accuracy the fluctuation relation for times larger than the mesoscopic time scale; the quantity \zL appears to begin a monotonic convergence after such times. This is consistent with the fact that \zW and \zL differ by a total time derivative, and that the tails of the probability distribution function of \zL are Gaussian.Comment: Major revision. Fig.10 was added. Version to appear in Journal of Statistical Physic

Testing the Gaussian expansion method in exactly solvable matrix models

The Gaussian expansion has been developed since early 80s as a powerful analytical method, which enables nonperturbative studies of various systems using `perturbative' calculations. Recently the method has been used to suggest that 4d space-time is generated dynamically in a matrix model formulation of superstring theory. Here we clarify the nature of the method by applying it to exactly solvable one-matrix models with various kinds of potential including the ones unbounded from below and of the double-well type. We also formulate a prescription to include a linear term in the Gaussian action in a way consistent with the loop expansion, and test it in some concrete examples. We discuss a case where we obtain two distinct plateaus in the parameter space of the Gaussian action, corresponding to different large-N solutions. This clarifies the situation encountered in the dynamical determination of the space-time dimensionality in the previous works.Comment: 30 pages, 15 figures, LaTeX; added references for section

Activated Random Walkers: Facts, Conjectures and Challenges

We study a particle system with hopping (random walk) dynamics on the integer lattice $\mathbb Z^d$. The particles can exist in two states, active or inactive (sleeping); only the former can hop. The dynamics conserves the number of particles; there is no limit on the number of particles at a given site. Isolated active particles fall asleep at rate $\lambda > 0$, and then remain asleep until joined by another particle at the same site. The state in which all particles are inactive is absorbing. Whether activity continues at long times depends on the relation between the particle density $\zeta$ and the sleeping rate $\lambda$. We discuss the general case, and then, for the one-dimensional totally asymmetric case, study the phase transition between an active phase (for sufficiently large particle densities and/or small $\lambda$) and an absorbing one. We also present arguments regarding the asymptotic mean hopping velocity in the active phase, the rate of fixation in the absorbing phase, and survival of the infinite system at criticality. Using mean-field theory and Monte Carlo simulation, we locate the phase boundary. The phase transition appears to be continuous in both the symmetric and asymmetric versions of the process, but the critical behavior is very different. The former case is characterized by simple integer or rational values for critical exponents ($\beta = 1$, for example), and the phase diagram is in accord with the prediction of mean-field theory. We present evidence that the symmetric version belongs to the universality class of conserved stochastic sandpiles, also known as conserved directed percolation. Simulations also reveal an interesting transient phenomenon of damped oscillations in the activity density

Twenty five years after KLS: A celebration of non-equilibrium statistical mechanics

When Lenz proposed a simple model for phase transitions in magnetism, he couldn't have imagined that the "Ising model" was to become a jewel in field of equilibrium statistical mechanics. Its role spans the spectrum, from a good pedagogical example to a universality class in critical phenomena. A quarter century ago, Katz, Lebowitz and Spohn found a similar treasure. By introducing a seemingly trivial modification to the Ising lattice gas, they took it into the vast realms of non-equilibrium statistical mechanics. An abundant variety of unexpected behavior emerged and caught many of us by surprise. We present a brief review of some of the new insights garnered and some of the outstanding puzzles, as well as speculate on the model's role in the future of non-equilibrium statistical physics.Comment: 3 figures. Proceedings of 100th Statistical Mechanics Meeting, Rutgers, NJ (December, 2008

Cherenkov radiation emitted by ultrafast laser pulses and the generation of coherent polaritons

We report on the generation of coherent phonon polaritons in ZnTe, GaP and LiTaO$_{3}$ using ultrafast optical pulses. These polaritons are coupled modes consisting of mostly far-infrared radiation and a small phonon component, which are excited through nonlinear optical processes involving the Raman and the second-order susceptibilities (difference frequency generation). We probe their associated hybrid vibrational-electric field, in the THz range, by electro-optic sampling methods. The measured field patterns agree very well with calculations for the field due to a distribution of dipoles that follows the shape and moves with the group velocity of the optical pulses. For a tightly focused pulse, the pattern is identical to that of classical Cherenkov radiation by a moving dipole. Results for other shapes and, in particular, for the planar and transient-grating geometries, are accounted for by a convolution of the Cherenkov field due to a point dipole with the function describing the slowly-varying intensity of the pulse. Hence, polariton fields resulting from pulses of arbitrary shape can be described quantitatively in terms of expressions for the Cherenkov radiation emitted by an extended source. Using the Cherenkov approach, we recover the phase-matching conditions that lead to the selection of specific polariton wavevectors in the planar and transient grating geometry as well as the Cherenkov angle itself. The formalism can be easily extended to media exhibiting dispersion in the THz range. Calculations and experimental data for point-like and planar sources reveal significant differences between the so-called superluminal and subluminal cases where the group velocity of the optical pulses is, respectively, above and below the highest phase velocity in the infrared.Comment: 13 pages, 11 figure