Researcher profile

Vesselin Tonchev

Vesselin Tonchev contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
12works
0followers
9topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

12 published item(s)

preprint2022arXiv

Predicted universality class of step bunching found on DC-heated Si(111) surfaces

Concerted experimental and numerical studies of step bunching on vicinal crystal surfaces resulting from step-down electromigration of partially charged adatoms, confirmed the theoretical prediction of scaling dependence of the minimal bunch distance $l_{\rm min}$ on the bunch size $N$: $l_{\rm min} \sim N^{-γ}$, with $γ= 2/3$. The value of the so called size-scaling exponent $γ$ was observed in experiments on vicinal surfaces of semiconducting, metallic, and dielectric materials. Careful theoretical investigations and numerical calculations predict a second value of $γ= 1/2$. However, this value is still not been reported from experiments. And we report here experimental observation of step bunching in the universality class relative to $γ= 1/2$. This is achieved by monitoring step flow during sublimation of Si(111)-vicinals heated by a direct step-down current at ~1200$^\circ$C. In the experiment we also measure other characteristic for the bunching quantities, such as the mean total number of steps in the bunch $N$ and the mean bunch width $W$. We then compare our findings with published experimental and numerical data to arrive at a theoretically consistent framework in terms of universality classes. The ultimate benefit of our study is not only to advance fundamental knowledge but also to provide further guidance for bottom-up synthesis of vicinal nanotemplates.

preprint2020arXiv

Bunch width versus macrostep height: A quantitative study of the effects of step-step repulsion

Bunching of steps at the surface of growing crystals can be induced by both directions of the driving force: step up and step down. The processes happen in different adatom concentrations and differ in character. In this study we show how the overall picture of the bunching process depends on the strength of short range step-step repulsion. The repulsive interaction between steps, controlled by an additional parameter, is introduced into the recently studied atomistic scale model of vicinal crystal growth, based on cellular automata. It is shown that the repulsion modifies bunching process in a different way, depending on the direction of the destabilizing force. In particular, bunch profiles, stability diagrams and time-scaling dependences of various bunch properties are affected when the step-step repulsion increases. The repulsion between steps creates a competition between two characteristic sizes - bunch width and macrostep height, playing the role of the second length scale that describes the step bunching phenomenon. A new characteristic time scale dependent on the step-step repulsion parameter emerges as an effect of interplay between (01) faceted macrosteps and (11) faceted bunches. The bunch height being the major characteristic size of the bunches is not influenced dramatically by the repulsion.

preprint2015arXiv

Maximal Density, Kinetics of Deposition and Percolation Threshold of Loose Packed Lattices

In many areas of research it is interesting how lattices can be filled with particles that have no nearest neighbors, or they are in limited quantities. Examples may be found in statistical physics, chemistry, materials science, discrete mathematics, etc. Using Monte Carlo (MC) simulation we study the kinetics of filling of square lattice (2D). Two complementary rules are used to fill the lattice. We study their influence on the kinetics of the process as well as on the properties of the obtained systems. According to the first rule the occupied sites may not share edges (nearest neighbors occupations are not permitted). Under this condition, the maximum possible concentration is 0.5, forming a checkerboard type structure. However, we found that if the deposition is done by random selection of sites the concentration of 0.5 is inaccessible and the maximum concentration is Cmax(2D)=0.3638 (0.0003) for 2D lattice. If the lattice is 3D we find that the maximal concentration is even lower Cmax(3D)=0.326 (0.001). The second rule establishes permission to break the first one with certain probability 0<=p<=1, thus the occupied sites can start to share edges when p>0. In this case higher then 0.3638 concentrations are accessible, even up to C=1. Therefore the percolation threshold Pc can be reached. Its value depends on the value of the probability p. Our model describes the kinetics of formation of thin films of particles attracted by the substrate but repulsing each other.

preprint2015arXiv

Time Scaling Relations for Step Bunches from Models with Step-Step Attractions (B1-Type Models)

