Kinetic lattice models of disorder

We study a class of stochastic Ising (or interacting particle) systems that exhibit a spatial distribution of impurities that change with time. It may model, for instance, steady nonequilibrium conditions of the kind that may be induced by diffusion in some disordered materials. Different assumptions for the degree of coupling between the spin and the impurity configurations are considered. Two interesting well-defined limits for impurities that behave autonomously are (i) the standard (i.e., quenched) bond-diluted, random-field, random-exchange, and spin-glass Ising models, and (ii) kinetic variations of these standard cases in which conflicting kinetics simulate fast and random diffusion of impurities. A generalization of the Mattis model with disorder that describes a crossover from the equilibrium case (i) to the nonequilibrium case (ii) and the microscopic structure of a generalized heat bath are explicitly worked out as specific realizations of our class of models. We sketch a simple classification of transition rates for the time evolution of the spin configuration based on the critical behavior that is exhibited by the models in case (ii). The latter are shown to have an exact solution for any lattice dimension for some special choice of rates.     

Decay to equilibrium in random spin systems on a lattice. II

We study the continuous spin systems on ad3-dimensional lattice with random ferromagnetic interactions of finite range. We show that, if the temperature is sufficiently high and the probability of interaction to be large is small enough, the almost sure decay to equilibrium has a subexponential upper bound.     

Generalized voter models

We examine a one-dimensional class of interacting particle systems which generalize some voter models. This class includes a particular case in the class of models of catalytic surfaces introduced by Swindle and Grannan. We show that this class has the clustering property of ordinary finite-range voter models, at least when one is concerned with translation-invariant measures on the state space.Research supported by P.M.S. 9157461.     

Correlation functions for simple fluids in a finite system under nonequilibrium constraints

The Landau-Lifschitz fluctuating hydrodynamics formalism is applied to study the statistical properties of simple fluids in a finite system under nonequilibrium constraints. The boundary conditions are explicitly taken into account so that the results can be compared with particle simulations. Two scenarios are investigated: a fluid subjected to a constant shear and a fluid subjected to a constant temperature gradient. By considering a fluid with vanishing thermal expansivity, exact results are obtained for the static and dynamic correlation functions.     

Long-time integration methods for mesoscopic models of pattern-forming systems

Spectral methods for simulation of a mesoscopic diffusion model of surface pattern formation are evaluated for long simulation times. Backwards-differencing time-integration, coupled with an underlying Newton–Krylov nonlinear solver (SUNDIALS-CVODE), is found to substantially accelerate simulations, without the typical requirement of preconditioning. Quasi-equilibrium simulations of patterned phases predicted by the model are shown to agree well with linear stability analysis. Simulation results of the effect of repulsive particle–particle interactions on pattern relaxation time and short/long-range order are discussed.     

Time-dependent density functional theory for quantum transport

The rapid miniaturization of electronic devices motivates research interests in quantum transport. Recently time-dependent quantum transport has become an important research topic. Here we review recent progresses in the development of time-dependent density-functional theory for quantum transport including the theoretical foundation and numerical algorithms. In particular, the reducedsingle electron density matrix based hierarchical equation of motion, which can be derived from Liouville–von Neumann equation, is reviewed in details. The numerical implementation is discussed and simulation results of realistic devices will be given.     

Field strength formulation,lattice Bianchi identities and perturbation theory for non-Abelian models

We develop an analytical approach for studying lattice gauge theories within the plaquette representation where the plaquette matrices play the role of the fundamental degrees of freedom. We start from the original Batrouni formulation and show how it can be modified in such a way that each non-Abelian Bianchi identity contains only two connectors instead of four. In addition, we include dynamical fermions in the plaquette formulation. Using this representation we construct the low-temperature perturbative expansion for U(1)$U(1)$ and SU(N)$\mathit{SU}(N)$ models and discuss its uniformity in the volume. The final aim of this study is to give a mathematical background for working with non-Abelian models in the plaquette formulation.     

