A simple acceleration algorithm for Ising systems with conserved magnetization (model B) is presented. The dynamical critical exponent and the domain growth exponent in spinodal decomposition are found to be equal to those observed for model A systems with no conservation law. Our results demonstrate that systems with global conservation laws are in the same dynamical universality class as systems with no conservation laws. PACS numbers: 64.60.Fr, 05.50.+q, 64.60.Ht, 75.40.MgThere are three principal themes in this Letter. First, a new acceleration algorithm is proposed for Ising spin systems having a conserved magnetization (model B) [1,2]. Physical arguments (for a one-dimensional system) and detailed computer simulations (in two dimensions) show that the dynamic critical exponent z (defined by the relation r-L z , where r is the correlation time and L the linear dimension of the system at criticality) for this algorithm is the same as that obtained for model A, a system with no conservation laws. Second, we have used our algorithm to study the dynamics of spinodal decomposition [3] of a 2D Ising system with conserved magnetization following the quench from a high-temperature homogeneous phase into the two-phase coexistence region. We find that in the late-time scaling regime, the exponent n characterizing the growth of the characteristic domain size R (R~~t n ) is the same as that for a
model A system («•{) and different from the conventional model B result (n = j ). Third, our results are consistent with recent predictions by Bray based on renormalization-group arguments [4-6]. They thus have a bearing on the recent controversy between Bray and Tamayo and Klein [7,8] (TK). We have studied one of the two versions of the TK model and find that the TK exponent is equal to the model A result, in contradiction to their 2D computersimulation results. Our results show that global conservation laws are irrelevant in determining the dynamical universality class, as argued recently by Bray [4-6].Our algorithm for the Ising system is a generalization of Glauber dynamics for model A-it consists of single spin flips governed by the usual Metropolis [9] rules. The global conservation of magnetization at a desired value Mo (model B) is enforced by a Creutz "demon" or bag [10]. Spin flips are only allowed if the total sample magnetization after the spin flip, Af, lies in the range A/o~£< M< Mo+8. Thus our algorithm smoothly extrapolates from model A (unbounded 8) to model B (5=0). In practice, all our calculations have been carried out with A/o=0, <5 = 2. Our principal result is that the exponents z and n are independent of whether 8 = 2 or 8 is unbounded, and correspond to the model A results. Note that 5 = 2 in the thermodynamic limit is a conserved