Mathematical Modeling of Phase Separation and Branching Process of the Film Structure during Binary Thin Film Deposition

: In this study, we applied a mathematical model to explore the mechanism and factors leading to phase separation and the formation of branching structures with nanocolumns extending from larger clusters formed on the substrate of a grown ﬁlm. The mathematical model simulated the growth of a thin ﬁlm over time by using partial differential equations, including the processes of adsorption, phase separation, and diffusion due to the curvature of the thin ﬁlm surface. The modeling results revealed the possible mechanism that could lead to the formation of the described branching structures. That mechanism can be divided into two main steps. The ﬁrst step is the growth of a relatively large cluster (of a component that makes up the branching phase) on the substrate during the initial growth stages. The second step is the division process of that large cluster into smaller clusters in the later growth stages. The model parameters inﬂuencing the growth conditions that lead to the formation mechanism of the branching structures were determined, and their inﬂuences on the phase structure were analyzed. the substrate material in the early growth stages. The further increase in the diffusion coefﬁcient D AS from 8.8 × 10 − 18 m 2 /s to 8.8 × 10 − 17 m 2 /s resulted in an increase in the average cluster size (of component A) from 9.5 to 18 nm; the average cluster sizes in that range were favorable for the branching process to occur.


Introduction
Nanocomposites, which are multiphase solid materials containing at least one phase with at least one dimension ranging from a few angstroms to a hundred nanometers, have attracted tremendous attention over the past few decades and have been recognized as one of the most promising materials systems for producing multifunctional properties [1,2]. A variety of nanocomposites containing or consisting of globular nanoparticles [3,4], nanocolumns [5,6], and layered structures [4,7] have already been fabricated and explored in many fields of application [2,5]. Chemico-physical properties of nanocomposites can be greatly influenced by nanoparticles size, shape, and mutual interactions [8], but the determination of two or more known substances and the preparation of nanostructures composed of those substances to optimally suit their applications, which usually encompasses the discovery of the optimal method and growth conditions, remain the challenging tasks. With the rapid development of nanoscience, encapsulation of nanoparticles has become an important issue. Encapsulation of nanoparticles in nanoshells or nanopores is expected to enable the physical isolation of nanoparticles, which could inhibit their migration and coalescence, making encapsulation a promising strategy to overcome various stability issues [9]. Various inorganic oxides, carbon, and organic cages can be used as encapsulating materials [9]. A better understanding of the mechanisms leading to the formation of a particular structure can help to obtain a better nanocomposite for a given application. A better understanding of nanocomposite growth mechanisms is

