Cluster growth in far-from-equilibrium particle models with diffusion, detachment, reattachment and deposition

Monolayer cluster growth in far-from-equilibrium systems is investigated by applying simulation and analytic techniques to minimal hard core particle (exclusion) models. The first model (I), for post-deposition coarsening dynamics, contains mechanisms of diffusion, attachment, and slow activated detachment (at rate epsilon<<1) of particles on a line. Simulation shows three successive regimes of cluster growth: fast attachment of isolated particles; detachment allowing further (epsilon t)^(1/3) coarsening of average cluster size; and t^(-1/2) approach to a saturation size going like epsilon^(-1/2). Model II generalizes the first one in having an additional mechanism of particle deposition into cluster gaps, suppressed for the smallest gaps. This model exhibits early rapid filling, leading to slowing deposition due to the increasing scarcity of deposition sites, and then continued power law (epsilon t)^(1/2) cluster size coarsening through the redistribution allowed by slow detachment. The basic (epsilon t)^(1/3) domain growth laws and epsilon^(-1/2) saturation in model I are explained by a simple scaling picture. A second, fuller approach is presented which employs a mapping of cluster configurations to a column picture and an approximate factorization of the cluster configuration probability within the resulting master equation. This allows quantitative results for the saturation of model I in excellent agreement with the simulation results. For model II, it provides a one-variable scaling function solution for the coarsening probability distribution, and in particular quantitative agreement with the cluster length scaling and its amplitude.Comment: Accepted in Phys. Rev. E; 9 pages with figure

Length and time scale divergences at the magnetization-reversal transition in the Ising model

The divergences of both the length and time scales, at the magnetization- reversal transition in Ising model under a pulsed field, have been studied in the linearized limit of the mean field theory. Both length and time scales are shown to diverge at the transition point and it has been checked that the nature of the time scale divergence agrees well with the result obtained from the numerical solution of the mean field equation of motion. Similar growths in length and time scales are also observed, as one approaches the transition point, using Monte Carlo simulations. However, these are not of the same nature as the mean field case. Nucleation theory provides a qualitative argument which explains the nature of the time scale growth. To study the nature of growth of the characteristic length scale, we have looked at the cluster size distribution of the reversed spin domains and defined a pseudo-correlation length which has been observed to grow at the phase boundary of the transition.Comment: 9 pages Latex, 3 postscript figure

Domain wall theory and non-stationarity in driven flow with exclusion

We study the dynamical evolution toward steady state of the stochastic non-equilibrium model known as totally asymmetric simple exclusion process, in both uniform and non-uniform (staggered) one-dimensional systems with open boundaries. Domain-wall theory and numerical simulations are used and, where pertinent, their results are compared to existing mean-field predictions and exact solutions where available. For uniform chains we find that the inclusion of fluctuations inherent to the domain-wall formulation plays a crucial role in providing good agreement with simulations, which is severely lacking in the corresponding mean-field predictions. For alternating-bond chains the domain-wall predictions for the features of the phase diagram in the parameter space of injection and ejection rates turn out to be realized only in an incipient and quantitatively approximate way. Nevertheless, significant quantitative agreement can be found between several additional domain-wall theory predictions and numerics.Comment: 12 pages, 12 figures (published version

Logarithmic coarsening and glassy behavior in a polymer model with mass-dependent diffusion

We present a model of polymer growth and diffusion with frustration mechanisms for density increase and with diffusion rates of Arrhenius form with mass-dependent energy barriers Gamma(m) ~ (m-1)^gamma. It shows non-universal logarithmic coarsening involving the exponent gamma. Strong-glass behavior is found in the typical times for disappearance of all polymers up to a given length, without reference to the equilibrium states of the macroscopic system. These features are predicted by numerical simulations, scaling theories and an analytic solution of the master equation within an independent interval approximation, which also provides the cluster size distribution.Comment: 16 pages, including 7 figures. To be published in Phys. Rev.

Correlation--function distributions at the Nishimori point of two-dimensional Ising spin glasses

The multicritical behavior at the Nishimori point of two-dimensional Ising spin glasses is investigated by using numerical transfer-matrix methods to calculate probability distributions $P(C)$ and associated moments of spin-spin correlation functions $C$ on strips. The angular dependence of the shape of correlation function distributions $P(C)$ provides a stringent test of how well they obey predictions of conformal invariance; and an even symmetry of $(1-C) P(C)$ reflects the consequences of the Ising spin-glass gauge (Nishimori) symmetry. We show that conformal invariance is obeyed in its strictest form, and the associated scaling of the moments of the distribution is examined, in order to assess the validity of a recent conjecture on the exact localization of the Nishimori point. Power law divergences of $P(C)$ are observed near C=1 and C=0, in partial accord with a simple scaling scheme which preserves the gauge symmetry.Comment: Final version to be published in Phys Rev

Non-universal coarsening and universal distributions in far-from equilibrium systems

Anomalous coarsening in far-from equilibrium one-dimensional systems is investigated by simulation and analytic techniques. The minimal hard core particle (exclusion) models contain mechanisms of aggregated particle diffusion, with rates epsilon<<1, particle deposition into cluster gaps, but suppressed for the smallest gaps, and breakup of clusters which are adjacent to large gaps. Cluster breakup rates vary with the cluster length x as kx^alpha. The domain growth law x ~ (epsilon t)^z, with z=1/(2+alpha) for alpha>0, is explained by a scaling picture, as well as the scaling of the density of double vacancies (at which deposition and cluster breakup are allowed) as 1/[t(epsilon t)^z]. Numerical simulations for several values of alpha and epsilon confirm these results. An approximate factorization of the cluster configuration probability is performed within the master equation resulting from the mapping to a column picture. The equation for a one-variable scaling function explains the above results. The probability distributions of cluster lengths scale as P(x)= 1/(epsilon t)^z g(y), with y=x/(epsilon t)^z. However, those distributions show a universal tail with the form g(y) ~ exp(-y^{3/2}), which disagrees with the prediction of the independent cluster approximation. This result is explained by the connection of the vacancy dynamics with the problem of particle trapping in an infinite sea of traps and is confirmed by simulation.Comment: 30 pages (10 figures included), to appear in Phys. Rev.

Scaling treatment of the random field Ising model

Analytic phenomenological scaling is carried out for the random field Ising model in general dimensions using a bar geometry. Domain wall configurations and their decorated profiles and associated wandering and other exponents $(\zeta,\gamma,\delta,\mu)$ are obtained by free energy minimization. Scaling between different bar widths provides the renormalization group (RG) transformation. Its consequences are (1) criticality at $h=T=0$ in $d \leq 2$ with correlation length $\xi(h,T)$ diverging like $\xi(h,0) \propto h^{-2/(2-d)}$ for $d<2$ and $\xi(h,0) \propto \exp[1/(c_1\gamma h^{\gamma})]$ for $d=2$, where $c_1$ is a decoration constant; (2) criticality in $d = 2+\epsilon$ dimensions at $T=0$, $h^{\ast}= (\epsilon/2c_1)^{1/\gamma}$, where $\xi \propto [(s-s^{\ast})/s]^{-2\epsilon/\gamma}$, $s \equiv h^{\gamma}$. Finite temperature generalizations are outlined. Numerical transfer matrix calculations and results from a ground state algorithm adapted for strips in $d=2$ confirm the ingredients which provide the RG description.Comment: RevTex v3.0, 5 pages, plus 4 figures uuencode