Kinetic modelling of epitaxial film growth with up- and downward step barriers

The formation of three-dimensional structures during the epitaxial growth of films is associated to the reflection of diffusing particles in descending terraces due to the presence of the so-called Ehrlich-Schwoebel (ES) barrier. We generalize this concept in a solid-on-solid growth model, in which a barrier dependent on the particle coordination (number of lateral bonds) exists whenever the particle performs an interlayer diffusion. The rules do not distinguish explicitly if the particle is executing a descending or an ascending interlayer diffusion. We show that the usual model, with a step barrier in descending steps, produces spurious, columnar, and highly unstable morphologies if the growth temperature is varied in a usual range of mound formation experiments. Our model generates well-behaved mounded morphologies for the same ES barriers that produce anomalous morphologies in the standard model. Moreover, mounds are also obtained when the step barrier has an equal value for all particles independently if they are free or bonded. Kinetic roughening is observed at long times, when the surface roughness w and the characteristic length $\xi$ scale as $w ~ t^\beta$ and $\xi ~ t^\zeta$ where $\beta \approx 0.31$ and $\zeta \approx 0.22$, independently of the growth temperature.Comment: 15 pages, 7 figure

Roughness exponents and grain shapes

In surfaces with grainy features, the local roughness $w$ shows a crossover at a characteristic length $r_c$, with roughness exponent changing from $\alpha_1\approx 1$ to a smaller $\alpha_2$. The grain shape, the choice of $w$ or height-height correlation function (HHCF) $C$, and the procedure to calculate root mean-square averages are shown to have remarkable effects on $\alpha_1$. With grains of pyramidal shape, $\alpha_1$ can be as low as 0.71, which is much lower than the previous prediction 0.85 for rounded grains. The same crossover is observed in the HHCF, but with initial exponent $\chi_1\approx 0.5$ for flat grains, while for some conical grains it may increase to $\chi_1\approx 0.7$. The universality class of the growth process determines the exponents $\alpha_2=\chi_2$ after the crossover, but has no effect on the initial exponents $\alpha_1$ and $\chi_1$, supporting the geometric interpretation of their values. For all grain shapes and different definitions of surface roughness or HHCF, we still observe that the crossover length $r_c$ is an accurate estimate of the grain size. The exponents obtained in several recent experimental works on different materials are explained by those models, with some surface images qualitatively similar to our model films.Comment: 7 pages, 6 figures and 2 table

Phase transitions in dependence of apex predator decaying ratio in a cyclic dominant system

Cyclic dominant systems, like rock-paper-scissors game, are frequently used to explain biodiversity in nature, where mobility, reproduction and intransitive competition are on stage to provide the coexistence of competitors. A significantly new situation emerges if we introduce an apex predator who can superior all members of the mentioned three-species system. In the latter case the evolution may terminate into three qualitatively different destinations depending on the apex predator decaying ratio $q$. In particular, the whole population goes extinct or all four species survive or only the original three-species system remains alive as we vary the control parameter. These solutions are separated by a discontinuous and a continuous phase transitions at critical $q$ values. Our results highlight that cyclic dominant competition can offer a stable way to survive even in a predator-prey-like system that can be maintained for large interval of critical parameter values.Comment: version to appear in EPL. 7 pages, 7 figure

Invasion controlled pattern formation in a generalized multi-species predator-prey system

Rock-scissors-paper game, as the simplest model of intransitive relation between competing agents, is a frequently quoted model to explain the stable diversity of competitors in the race of surviving. When increasing the number of competitors we may face a novel situation because beside the mentioned unidirectional predator-prey-like dominance a balanced or peer relation can emerge between some competitors. By utilizing this possibility in the present work we generalize a four-state predator-prey type model where we establish two groups of species labeled by even and odd numbers. In particular, we introduce different invasion probabilities between and within these groups, which results in a tunable intensity of bidirectional invasion among peer species. Our study reveals an exceptional richness of pattern formations where five quantitatively different phases are observed by varying solely the strength of the mentioned inner invasion. The related transition points can be identified with the help of appropriate order parameters based on the spatial autocorrelation decay, on the fraction of empty sites, and on the variance of the species density. Furthermore, the application of diverse, alliance-specific inner invasion rates for different groups may result in the extinction of the pair of species where this inner invasion is moderate. These observations highlight that beyond the well-known and intensively studied cyclic dominance there is an additional source of complexity of pattern formation that has not been explored earlier.Comment: 8 pages, 8 figures. To appear in PR

Scaling Invariance in a Time-Dependent Elliptical Billiard

We study some dynamical properties of a classical time-dependent elliptical billiard. We consider periodically moving boundary and collisions between the particle and the boundary are assumed to be elastic. Our results confirm that although the static elliptical billiard is an integrable system, after to introduce time-dependent perturbation on the boundary the unlimited energy growth is observed. The behaviour of the average velocity is described using scaling arguments

Fermi acceleration and suppression of Fermi acceleration in a time-dependent Lorentz Gas

We study some dynamical properties of a Lorentz gas. We have considered both the static and time dependent boundary. For the static case we have shown that the system has a chaotic component characterized with a positive Lyapunov Exponent. For the time-dependent perturbation we describe the model using a four-dimensional nonlinear map. The behaviour of the average velocity is considered in two situations (i) non-dissipative and (ii) dissipative. Our results show that the unlimited energy growth is observed for the non-dissipative case. However, when dissipation, via damping coefficients, is introduced the senary changes and the unlimited engergy growth is suppressed. The behaviour of the average velocity is described using scaling approach