Materials and Methods
We employed a three-dimensional grid to simulate the growth processes of a thin film in three spatial dimensions and time. The composition of a growing thin film in a modeling domain is described by using the relative concentrations of thin film components. associated with each grid point, and they represented the concentrations of components A, B, and substrate material, respectively, at grid point i, j, k. The indices i, j specified the positions of a respective grid cell in the horizontal directions, k indicating the position in the vertical direction that was the film growth direction. The condition c i,j,k A + c i,j,k B + c i,j,k S ≤ 1 was always met in any grid cell. A modification of previously published models [20] was employed in this work to simulate the growth processes of thin films. The modified model still includes the processes of phase separation, diffusion of both thin components due to surface curvature, and adsorption. The mathematical description of the phase separation of both components occurring through the surface diffusion was based on the Cahn-Hilliard equation [24] and on the idea presented in our previous work [20]. The present model uses the same equation and variables as those used in [20] to obtain a mathematical description for the phase separation; however, the present model is constructed slightly differently. Changes in the concentrations of component A due to phase separation in the surface layer of a growing film were described by the following equation [20]: where c AS (c i,j AS , i = 1, . . . , I, j = 1, . . . , J) denotes the concentrations of component A in the film surface cells, D A is the diffusion coefficient, f (ϕ) = 2 2ϕ 10 − ϕ 8 − ϕ 6 − ϕ 4 − ϕ 2 + 2 is the function of the free-energy density of a homogenous system (the even function f (ϕ) used in this work is same as that used in [25], which is related to the solubilities of components in a binary system), and γ is the power coefficient of the phase gradient. The concentrations of component A in the surface layer, c i,j AS (this is the concentration variable used in Equation (1)), is defined through by following relation [13,20]: where D AA , D AB , and D AS are the diffusion coefficients of component A on pure phase A, phase B, and substrate material, respectively. The weights c AF , c BF , c SF are the surface relative concentrations (they are assessed considering the one full surface layer) of component A, component B, and substrate material, respectively, and they are defined for each grid surface cell separately by the following relationships [20]: From the relationships given in Equation (4) the relative concentration of component AF , is defined as the sum of the concentrations of component A from the respective highest occupied grid cell denoted by i, j, k * and the one grid cell below it indicated by i, j, k * − 1.
The process of phase separation described by Equation (1) (2) and (3) by adjusting them to be suitable for component B: where k * is the highest position (in the growth direction) at which c i,j,k  (4).
The diffusion of atoms on nanoislands due to their surface curvature is another process included in the model [26]. Changes in the concentrations of component A due to the surface curvature of a growing film were defined by using the following equation [20] (the mathematical description of this process is the same as that in [20]): The term ∇ 2 h denotes that the surface curvature, h, was the position of the thin film surface, and p A was the proportionality coefficient. Equation (6) ensures that atoms of respective type diffuse from crests to valleys of nanoislands due to the surface curvature. Equation (6) also defines changes in the concentrations of component B due to the surface curvature, but c BS , p B , and D B have to be used instead of c AS , p A , and D A , respectively, in Equation (6).
Changes in the concentrations of both depositing species (at the highest occupied ) caused by the adsorption process were described by using the following equations [20] (the mathematical description of this process is same as in [20]): where i = 1,..., I, j = 1,..., J, k AA , k AB , and k AS are the sticking coefficients (k AB denotes the sticking coefficient of A-type atoms to the surface made of component B); and i A and i B are the relative fluxes of both film components; c i,j AF , c i,j BF , and c i,j SF are the surface concentrations of components A, B, and substrate material, respectively. The sum, ∂c /∂t, amounted to the total growth rate of the thin film.
All the three equations describing the three different processes with respect to the film component were added together to obtain the final model for that component. Because we investigated two growing phases, the two separate models dealing with each growing phase and consisting of Equations (1), (6) and (7) were employed. Our previous model [20] used the two Cahn-Hilliard equations for each component, so it consisted of the four summands.