The step bunching instability is studied in three models of step motion defined in terms of ordinary differential equations (ODE). The source of instability in these models is step-step attraction, it is opposed by step-step repulsion and the developing surface patterns reflect the balance between the two. The first model, TE2, is a generalization of the seminal model of Tersoff et al. (1995). The second one, LW2, is obtained from the model of Liu and Weeks (1998) using the repulsions term to construct the attractions one with retained possibility to change the parameters in the two independently. The third model, MM2, is a minimal one constructed ad hoc and in this article it plays a central role. New scheme for scaling the ODE in vicinal studies is applied towards deciphering the pre-factors in the time-scaling relations. In all these models the patterned surface is self-similar - only one length scale is necessary to describe its evolution (hence B1-type). The bunches form finite angles with the terraces. Integrating numerically the equations for step motion and changing systematically the parameters we obtain the overall dependence of time-scaling exponent \b{eta} on the power of step-step attractions p as \b{eta} = 1/(3+p) for MM2 and hypothesize based on restricted set of data that it is \b{eta} = 1/(5+p) for LW2 and TE2.

preprint2015arXiv

Unstable vicinal crystal growth from cellular automata

In order to study the unstable step motion on vicinal crystal surfaces we devise vicinal Cellular Automata. Each cell from the colony has value equal to its height in the vicinal, initially the steps are regularly distributed. Another array keeps the adatoms, initially distributed randomly over the surface. The growth rule defines that each adatom at right nearest neighbor position to a (multi-) step attaches to it. The update of whole colony is performed at once and then time increases. This execution of the growth rule is followed by compensation of the consumed particles and by diffusional update(s) of the adatom population. Two principal sources of instability are employed: biased diffusion and infinite inverse Ehrlich-Schwoebel barrier (iiSE). Since these factors are not opposed by step-step repulsion the formation of multi-steps is observed but in general the step bunches preserve a finite width. We monitor the developing surface patterns and quantify the observations by scaling laws with focus on the eventual transition from diffusion-limited to kinetics-limited phenomenon. The time-scaling exponent of the bunch size N is 1/2 for the case of biased diffusion and 1/3 for the case of iiSE. Additional distinction is possible based on the time-scaling exponents of the sizes of multi-steps, these are 0.36-0.4 (for biased diffusion) and 1/4 (iiSE).

preprint2012arXiv

Classification of step bunching phenomena

The classification of bunching of straight steps on vicinal crystal surfaces identifies two types according to the behavior of the minimal step-step distance in the bunch lmin with increasing the number of steps N in it. In the B1-type lmin remains constant while in the B2-type it decreases. Both types are illustrated by new results for well-known models. The precise numerical analysis is aimed at the intermediate asymptotic regime where self-similar spatiotemporal patterns develop. In the model of Tersoff et al. the regular step train is destabilized by step-step attraction of infinite range. It is shown that this model belongs to the B1-type and the same time-scaling exponent of 1/5 for N, terrace width and bunch width is obtained. An extended set of scaling exponents is obtained from the model of S.Stoyanov of diffusion-limited evaporation affected by electromigration of the adatoms. This model is of B2-type and shows a systematic shift of the exponents with respect to the predictions of the hypothesis for universality classes in bunching thus requiring further modification of it.

preprint2012arXiv

Mesoscaled Patterning In Nanostructured Thin Films/Silicate Glass Systems

Various nanostructured thin solid films on silicate based glass plates commonly used as substrates in contemporary vacuum technologies are studied. Evolution of mesostructured spiral and spiral crack patterns as a result of shrinkage or crystalliz4 tion processes is demonstrated. It is established that the surface hydrated layer of the glass is liable for the patterning observed. Evidence is presented supporting a mechanism of selforganization phenomena in the systems studied. It is found that nano-sized sol gel silica thin film can be successfully used as a model system of the surface hydrated layer. Comparison of different thin solid films on silicate glass plates shows contour and dimensional similarity of memory patterns thus pointing out the occurrence of a general phenomenon.

preprint2012arXiv

Step Bunching In Conserved Systems: Scaling And Universality

We study the step bunching process in three different 1D step flow models and obtain scaling relations for the step bunches formed in the long times limit. The first one was introduced by S.Stoyanov [Jap. J.Appl. Phys. 29, (1990) L659] as the simplest &#39;realistic&#39; model of step bunching due to drift of the adatoms. Here we show that it could lead to (at least) two different types of step bunching, depending on the magnitude of the drift. The other two models are minimal models: the equations for step velocity are constructed ad hoc from two terms with opposite effects - destabilizing and, respectively, stabilizing the regular step train.

