Comment on ``Scaling Laws for a System with Long-Range Interactions within Tsallis Statistics''

In their recent Letter [Phys. Rev. Lett. 83, 4233 (1999)], Salazar and Toral (ST) study numerically a finite Ising chain with non-integrable interactions decaying like 1/r^(d+sigma) where -d <= sigma <= 0 (like ST, we discuss general dimensionality d). In particular, they explore a presumed connection between non-integrable interactions and Tsallis's non-extensive statistics. We point out that (i) non-integrable interactions provide no more motivation for Tsallis statistics than do integrable interactions, i.e., Gibbs statistics remain meaningful for the non-integrable case, and in fact provide a {\em complete and exact treatment}; and (ii) there are undesirable features of the method ST use to regulate the non-integrable interactions.Comment: Accepted for publication in Phys. Rev. Let

Single-Species Three-Particle Reactions in One Dimension

Renormalization group calculations for fluctuation-dominated reaction-diffusion systems are generally in agreement with simulations and exact solutions. However, simulations of the single-species reactions 3A->(0,A,2A) at their upper critical dimension d_c=1 have found asymptotic densities argued to be inconsistent with renormalization group predictions. We show that this discrepancy is resolved by inclusion of the leading corrections to scaling, which we derive explicitly and show to be universal, a property not shared by the A+A->(0,A) reactions. Finally, we demonstrate that two previous Smoluchowski approaches to this problem reduce, with various corrections, to a single theory which yields, surprisingly, the same asymptotic density as the renormalization group.Comment: 8 pages, 5 figs, minor correction

Fast and Accurate Coarsening Simulation with an Unconditionally Stable Time Step

We present Cahn-Hilliard and Allen-Cahn numerical integration algorithms that are unconditionally stable and so provide significantly faster accuracy-controlled simulation. Our stability analysis is based on Eyre's theorem and unconditional von Neumann stability analysis, both of which we present. Numerical tests confirm the accuracy of the von Neumann approach, which is straightforward and should be widely applicable in phase-field modeling. We show that accuracy can be controlled with an unbounded time step Delta-t that grows with time t as Delta-t ~ t^alpha. We develop a classification scheme for the step exponent alpha and demonstrate that a class of simple linear algorithms gives alpha=1/3. For this class the speed up relative to a fixed time step grows with the linear size of the system as N/log N, and we estimate conservatively that an 8192^2 lattice can be integrated 300 times faster than with the Euler method.Comment: 14 pages, 6 figure

Anisotropic Coarsening: Grain Shapes and Nonuniversal Persistence

We solve a coarsening system with small but arbitrary anisotropic surface tension and interface mobility. The resulting size-dependent growth shapes are significantly different from equilibrium microcrystallites, and have a distribution of grain sizes different from isotropic theories. As an application of our results, we show that the persistence decay exponent depends on anisotropy and hence is nonuniversal.Comment: 4 pages (revtex), 2 eps figure

Applications of Field-Theoretic Renormalization Group Methods to Reaction-Diffusion Problems

We review the application of field-theoretic renormalization group (RG) methods to the study of fluctuations in reaction-diffusion problems. We first investigate the physical origin of universality in these systems, before comparing RG methods to other available analytic techniques, including exact solutions and Smoluchowski-type approximations. Starting from the microscopic reaction-diffusion master equation, we then pedagogically detail the mapping to a field theory for the single-species reaction k A -> l A (l < k). We employ this particularly simple but non-trivial system to introduce the field-theoretic RG tools, including the diagrammatic perturbation expansion, renormalization, and Callan-Symanzik RG flow equation. We demonstrate how these techniques permit the calculation of universal quantities such as density decay exponents and amplitudes via perturbative eps = d_c - d expansions with respect to the upper critical dimension d_c. With these basics established, we then provide an overview of more sophisticated applications to multiple species reactions, disorder effects, L'evy flights, persistence problems, and the influence of spatial boundaries. We also analyze field-theoretic approaches to nonequilibrium phase transitions separating active from absorbing states. We focus particularly on the generic directed percolation universality class, as well as on the most prominent exception to this class: even-offspring branching and annihilating random walks. Finally, we summarize the state of the field and present our perspective on outstanding problems for the future.Comment: 10 figures include

Anomalous dimension in a two-species reaction–diffusion system

We study a two-species reaction–diffusion system with the reactions and , with general diffusion constants D A and D B . Previous studies showed that for dimensions the B particle density decays with a nontrivial, universal exponent that includes an anomalous dimension resulting from field renormalization. We demonstrate via renormalization group methods that the scaled B particle correlation function has a distinct anomalous dimension resulting in the asymptotic scaling , where the exponent results from the renormalization of the square of the field associated with the B particles. We compute this exponent to first order in , a calculation that involves 61 Feynman diagrams, and also determine the logarithmic corrections at the upper critical dimension . Finally, we determine the exponent numerically utilizing a mapping to a four-walker problem for the special case of A particle coalescence in one spatial dimension