Numerical calculation of domains of analyticity for perturbation theories in the presence of small divisors

We study numerically the complex domains of validity for KAM theory in generalized standard mappings. We compare methods based on Padé approximants and methods based on the study of periodic orbits.     

On the corrections to scaling in three-dimensional Ising models

The leading correction-to-scaling amplitudes for the spin-1/2, nearest-neighbor sc, bcc, and fee Ising models are considered with the particular aim of determining their signs. On the basis of previous two-variable series analyses by Chen, Fisher, and Nickel and renormalization group=4–d expansions, it is concluded that the correction amplitudes for the susceptibility, correlation length, specific heat, and spontaneous magnetization arenegative for all three lattices. Thus, for example, the effective exponent _eff(T) asymptotically approaches the true susceptibility exponent fromabove. Other earlier and more recent high-temperature series and field-theoretic analyses are seen to be consistent with this result. However, the usual nonasymptotic, perturbative field-theoretic approaches are essentially committed to positive correction amplitudes. The question of the signs therefore relates directly to the applicability of these non-asymptotic field-theoretic calculations to three-dimensional Ising models as well as to different experimental systems.     

High-temperature series for scalar-field Lattice models: Generation and analysis

An implementation of the free-embedding scheme for high-temperature series generation on the body-centered cubic family of lattices in arbitrary dimensiond is, described. Series to order 21 in inverse temperature are tabulated for several scalar field models, both for the magnetic susceptibility and for the second moment of the spin correlation function. The critical behavior of a family of 3-dimensional double Gaussian models, which interpolate continuously between the spin-1/2 Ising model and the Gaussian model, is analyzed in detail away from the Gaussian model limit using confluent inhomogeneous secondorder differential approximants. With our best estimate of the correction-to-scaling exponent,=0.52±0.03, the leading exponents for the susceptibility and correlation length for this family are consistent with universality and are given by=1.237±0.002 and =0.630±0.0015, respectively, and=2–/=0.0359±0.0007.     

The continuous-spin Ising model,g0∶φ4∶d field theory,and the renormalization group

We have used the method of high-temperature series expansions to investigate the critical point properties of a continuous-spin Ising model and g₀⁴_d Euclidean field theory. We have computed through tenth order the hightemperature series expansions for the magnetization, susceptibility, second derivative of the susceptibility, and the second moment of the spin-spin correlation function on eight different lattices. Our analysis of these series is made using integral and Padé approximants. In three dimensions we find that hyperscaling fails for sufficiently Ising-like systems; the strong coupling limit of g₀⁴₃ depends on how the ultraviolet cutoff is removed. The level contours of the renormalized coupling constant for this model in theg ₀, correlation-length plane exhibit a saddle point. If the ultraviolet cutoff is removed beforeg ₀ , the usual field theory results and the renormalization-group fixed point with hyperscaling is obtained. If the order of these limits is reversed, the Ising model limit where hyperscaling fails and the field theory is trivial is obtained. In four dimensions, we find that hyperscaling fails completely; g₀⁴₄ is trivial for all g₀ when the ultraviolet cutoff is removed.Work supported in part by the U.S. Department of Energy.     

Transient bimodality in interacting particle systems

