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

Driven flow with exclusion and transport in graphene-like structures

The totally asymmetric simple exclusion process (TASEP), a well-known model in its strictly one-dimensional (chain) version, is generalized to cylinder (nanotube) and ribbon (nanoribbon) geometries. A mean-field theoretical description is given for very narrow ribbons ("necklaces"), and nanotubes. For specific configurations of bond transmissivity rates, and for a variety of boundary conditions, theory predicts equivalent steady state behavior between (sublattices on) these structures and chains. This is verified by numerical simulations, to excellent accuracy, by evaluating steady-state currents. We also numerically treat ribbons of general width. We examine the adequacy of this model to the description of electronic transport in carbon nanotubes and nanoribbons, or specifically-designed quantum dot arrays.Comment: RevTeX, 13 pages, 9 figures (published version

Mean field and Monte Carlo studies of the magnetization-reversal transition in the Ising model

Detailed mean field and Monte Carlo studies of the dynamic magnetization-reversal transition in the Ising model in its ordered phase under a competing external magnetic field of finite duration have been presented here. Approximate analytical treatment of the mean field equations of motion shows the existence of diverging length and time scales across this dynamic transition phase boundary. These are also supported by numerical solutions of the complete mean field equations of motion and the Monte Carlo study of the system evolving under Glauber dynamics in both two and three dimensions. Classical nucleation theory predicts different mechanisms of domain growth in two regimes marked by the strength of the external field, and the nature of the Monte Carlo phase boundary can be comprehended satisfactorily using the theory. The order of the transition changes from a continuous to a discontinuous one as one crosses over from coalescence regime (stronger field) to nucleation regime (weaker field). Finite size scaling theory can be applied in the coalescence regime, where the best fit estimates of the critical exponents are obtained for two and three dimensions.Comment: 16 pages latex, 13 ps figures, typos corrected, references adde

Two-dimensional Site-Bond Percolation as an Example of Self-Averaging System

The Harris-Aharony criterion for a statistical model predicts, that if a specific heat exponent $\alpha \ge 0$, then this model does not exhibit self-averaging. In two-dimensional percolation model the index $\alpha=-{1/2}$. It means that, in accordance with the Harris-Aharony criterion, the model can exhibit self-averaging properties. We study numerically the relative variances $R_{M}$ and $R_{\chi}$ for the probability $M$ of a site belongin to the "infinite" (maximum) cluster and the mean finite cluster size $\chi$. It was shown, that two-dimensional site-bound percolation on the square lattice, where the bonds play the role of impurity and the sites play the role of the statistical ensemble, over which the averaging is performed, exhibits self-averaging properties.Comment: 15 pages, 5 figure

Duality and phase diagram of one dimensional transport

The observation of duality by Mukherji and Mishra in one dimensional transport problems has been used to develop a general approach to classify and characterize the steady state phase diagrams. The phase diagrams are determined by the zeros of a set of coarse-grained functions without the need of detailed knowledge of microscopic dynamics. In the process, a new class of nonequilibrium multicritical points has been identified.Comment: 6 pages, 2 figures (4 eps files

Smoothly-varying hopping rates in driven flow with exclusion

We consider the one-dimensional totally asymmetric simple exclusion process (TASEP) with position-dependent hopping rates. The problem is solved,in a mean field/adiabatic approximation, for a general (smooth) form of spatial rate variation. Numerical simulations of systems with hopping rates varying linearly against position (constant rate gradient), for both periodic and open boundary conditions, provide detailed confirmation of theoretical predictions, concerning steady-state average density profiles and currents, as well as open-system phase boundaries, to excellent numerical accuracy.Comment: RevTeX 4.1, 14 pages, 9 figures (published version

Diffusion Enhances Chirality Selection

Diffusion effect on chirality selection in a two-dimensional reaction-diffusion model is studied by the Monte Carlo simulation. The model consists of achiral reactants A which turn into either of the chiral products, R or S, in a solvent of chemically inactive vacancies V. The reaction contains the nonlinear autocatalysis as well as recycling process, and the chiral symmetry breaking is monitored by an enantiomeric excess $\phi$. Without dilution a strong nonlinear autocatalysis ensures chiral symmetry breaking. By dilution, the chiral order $\phi$ decreases, and the racemic state is recovered below the critical concentration $c_c$. Diffusion effectively enhances the concentration of chiral species, and $c_c$ decreases as the diffusion coefficient $D$ increases. The relation between $\phi$ and $c$ for a system with a finite $D$ fits rather well to an interpolation formula between the diffusionless(D=0) and homogeneous ($D=\infty$) limits.Comment: 7 pages, 6 figure