Results and Discussion
To investigate and determine factors resulting in the growth of branching structures, a series numerical experiments was performed. The lateral dimensions of the simulation domain were 36 nm × 36 nm, and the spacing of the computational grid was 1 nm. The boundaries of the simulation domain parallel to the growth direction were used as the periodic boundary conditions. In the simulation domain, the substrate was represented as five perfectly flat layers occupying those layers counting from the bottom of the simulation domain. The basic values of the parameters used in calculations are given in Table 1. The parameters D AA , D BB , D AB , D BA , k BA , k BB , i A , i B , p A , p B , and γ were kept constant in all calculations. Some of those parameter values are consistent with those used in our previous works [13,20,27]. As shown by [27], the diffusion values similar to those of D AB and D BA given in Table 1 resulted in the formation of columnar structures throughout the whole film thickness (diffusion coefficients D AB and D BA determine the phase separation rate in the later growth stages when the substrate is covered with depositing species). The presence of nanocolumns was also observed in the branching structures being considered. The diffusion coefficients D AA and D BB were chosen of the same order as D AB , D BA . The values of p A and p B used in this work were lower than those used in our previous works [13,20]. Those parameters determine the diffusion on nanoislands due to their surface curvature. The influence of that effect was not investigated in this work, so the relatively smaller values of p A and p B were used. The parameters D AS , D BS , k AA , k AB , k AS , and k BS were varied in respective numerical experiments. Any changes in their values will be specified.   Figure 1a, those boundaries only visually divide all the clusters formed on the substrate. The relationship between lateral dimensions of nanoparticles made of some component and diffusion coefficient of that component was also observed in other works [13,20,28,29]. From the images in Figure 1a-d it can be seen that the relatively large nanoclusters made of component A (with respect to the lateral dimensions) formed during the initial stages of the growth do not keep their lateral sizes throughout the whole film thickness, and the large clusters of component A formed on the substrate material start to branch into smaller ones. In the vertical cross section in Figure 1d, we observe relatively narrow elongated (in the growth direction) nanoparticles made of component A grown from a basic relatively large nanocluster formed in the initial stages of film growth. As can be seen from Figure 1d, the smaller nanoparticles of pure phase A remained connected to the large cluster made of phase A. As the image in Figure 1d suggests, such a branching process started to occur right after the majority part of the substrate had been covered with either phase A or B, so this branching process must be related to some different growth mechanisms resulting in the growth of different structures in the early and later growth stages. Judging from the values of the parameters used in the calculations presented in Figure 1, besides the relatively high value of D AS (the diffusion coefficient of component A on the substrate material is two orders of magnitude higher than the diffusion coefficients of component A on pure phases of A and B) used to generate the images in Figure 1a Figure 1e-h, the four-times lower diffusion coefficient, D AS , reduces the dimensions of nanoparticles of component A formed on the substrate material by two to three times. However, as can be seen from the images in Figure 1g,h, the branching process still occurred. Figure 1h suggests that smaller lateral dimensions of initial clusters of component A formed on the substrate result in a smaller number of elongated nanoparticles of component A extending from those initial clusters to the surface of the grown film. Figure 1i shows that the further decrease in the diffusion coefficient, D AS , to 8.8 × 10 −19 m 2 /s changed the structure formed on the substrate very significantly. As can be seen in Figure 1i, there are no encapsulated nanoparticles of component A in the given horizontal cross section through the plane z = 6 (with respect to the horizontal directions); the clusters made of component A are connected to one another. This behavior is caused by a relatively high concentration of phase A grown on the substrate material and by a relatively low value of the diffusion coefficient of component A on the substrate material, D AS (D AS = 8.8 × 10 −19 m 2 /s used to calculate Figure 1i-l was lower than the value of D BS , which was kept constant in the presented calculations). As shown by the images in Figure 1j,k, although the composition of the thin film obtained with the lowest diffusion coefficient D AS in the mid-depth and surface layers remain similar to the cases obtained by using the higher diffusion coefficient D AS , the phase structures of layers grown on the substrate are quite different. As shown by the vertical cross section given in Figure 1l, no distinguishable branching structures can be observed at the relatively low value of D AS . given in Figure 1l, no distinguishable branching structures can be observed at the relatively low value of DAS.  As previously mentioned, the formation of the branching structures observed in Figure 1a-h was not only influenced by the diffusivities of deposited components on the substrate material but also by the values of the sticking coefficients. Figure 2 shows horizontal and vertical cross sections obtained by using different sticking coefficients of components A and B to the substrate material. All the other model parameters were the same as those used to calculate the structure in Figure 1a Comparing the images in Figures 1a and 2a,e, we notice that such a change in the value of the sticking coefficient of component, A, to the substrate , kAS (keeping the sum kAS + kBS constant), clearly has an influence on the phase structures formed on the substrate. In all the three cases shown in Figures 1a and 2a,e, we observe groups of clusters of com-  Comparing the images in Figures 1a and 2a,e, we notice that such a change in the value of the sticking coefficient of component, A, to the substrate, k AS (keeping the sum k AS + k BS constant), clearly has an influence on the phase structures formed on the substrate. In all the three cases shown in Figures 1a and 2a Figure 2e). Therefore, the average lateral sizes of clusters of component A formed on the substrate (and calculated in the plane z = 6) decrease from 17.8 nm to 10 nm as the value of the sticking coefficient, k AS , is lowered from 0.75 to 0.45 (keeping the sum k AS + k BS constant). As can be seen in Figure 1b,c and Figure 2b,c,f,g, the changes in the sticking coefficients related to the substrate do not significantly influence the growth dynamics and phase structures obtained in the later growth stages. In the vertical cross sections in Figures 1d and 2d,h, we notice that the branching process of the relatively large clusters of component A (formed on the substrate in the initial growth stages) occurred in all three cases presented. The only difference between those cases is the number of elongated nanoparticles of component A connected to the large basic clusters and extending from those clusters to the surfaces of the grown films. If the lateral dimensions of any basic cluster of component A are smaller, then fewer elongated nanoparticles can grow from that basic large cluster because, as can be seen in Figures 1c and 2c,g, the lateral dimensions of nanoparticles of component A in the plane z = 18 (the layer near the surface of the film) did not change significantly as the value of the sticking coefficient k AS was varied.
We can conclude that the branching mechanism observed in the cases in Figures 1 and 2 consists of two main steps: the first is the formation of a relatively large cluster of a component that makes up the branching phase. The second step is the division process of that large cluster into smaller clusters (with respect to their lateral dimensions) as a film thin grows thicker, as well as the formation of nanocolumns, elongated in the growth direction, extending from the initial large cluster to the surface of the grown film. We can presume that the first step should be achieved by any process or effect of interaction between the phase (branching phase) and the substrate material, which can result in relatively large lateral dimensions of clusters formed on the substrate. However, that effect must disappear in the later growth stages when the substrate is covered with a deposited film, and the growth conditions in the later stages must be favorable to the formation of nanoparticles with a smaller lateral size to achieve the branching process. The important role of the substrate material is also supported by the fact that the formation of such branching structures in yttrium-stabilized zirconia (YSZ) thin films grown by using e-beam deposition technique [22] was only observed when the film was deposited on a porous Ni-YSZ substrate.
As previously revealed, the formation of relatively large nanoclusters on the substrate during early growth stages is the first necessary condition for the branching process to occur through the mechanism described. Figure 3 Figure 1). According to Figure 3, the average lateral grain size near the substrate monotonically increases with the increase in the diffusion coefficient, D AS . The relationship between the average lateral grain size and the diffusion coefficient given in Figure 3 can be approximated by a logarithmic trendline, which, as shown by the value of R 2 in Figure 3 the substrate material in the early growth stages. The further increase in the diffusion coefficient D AS from 8.8 × 10 −18 m 2 /s to 8.8 × 10 −17 m 2 /s resulted in an increase in the average cluster size (of component A) from 9.5 to 18 nm; the average cluster sizes in that range were favorable for the branching process to occur. mated by a logarithmic trendline, which, as shown by the value of R 2 in Figure 3, explains 95.6 of the variation of the average lateral grain size related to the changes in the diffusion coefficient, DAS, in the given range. An analysis of the structures obtained with all five diffusion coefficient DAS reveals that, only at the lowest value of DAS (at 8.8 × 10 −19 m 2 /s), no branching process was observed. At the lowest value of DAS used, the formation of clusters with a size of only 5.2 nm, on average, was observed on the substrate material in the early growth stages. The further increase in the diffusion coefficient DAS from 8.8 × 10 −18 m 2 /s to 8.8 × 10 −17 m 2 /s resulted in an increase in the average cluster size (of component A) from 9.5 to 18 nm; the average cluster sizes in that range were favorable for the branching process to occur.   Figure 4 shows the plot of average lateral cluster size (of component A) near the substrate versus difference k AS − k BS . The condition k AS + k BS = 0.9 was maintained when the difference between k AS and k BS was varied. The average lateral grain sizes were assessed at the following values of the differences between k AS and k BS : −0.6, −0.3, 0, 0.3, and 0.6. Therefore, the difference k AS − k BS = −0.6 corresponds to the pair of k AS = 0.15, k BS = 0.75, k AS − k BS = 0 corresponds to k AS = 0.45 and k BS = 0.45, k AS − k BS = 0.6 corresponds to k AS = 0.75 and k BS = 0.15, etc. The other model parameters were the same as those used to calculate the structure given in Figure 1a-d. Figure 4 shows a monotonic relationship between the average grain size near the substrate and the difference k AS − k BS ; such a relationship confirms that the lateral sizes of nanoclusters near the substrate can be influenced by the values of the sticking coefficients. According to Figure 4, the increase in the difference k AS − k BS from −0.6 to 0 increases the average lateral size of nanoclusters from 8.6 to 10 nm, but the further increase in the difference k AS − k BS from 0 to 0.6 results in a significantly greater increase in the average cluster size from 10 to 18 nm. Such an asymmetric relationship should be highly conditioned by the use of different relative fluxes of depositing species (i A = 0.2666 s −1 , i B = 0.5333 s −1 ) and different diffusion coefficients on the substrate material, D AS and D BS . From an analysis of the structures obtained with all five differences k AS − k BS considered, we noticed that the branching process was observed in all the cases, but an increase in the average lateral cluster size near the substrate caused changes in the patterns of the branching structures obtained. A larger lateral cluster size near the substrate usually had more elongated nanocolumns extending from that cluster to the surface of the film.  Figure 5a-d,i-l the sticking coefficients, k AA and k AB , were reduced monotonically. In Figure 5a,c,i, where the concentration plots of component A near the substrate are presented, we notice that the changes in the sticking coefficients, k AA and k AB , do not have much of an influence on the structures formed on the substrate material in the initial growth stages. As shown by Figure 5a,e,i, the laterally large clusters of phase A are formed on the substrate in all three cases given, which is not an unexpected result because the sticking coefficients, k AA and k AB , influence the adsorption rates of type-A particles on only pure phase A and pure phase B (and not on the substrate material), respectively, and those coefficients should have a relatively low influence on the structures formed on the substrate. The effect of the sticking coefficients, k AA and k AB , on the phase structures formed is noticeable in the images given in Figure 5b,c,f,g,j,k, where the distributions of component A in the mid-depth and surface layers of the films, respectively, are provided.
As can be seen in those images, the fractions of component A in the given layers visually decrease as the values of the two sticking coefficients are lowered from 0.75 to 0.5 and from 0.5 to 0.25, respectively. As can be seen in Figure 5b,c,f,g,j,k, the lateral dimensions and aspect ratios of nanoparticles made of component A (with respect to the planes considered) are highly influenced by the values of the sticking coefficients, k AA and k AB . In Figure 5j,k, we notice that at the lowest sticking coefficients used, the concentrations of component A in the layers considered are relatively low, and the aspect ratios of nanoparticles made of component A with respect to the planes considered are near the ratio 1:1 for almost all nanoparticles. In Figure 5b,c,f,g we notice that at the higher values of the sticking coefficients used, the dimensions and aspect ratios of nanoparticles of component A with respect to the planes considered are significantly different than in the case with the lowest values of the sticking coefficients. In the case with the highest two values of both sticking coefficients, k AA and k AB , we observe an increase in the lateral dimensions of nanoparticles (with respect to planes considered) of component A, and those increases in the lateral dimensions of nanoparticles are mostly seen in the only one direction, which results in higher deviations of the aspect ratios of nanoparticles made of component A from the ratio 1:1. As can be seen in Figure 5a,e,i, the formation of relatively large (in their lateral dimensions) clusters made of component A on the substrate in the initial growth stages should be favorable to the development of branching structures in all the cases presented; however, as can be seen from Figure 5, the branching structures with elongated nanoparticles made of component A, extending from basic large clusters of component A lying on the substrate to the surface of the grown film, are observed in the cases with lower values of sticking coefficients k AA and k AB . Therefore, a relatively low amount of component A being adsorbed in the later growth stages when the substrate is fully covered with the depositing species is an important condition for the branching process to occur. If the amount of component A adsorbed in the later growth stages is too high, we do not observe any elongated nanoparticles of component A extending from the substrate to the surface of the grown film, which would be an essential feature of the branching structures being considered; however, in those cases, during the later growth stages, we observe the formation of phase structures containing some higher-dimension nanoparticles of component A. The fraction of particles of component A adsorbed in the later growth stages should not be too low because a lower amount of component A would result in a lower number of branching nanocolumns, which is also an important feature of branching structures. s −1 , iB = 0.5333 s −1 ) and different diffusion coefficients on the substrate material, DAS and DBS. From an analysis of the structures obtained with all five differences kAS − kBS considered, we noticed that the branching process was observed in all the cases, but an increase in the average lateral cluster size near the substrate caused changes in the patterns of the branching structures obtained. A larger lateral cluster size near the substrate usually had more elongated nanocolumns extending from that cluster to the surface of the film.   Figure 5a,c,i, where the concentration plots of component A near the substrate are presented, we notice that the changes in the sticking coefficients, kAA and kAB, do not have much of an influence on the structures formed on the substrate material in the initial growth stages. As shown by Figure 5a,e,i, the laterally large clusters of phase A are formed on the substrate in all three cases given, which is not an unexpected result because the sticking coefficients, kAA and kAB, influence the adsorption rates of type-A particles on only pure phase A and pure phase B (and not on the substrate material), respectively, and those coefficients should have a relatively low influence on the structures formed on the substrate. The effect of the sticking coefficients, kAA and kAB, on the phase structures formed is noticeable in the images given in Figure 5b,c,f,g,j,k, where the distributions of component A in the mid-depth and surface layers of the films, , respectively, are provided.  As can be seen in those images, the fractions of component A in the given layers visually decrease as the values of the two sticking coefficients are lowered from 0.75 to 0.5 and from 0.5 to 0.25, respectively. As can be seen in Figure 5b,c,f,g,j,k, the lateral dimensions and aspect ratios of nanoparticles made of component A (with respect to the planes considered) are highly influenced by the values of the sticking coefficients, kAA and kAB. In Figure 5j,k, we notice that at the lowest sticking coefficients used, the concentrations of component A in the layers considered are relatively low, and the aspect ratios of nanoparticles made of component A with respect to the planes considered are near the ratio 1:1 for almost all nanoparticles. In Figure 5b,c,f,g we notice that at the higher values of the sticking coefficients used, the dimensions and aspect ratios of nanoparticles of component A with respect to the planes considered are significantly different than in the case with the lowest values of the sticking coefficients. In the case with the highest two values of both sticking coefficients, kAA and kAB, we observe an increase in the lateral dimensions of nanoparticles (with respect to planes considered) of component A, and those increases in To provide a better representation of the nanostructures modeled, Figure 6 Figure 1; the images given in Figures 1 and 6 were drawn from the same modeling data. Component B and the substrate were removed from the given plots to obtain a better visual representation of nanoclusters made of component A, as we are focused mainly on the form of those clusters. As can be seen in the three-dimensional plots given in Figure 6, some areas near the substrates (the top of the substate is always located at z = 5) are obscured by other parts of nanoclusters, so those three-dimensional plots should be examined together with the respective cross-sectional images given in Figure 1. Figure 6 confirms that the growth of branching structures made of component A was observed only at the two highest diffusion coefficient D AS (8.8 × 10 −17 m 2 /s (a) and 2.2 × 10 −17 m 2 /s), whereas at the lowest value of D AS , the growth of branching structures was inhibited.
diffusion coefficients of component A on the substrate material, DAS. The values of DAS were as follows: DAS = 8.8 × 10 −17 m 2 /s (a), DAS = 2.2 × 10 −17 m 2 /s (b), and DAS = 8.8 × 10 −19 m 2 /s (c). The cross-sectional views obtained with those diffusion coefficients are presented in Figure 1; the images given in Figures 1 and 6 were drawn from the same modeling data. Component B and the substrate were removed from the given plots to obtain a better visual representation of nanoclusters made of component A, as we are focused mainly on the form of those clusters. As can be seen in the three-dimensional plots given in Figure 6, some areas near the substrates (the top of the substate is always located at z = 5) are obscured by other parts of nanoclusters, so those three-dimensional plots should be examined together with the respective cross-sectional images given in Figure 1. Figure 6 confirms that the growth of branching structures made of component A was observed only at the two highest diffusion coefficient DAS (8.8 × 10 −17 m 2 /s (a) and 2.2 × 10 −17 m 2 /s), whereas at the lowest value of DAS, the growth of branching structures was inhibited. To gain a better understanding of the mechanism and growth of branching structures in time, Figure 7 provides three-dimensional plots of nanoclusters made of component A observed at different growth times. The plots are represent the following growth times: 0.1T (a), 0.2T (b), 0.3T (c), and 0.5T (d), where T is the total growth time modeled. The plots given in Figure 1e-h, 6b, and 7 were taken from the same numerical experiment, but they show the structures observed at different growth times. Figure 6b represents the structure observed at the growth time, T. Figure 7a shows a group of nanoislands of component A formed on the substrate material at the growth time 0.1T. The heights of those nanoislands Such a process of nucleation can only occur when the formation of relatively large clusters is favorable in the initial growth stages when nanoislands grow on the substrate, whereas the formation of nanoparticles with smaller lateral dimensions is favorable in the later growth stages. Such a change in the lateral dimensions of nanoclusters is influenced by the differences between the diffusivities (and sticking coefficients) of depositing species on the substrate material and pure their phases. As shown in Figure 7c, considering the relatively wide nanoislands of component A formed on the substrate, phase A continues growing primarily on outer parts of those nanoislands (although, as shown by Figure 6a, phase A may continue growing not only on outer parts of nanoislands but also on inner parts if the lateral dimensions of a nanoisland are large enough). In those areas, the growth of phase A can be still favorable because any particle of type B adsorbed to outer areas of a nanoisland of component A can easily diffuse to reach clusters made of component B, which are formed all around any nanoisland of component A or are nucleated on the inner areas of any sufficiently large nanoisland of component A. In Figure 7c, we notice that on some outer parts of the nanoislands of component A (where the growth of phase A was favorable), phase A grows slower than on the others, which, as shown by Figure 7d, leads to the division of the nanoclusters made of phase A into smaller cluster. Those smaller cluster continue growing as nanocolumns extending from the large initial clusters to the surface of the grown film. Such a branching process is promoted (in the later growth stages) by decreasing the amount of component A adsorbed in the later growth stages. In our model, the amount of component A adsorbed in the later growth stages is simply determined by the sticking coefficients k AA and k AB (the sticking coefficients of component A to pure phase A and B, respectively). The influence of substrate temperature on the branching process can be discussed on the basis of this study, as well as our previous works [13,20]. If we consider the temperature range 0.2 < ( is melting temperature), diffusion coefficients increase exponentially with the substrate temperature, TS, in that range. According to the Langmurian temperature dependence of adsorption processes, the growth rate is proportional to The influence of substrate temperature on the branching process can be discussed on the basis of this study, as well as our previous works [13,20]. If we consider the temperature range 0.2T m < T (T m is melting temperature), diffusion coefficients increase exponentially with the substrate temperature, T S , in that range. According to the Langmurian temperature dependence of adsorption processes, the growth rate is proportional to 1/ √ T s . Therefore, diffusion coefficients increase with an increase in the substrate temperature, and the growth rate decreases with an increase in the substrate temperature. In [13], it is shown that the average lateral size of nanoclusters increases with an increase in the substrate temperature. In [13], it is also shown that the average lateral size of nanoclusters increases with a decrease in the growth rate. Therefore, based on the findings of this and other works [13,20], a relatively high substrate temperature (assuming both dependencies of diffusivities and the growth rate on substrate temperature specified) should be favorable for the formation of laterally large clusters in initial growth stages, but high substrate temperatures are not favorable for the formation of nanocolumns with relatively small lateral dimensions. This leads to the conclusion that relatively high substrate temperatures in initial growth stages and lower substrate temperatures in later growth stages should be favorable for the formation of branching structures.

Conclusions
(1). The modeling results revealed the possible mechanism through which the formation of branching structures with nanocolumns elongated in the growth direction and extending from a relatively large cluster formed on the substrate to the surface of a grown film could be achieved. (2). The growth mechanism of branching structures can be divided into the two main steps. The first step is the growth of a relatively large cluster (of a component that makes up the branching phase) on the substrate. The second step is the division process of that large cluster into smaller clusters as a film thin grows thicker, as well as the formation of nanocolumns elongated in the growth direction and extending from the initial large cluster to the surface of the grown film.