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

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

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

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

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

Modeling one-dimensional island growth with mass-dependent detachment rates

We study one-dimensional models of particle diffusion and attachment/detachment from islands where the detachment rates gamma(m) of particles at the cluster edges increase with cluster mass m. They are expected to mimic the effects of lattice mismatch with the substrate and/or long-range repulsive interactions that work against the formation of long islands. Short-range attraction is represented by an overall factor epsilon<<1 in the detachment rates relatively to isolated particle hopping rates [epsilon ~ exp(-E/T), with binding energy E and temperature T]. We consider various gamma(m), from rapidly increasing forms such as gamma(m) ~ m to slowly increasing ones, such as gamma(m) ~ [m/(m+1)]^b. A mapping onto a column problem shows that these systems are zero-range processes, whose steady states properties are exactly calculated under the assumption of independent column heights in the Master equation. Simulation provides island size distributions which confirm analytic reductions and are useful whenever the analytical tools cannot provide results in closed form. The shape of island size distributions can be changed from monomodal to monotonically decreasing by tuning the temperature or changing the particle density rho. Small values of the scaling variable X=epsilon^{-1}rho/(1-rho) favour the monotonically decreasing ones. However, for large X, rapidly increasing gamma(m) lead to distributions with peaks very close to and rapidly decreasing tails, while slowly increasing gamma(m) provide peaks close to /2$ and fat right tails.Comment: 16 pages, 6 figure