We consider a system of spins which have values ±1 and evolve according to a jump Markov process whose generator is the sum of two generators, one describing a spin-flipGlauber process, the other aKawasaki (stirring) evolution. It was proven elsewhere that if the Kawasaki dynamics is speeded up by a factor ^–2, then, in the limit 0 (continuum limit), propagation of chaos holds and the local magnetization solves a reaction-diffusion equation. We choose the parameters of the Glauber interaction so that the potential of the reaction term in the reaction-diffusion equation is a double-well potential with quartic maximum at the origin. We assume further that for each the system is in a finite interval ofZ with ^–1 sites and periodic boundary conditions. We specify the initial measure as the product measure with 0 spin average, thus obtaining, in the continuum limit, a constant magnetic profile equal to 0, which is a stationary unstable solution to the reaction-diffusion equation. We prove that at times of the order ^–1/2 propagation of chaos does not hold any more and, in the limit as 0, the state becomes a nontrivial superposition of Bernoulli measures with parameters corresponding to the minima of the reaction potential. The coefficients of such a superposition depend on time (on the scale ^–1/2) and at large times (on this scale) the coefficient of the term corresponding to the initial magnetization vanishes (transient bimodality). This differs from what was observed by De Masi, Presutti, and Vares, who considered a reaction potential with quadratic maximum and no bimodal effect was seen, as predicted by Broggi, Lugiato, and Colombo.     

Low-temperature phase diagrams of quantum lattice systems. I. Stability for quantum perturbations of classical systems with finitely-many ground states

Starting from classical lattice systems ind2 dimensions with a regular zerotemperature phase diagram, involving a finite number of periodic ground states, we prove that adding a small quantum perturbation and/or increasing the temperature produce only smooth deformations of their phase diagrams. The quantum perturbations can involve bosons or fermions and can be of infinite range but decaying exponentially fast with the size of the bonds. For fermions, the interactions must be given by monomials of even degree in creation and annihilation operators. Our methods can be applied to some anyonic systems as well. Our analysis is based on an extension of Pirogov-Sinai theory to contour expansions ind+1 dimensions obtained by iteration of the Duhamel formula.     

Nonequilibrium lattice models: Series analysis of steady states

A perturbation theory for steady states of interacting particle systems is developed and applied to several lattice models with nonequilibrium critical points near an absorbing state. The expansion is expressed directly in terms of the kinetic parameter (creation rate), rather than in powers of the interaction. An algorithm for generating series expansions for local properties is described. Order parameter series (16 terms) and precise estimates of critical properties are presented for the one-dimensional contact process and several related models.     

Semiclassical theory for many-body fermionic systems

We present a treatment of many-body fermionic systems that facilitates an expression of well-known quantities in a series expansion inħ. The ensuing semiclassical result contains, to a leading order of the response function, the classical time correlation function of the observable followed by the Weyl-Wigner series; on top of these terms are the periodic-orbit correction terms. The treatment given here starts from linear response assumption of the many-body theory and in its connection with semiclassical theory, it assumes that the one-body quantal system has a classically chaotic dynamics. Applications of the framework are also discussed.     

Natural boundaries for area-preserving twist maps

We consider KAM invariant curves for generalizations of the standard map of the form (x, y)=(x+y, y+f(x)), wheref(x) is an odd trigonometric polynomial. We study numerically their analytic properties by a Padé approximant method applied to the function which conjugates the dynamics to a rotation +. In the complex plane, natural boundaries of different shapes are found. In the complex plane the analyticity region appears to be a strip bounded by a natural boundary, whose width tends linearly to 0 as tends to the critical value.     

Continuous spectrum in the ground state of two spin-1/2 models in the infinite-volume limit

We show that in the ground states of the infinite-volume limits of both the spin-1/2 anisotropic antiferromagnetic Heisenberg model (in dimensions d2), and the ferromagnetic Ising model in a strong transverse field (in dimensions d1) there is an interval in the spectrum above the mass gap which contains a continuous band of energy levels. We use the methods of Bricmont and Fröhlich to develop our expansions, as well as a method of Kennedy and Tasaki to do the expansions in the quantum mechanical limit. Where the expansions converge, they are then shown to have spectral measures which have absolutely continuous parts on intervals above the mass gaps.     

Super-effective-field theory and coherent-anomaly method in cooperative phenomena

Conceptual arguments on the coherent-anomaly method (CAM) and on the super-effective-field theory are presented to explain the basic ideas of these theories. Some possible applications are also suggested.