preprint2011arXiv

Diffusion Limited Aggregation with modified local rules

Results from a modified Diffusion Limited Aggregation (DLA) model are presented. The modifications of the classical DLA model are in the attachment to the cluster rules and in the scheme of particle generation/killing. In the classical DLA model if a particle reaches the growing cluster it sticks to it immediately and irreversibly and then the next particle is released. We will abandon this original prerequisite, and by changing the sticking probability to the cluster we will change the diffusion regime towards more kinetic one. For a growing cluster consisting of only one type of particles this variation in the sticking probability is (more or less) a rude violation of the hypothesis for diffusion limitation in the DLA model. Since in a lot of experiments different types of particles are used with different sticking probabilities (e.g. different regimes of attachment), we develop a modified DLA model with two types of particles. The second modification we introduce at that point is a scheme for particle generation/killing we call &#34;second chance&#34; - when a particle is killed after reaching a given limiting distance from the cluster, it is killed and then returned to the point it was originally generated. Thus the model is capable to produce a great variety of growing patterns (fractals, spirals) by changing only a single parameter and we are able to construct a morphological diagram of our generalized DLA model with two different types of particles.

preprint2011arXiv

Hybrid Models of Step Bunching

We introduce two hybrid models of step bunching on vicinal crystal surfaces. The model equations for step velocity are constructed by the two possible exchanges of terms between the equations of two primary models MM2 and LW2 [arXiv:1011.1863], both showing the specific type of bunching with minimal step-step distance lmin in the bunch independent of the number of steps N in it. This feature is preserved only in the hybrid model LW2MM (the first term in the model equation comes from LW2 and the second one - from MM2) but in a rather complex fashion -- the surface slope is largest in the both ends of the bunch and after a sharp decrease jumps again to become constant in the inner part. We restrict our considerations to the simplest case of p = 0, p being the exponent in the destabilizing term in the velocity equations. The time-scaling exponent of N in LW2MM is ~1/3 and is independent of n, the exponent in the stabilizing term of the velocity equations. The other model, MM2LW, shows an interesting type of step bunching -- some bunches grow to a certain size and then decay emitting steps towards the two adjacent bunches. The bunch compression with the increase of N is pronounced.

preprint2010arXiv

Models of step bunching: Turning repulsion into attraction

We report numerical results for two models of vicinal motion. The first, LW, aims at crystal evaporation when the detachment from steps is slow [Liu and Weeks, PRB 57, 23 (1998) 14891]. The source of destabilization is electromigration force acting on the adatoms. The destabilizing part of equation(s) of the step velocity is linear in the widths of the adjacent terraces with larger contribution of the terrace behind. This asymmetry is controlled by a single parameter b. The stabilization part accounts for the tendency to equidistant spacing dictated by the interstep repulsions. We construct the second model, LW2, from LW in the same manner as was constructed Minimal Model 2 (MM2) from another minimal model [B.Ranguelov et al., Nanoscience and Nanotechnology 6, (2006) 31] - keeping the &#39;repulsions term&#39; from LW and introducing a similar one with opposite sign as &#39;attractions term&#39;. For LW we obtain for first time that in the pre-factor of the time scaling of the number of steps in the bunch N enters only the parameter b. In LW2 we find the same type of step bunching as in MM2 - the surface slope in the bunch is constant and not function of N. Further, we obtain the time scaling of N with exponent ~0.18, found also in experiments, and the minimal distance in the bunch as function of the model parameters.

preprint2010arXiv

Self-Organization on Unstable Vicinal Surfaces with Competing Interactions

Long-ranged step-step attractions destabilize the vicinal crystal surfaces. Their competition with shorter-ranged step-step repulsions results in self-organized patterns. The exponent in the time-scaling of their characteristic size is influenced only by the range of the attractions but not by the range of the repulsions as we show based on precise numerical analysis of two different minimal models. Another anisotropy we identify is that the vicinal surface is not destabilized by shorter-ranged attractions.