Abstract
This paper presents recent advances in the theory of multicomponent crystal growth from gas or solution, focusing on the most common step-flow mechanisms: Burton-Cabrera-Frank, Chernov, and Gilmer-Ghez-Cabrera. Analytical expressions for the spiral crystal growth rate are presented, taking into account the properties of all species involved in the growth process. The paper also outlines theoretical approaches to consider these mechanisms in multicomponent systems, providing a foundation for future developments and exploration of previously unexplored effects. Some special cases are discussed, including the formation of nanoislands of pure components on the surface and their self-organization, the impact of applied mechanical stress on the growth rate, and the mechanisms of its influence on growth kinetics. The growth due to chemical reactions on the surface is also considered. Possible future directions for developing the theory are outlined. A brief overview of numerical approaches and software codes that are useful in theoretical studies of crystal growth is also given.
1 Introduction
In the past century, the field of crystal growth has emerged as a distinct and important branch of materials science, driven by the increasing demand for crystals in various technological applications such as electronics, optics, semiconductors, and more. While silicon formed the basis of modern electronics in the 20th century, high-tech industries are gradually transitioning towards more complex crystalline materials that can facilitate the creation of devices with superior characteristics. These materials include widebandgap semiconductors such as gallium and aluminum nitrides (GaN, AlN), which are widely used in optoelectronics (), perovskite solar cells (; ), prospective superconductors (), thermoelectric materials (), and other functional compounds (; ), most of which consist of multiple components. A rapidly growing area of research is the development of metal-organic and other types of frameworks (; ) for gas capture, separation, storage (), catalysis (), and various other applications. These materials are crystalline-like in structure and their growth mechanisms are similar to those of crystals, albeit more complex (See Figure 1B). They exhibit the same growth phenomena and in certain cases can be described by the same models and equations.
FIGURE 1
Note that there are various regimes of crystal growth, such as step-flow/spiral, nucleation, dendritic, and rough growth (refer to Figure 1A), and the realization of each depends on crystal properties and growth conditions. The growth regime determines the growth rate, crystal purity, defect concentration, and homogeneity, and a deep understanding of all the processes that occur during growth in different regimes is crucial for controllable synthesis of crystals with desired properties. This necessity has led to the rapid development of crystal growth theory and the detailed study of the growth regimes, especially step-flow/spiral and nucleation, since they are mostly used in the industry. The classical work by Burton, Cabrera, and Frank (BCF) () is a cornerstone of the theory that explains step-flow growth from vapors at low supersaturations due to screw dislocations. It revolutionized the understanding of crystal growth and allowed for the quantification of theoretical models and experiments. Since then, the theory has progressed, and various phenomena and effects have been investigated, such as the influence of elastic stresses (), advacancies (; ; ), and impurities (), among others. Analytical descriptions of different growth regimes and instabilities () have been proposed, including the formation of dendrites (), the impact of island formation on terraces (), various types of self-organization, and other factors. There are numerous reviews and books which summarize the results (; ; ) and different aspects of the theory, evolution of steps () and kinks () on the surface. The BCF model has stimulated the development of other important growth models such as the Chernov model for growth from solutions or gas phase in the presence of a carrier gas () and the mixed theory of Gilmer-Ghez-Cabrera (). Together with nucleation theory (), these models have facilitated the quantitative description of most gas or solution growth processes used in industry and accelerated the development of technologies and methods for obtaining low-defect crystalline materials like silicon. However, the majority of theories were developed using the classical model of the Kossel crystal. While it works well in single-component cases, it cannot always describe the growth of multicomponent crystals due to the different properties of the atoms (diffusion coefficients, atomic volumes, etc.). This calls for the extension of existing growth models and theory to describe crystal growth in multicomponent systems, which is in high demand in modern technology.
This perspective aims to discuss our recent developments in the theory of crystal growth in multicomponent systems via step-flow/spiral mechanisms, specifically focusing on the BCF (), Chernov (), and Gilmer-Ghez-Cabrera () models. We note that the multicomponent growth via the second important mechanism - nucleation, both classical and non-classical, is summarized elsewhere (; ). The article is organized as follows: In Section 2, a description of the mechanisms and their applicability limits is provided to help readers better understand which model is best suited for their experiments. Section 3 covers various aspects and effects that arise during growth via the step-flow mechanism. In Section 4, a discussion of growth mechanisms and effects is presented. The final expressions for the crystal growth rates are summarized in Supplementary Table S1. For the convenience of the readers, classical single-component expressions are also included. Additionally, the article provides a rigorous mathematical formulation of the models and a brief course on their solution in the Supplementary Materials, allowing readers to modify the models for their specific purposes if necessary. It should be noted that examples of the application of the proposed models to describe the growth of specific materials are provided in our referenced papers.
2 Crystal growth mechanisms
The problem of multicomponent spiral growth was formulated in () based on the classical BCF model. In this formulation, the crystal growth process involves the incorporation of building units into kinks at the spiral steps (see Figure 1C) on the surface in accordance with the following reaction:where represents the pure components (or building units), is the solid crystalline phase; is the stoichiometric coefficient. The subscripts g) and s) denote the gas phase and solution, respectively, and N is the number of components. As previously mentioned, the mass transfer of components from the media to the kinks on the steps can occur through various mechanisms (as illustrated in Figure 1D), which will be discussed in detail later.
2.1 Burton-Cabrera-Frank mechanism
The Burton-Cabrera-Frank (BCF) model is commonly used to describe step-flow and spiral crystal growth from vapor. In the multicomponent version of this model (), all components that constitute the crystal are delivered to the terraces via uniform gas fluxes over the entire crystal surface. These components adhere to the surface and become adatoms (or admolecules), which then diffuse over the terraces towards existing steps and incorporate into them, thereby ensuring step motion and crystal growth. Adatoms can also evaporate back to the gas phase. Each type of adatom on the surface has unique properties, such as the surface diffusion coefficient and lifetime before evaporation. The strict mathematical formulation of this and the following models is presented in the Supplementary Materials. Here, we do not focus on mathematics, which is presented in detail in (), and instead provide only the resulting analytical expression for the growth rate R (Supplementary Table S1). It turns out that the process of growth of a multicomponent crystal may be described by the same equation as a single-component crystal, using the so-called generalized coefficients of diffusion and incorporation , and generalized supersaturation . These coefficients are determined by the properties of the individual components and their surface concentration ratios. It can be seen from the definition of the generalized coefficients that the growth process is typically limited by the component that has the smallest product of diffusion or incorporation coefficient and concentration. It should be noted that this model assumes that: 1) reaction 1) occurs only and directly at the kinks, and not in the gas phase or on the terrace, which is valid for many growth processes; 2) the surface coverage is low, and the adatoms do not affect the diffusion or deposition of other atoms; 3) the step is always fully covered with kinks and is therefore a good sink for adatoms; note that the mechanisms of incorporation into kinks are presented in detail in the review (); and 4) the rate of advancement of the step is sufficiently slow to neglect the convective flux of the adatoms towards the step.
2.2 Chernov mechanism
The Chernov model () is commonly used to describe crystal growth from solution, but it is also applicable for growth from intrinsic vapors highly diluted by the carrier gas. Compared with the BCF model, the Chernov model has two key differences. First, the mass transfer in the volume is insufficient to maintain constant pressures/concentrations directly at the crystal surface, resulting in a depleted diffusion layer of some thickness δ in the mother liquid (or gas) phase. Second, incorporation occurs directly from the volume into the steps, bypassing the intermediate state on the terrace; thus, there is no surface diffusion. In the multicomponent version of this model, the crystal still grows in accordance with Reaction 1). At one of the boundaries of the depleted diffusion layer, constant concentrations of the components are maintained due to mixing or stirring (see Figure 1D). The components then diffuse through this layer and are incorporated directly into the steps one after another, resulting in crystal growth. The assumptions behind this model are almost the same as those in the BCF model (see section 2.1). However, it is important to note that this model fully neglects surface diffusion and is therefore only valid for systems with slow surface kinetics (e.g., owing to high surface diffusion barriers, surface reconstruction, coverage of the surface by an arbitrary adsorbate layer interfering with the surface diffusion of components, or other reasons). Furthermore, the thickness of the boundary layer is typically of the order of 10−3 to 10−5 cm, which is much smaller than the typical size of the growing crystal. As a result, some macro-diffusional fields may arise during growth, which must be considered to correctly describe the crystal growth over the entire surface.
2.3 Gilmer-Ghez-Cabrera mechanism
The Gilmer-Ghez-Cabrera model () is a significant generalization of the Burton-Cabrera-Frank and Chernov models. While the BCF and Chernov models consider only one type of diffusion, the Gilmer-Ghez-Cabrera mechanism takes into account both surface diffusion and volume diffusion of substances to the steps in the presence of a solvent or carrier gas. This is particularly important for multicomponent systems, where the flux of some components to the steps may be limited by volume diffusion, whereas for others, it may be limited by surface diffusion. The BCF and Chernov models may yield incorrect results in such cases because they consider only one type of flux (refer to Figure 1D). Notably, in the Gilmer-Ghez-Cabrera mechanism, adatom incorporation occurs through an intermediate adsorbed state on the terrace because there is no direct incorporation, as in the Chernov mechanism. It is worth mentioning that the Van Der Eerden model () is the most complex model that accounts for both direct and indirect incorporation. However, this model has not yet been extended to multicomponent systems.
3 Different effects inherent to multicomponent systems
3.1 Nucleation of pure components on the crystal surface
In (), the authors demonstrated that nucleation of islands consisting of pure components on the terraces between the steps is possible during step-flow growth, similar to the process of dew formation (refer to Figure 2A). In a subsequent study (), it was shown that the precipitation of such islands can significantly impact the growth process, resulting in the formation of undesirable defects and inclusions in the crystal, and even cause different types of morphological instability (). This phenomenon may occur when the surface concentration of a particular component on the surface between the steps exceeds the equilibrium single-component concentration required for nucleation. The authors also developed a dynamical criterion for determining whether nucleation of the ith component will occur during the BCF growth mode in a train of advancing steps.where represents the time required for the formation of a single island at the surface concentration maintained at the terrace between subsequent steps, and , is the time required for the multicomponent step to traverse the terrace width and “sweep” it from adatoms, allowing for the nucleation process to restart.
FIGURE 2
3.2 Impact of chemical reactions
The growth of crystals in the presence of chemical reactions is common in modern synthesis techniques. For example, the metal-organic chemical vapor deposition (MOCVD) process used for GaN crystal growth typically involves surface chemistry (). In (), the authors investigated BCF-type growth in the most general case:when both the initial reagents and reaction products , are present and do not correspond to the building units of the crystal (see Figure 2B). It was shown, that in this complex case, the reaction must be broken down into sub-reactions, each bringing the building units/adatoms of the crystal to the surface. The growth in this case can still be described by the same single-component formula (see Supplementary Table S1), but if the thermal desorption of adatoms is slow and can be neglected in comparison to the fast reaction rate, the mean free paths and lifetimes of adatoms on the surface are primarily determined by reverse chemical sub-reactions, that is, by the pressures of the reaction products. Simultaneously, supersaturation will depend on the chemical reaction affinity, , where is the equilibrium constant of the reaction, and . This provides additional control levers for the crystal growth process.
3.3 Impact of mechanical stress and advacancies
Previous studies have established that mechanical stress may have a significant effect on crystal growth (). In (), it was demonstrated from a thermodynamic point of view that stress affects the growth of multicomponent crystals in the same way by reducing the effective supersaturation according to: , where is the isotropic elastic stress and is the volume of the crystalline cell. This reduction factor can be crucial for the epitaxial growth of semiconductor or MOF-on-MOF heterostructures. Another important factor is the emergence and diffusion of vacancies on the crystal surface, which can recombine with adatoms and significantly affect the crystal growth or evaporation (; ; ; ). In (; ), the authors extended the BCF model to the multicomponent case and accounted for the formation of additional vacancies resulting from applied stress (see Figure 2C). In this case the characteristic diffusion length, which determines the growth rate, is limited by the lowest diffusion lengths of adatoms and advacancies due to evaporation and recombination (see Figure 2C). It was demonstrated that the processes of adatom deposition, diffusion, and incorporation into the steps due to supersaturation in the gas phase are similar to the same processes for advacancies but due to applied tensile loads.
3.4 Morphological instability in multicomponent systems
The phenomenon of morphological instability was first introduced by Mullins and Sekerka (). It typically results in a change in the shape of a growing crystal or thin film, leading to irregular dendritic shapes or other inhomogeneities at a specific spatial wavelengths. In (; ) the authors expanded the Mullins-Sekerka morphological stability theory to the growth of multicomponent strained thin films and spherical multicomponent particles. Their research showed that once a certain level of supersaturation is exceeded, the aforementioned shapes become unstable, leading to the emergence of irregular dendritic or undulating shapes. Analytical criteria were developed for this phenomenon, which connect the properties of different components (e.g., diffusion coefficients and atomic volumes), ratio of concentrations, and total supersaturation. This allows the selection of stable growth conditions for the crystal surface.
4 Discussion
The models presented here allow for the description of a wide range of growth processes from gases or solutions. They also provide a means for the experimental determination of all the kinetic coefficients for each component in complex multicomponent systems. As demonstrated in (), the reaction 1) between adatoms at the kinks places a constraint only on the product of equilibrium concentrations which should equal to the equilibrium constant ; , and on the stoichiometry of the fluxes of the components from terrace to kink and back, but it does not determine themselves. This is in contrast to single-component systems, in which the equilibrium concentration is fixed at a constant temperature. Therefore, there are many possible pressure (or concentration) equilibrium sets at which a multicomponent system can exist. By selecting a set in which the pressure of one component is much lower than that of the others, one can ensure that this component limits the generalized coefficients and , which describe the dependence of the growth rate on supersaturation (see Supplementary Table S1). This means that by measuring this dependence at a series of different equilibrium pressure sets, it is possible to determine the individual kinetic coefficients of each component and use them to describe the growth at any other combination of pressures.
The criterion 3) presented in Sec. 3.1 shows that to eliminate the undesirable nucleation of pure components, the crystal should be grown at a pressure set where the surface concentrations of adatoms of each component do not exceed its critical value. This phenomenon can also be used for the controlled formation of self-organized quantum dots or nanowire growth in certain systems (). In Sec. 3.2, BCF-growth when different types of adatoms are delivered to the surface by chemical reactions is considered. It differs from simple growth from vapor because chemical reactions alter the surface diffusion lengths of each component in comparison to thermal desorption, providing additional means to control the growth process. During a chemical reaction, the partial pressures of the reagents can be changed independently, allowing for individual control of supersaturation (chemical affinity), the mean free path of adatoms of each type, and their equilibrium concentrations. This enables control of the impact of each component on the generalized coefficients to ensure the maximum growth rate or to avoid undesirable phenomena, such as nucleation or instabilities. In Sec. 3.3, the impact of mechanical stress on the growth is considered. It is worth noting that mechanical stress are often present during the epitaxial growth of semiconductor thin films and heterostructures. The stress not only affect thermodynamics but also impact the kinetics of growth by determining the concentration of vacancies on the surface.
5 Future work
Although the mechanisms considered above are capable of describing numerous growth processes, significant gaps in knowledge and issues must be resolved to fully understand the multicomponent growth process. Further theoretical investigations could focus on several areas, including the impact of surface coverage of different components on mass transport (diffusion) both on the surface and in solution; the use of different adsorption isotherms for components to modify the relationship between surface supersaturation and pressure/concentration; the variation of Schwoebel barriers for each component, potentially affecting the distribution of adatoms; the detailed study of adatom incorporation into kinks, including the statistics of multicomponent kinks; the effect of surface reconstruction on diffusion mechanisms, which may depend on component concentration and growth rate; various morphological instabilities in multicomponent systems; the influence of impurities on incorporation or mass transfer of different components; anisotropic phenomena; the impact of convective flow of components during growth at higher supersaturations, when the step advancement rate cannot be ignored; and many other effects and their combinations.
6 Computer modeling
We note that in recent years, significant advancements have been made in high-performance modeling methods for multicomponent crystal structures and their growth. Some of these techniques include ab initio quantum mechanical (QM) modeling and stable structure prediction (), molecular dynamics (MD) simulations (; ), kinetic Monte Carlo (kMC) simulations (), phase-field modeling (PF) (; ), and the level-set approach (; ), among others. Each of these methods has its own temporal and spatial scales, strengths, and weaknesses (), but all can be used to deepen our understanding of the growth process in multicomponent systems. For example, the quantum mechanical approach can be used to find the stable crystalline structure with a given composition and its surface reconstruction under specific growth conditions. See, e.g., the USPEX code (; ). The latter may significantly affect the diffusion of the different components on the surface. The QM approach is also commonly used to find in silico the activation energy of surface diffusion and desorption of atoms of various components (), which is needed to use the formulas provided in the previous sections. QM can also be used to analyze the formation energies of kinks of different types (). Note that the concentration of kinks mainly determines the incorporation coefficients used in the theoretical models. The next scale is molecular dynamics, which may be used to analyze the statistics of kinks, steps, nucleation (; ), and the rates of incorporation of different components into the crystalline cell. We note that such phenomena can now be measured using experimental techniques (). MD is also often used for the prediction of diffusion coefficients of components in solution (), which may be useful in the case of growth from solution by the Gilmer-Ghez-Cabrera and Chernov mechanisms. Note that for an accurate description of the processes using the MD method, machine learning potentials are often used, which are trained on quantum modeling data, even for crystals containing several components. See, for example, the MLIP code (). The Monte-Carlo method is capable of describing larger time and spatial scales, up to micrometers and milliseconds, respectively, even on personal computers. Precise modeling requires knowledge of all the probabilities of atomistic events occurring on a surface, which can be obtained through quantum mechanics or molecular dynamics. Recently, the authors of the CrystalGrower package () implemented a unified kinetic three-dimensional partition model that can model the crystal habit and surface topology of any crystalline structure or MOF in different growth regimes, including nucleation and spiral growth. Note that various method combinations have also been developed, including a combination of Monte Carlo simulations with cellular automaton (), which, while less accurate, can significantly speed up the performance and detect various macro-phenomena on the surface, such as different types of step instabilities and growth regimes. The phase-field approach is also noteworthy, as it allows for the effective study of macrosurface kinetics and spiral growth (; ; ). To provide a complete in silico description of macroscopic crystal growth or dissolution, a multiscale approach is often used, involving all the aforementioned techniques (; ) from nano-to microscale, simultaneously coupled with finite-element models. In (), this approach was applied to the growth of multicomponent organic crystals. Furthermore, machine learning methods are being actively developed and hold promise for various applications, such as the search for stable crystalline structures (), determination of kinetic parameters of adatoms on the surface (), and optimization of growth chambers (). These new digital tools provide extensive information for the theoretical understanding of growth in multicomponent systems while also offering opportunities to discover new collective phenomena and effects and study them. They can also be successfully applied to solve the problems mentioned above.
7 Conclusion
In conclusion, we have summarized several growth mechanisms and effects involved in the growth of multicomponent crystals from vapors or solutions. Our findings demonstrate that growth can be described by the same equations as in the single-component case, with the use of generalized diffusion and incorporation coefficients based on the properties of each component. The results presented herein can be applied to a wide range of complex multicomponent systems, allowing the prediction of the growth rate based on growth conditions or, vice versa, the determination of individual kinetic properties of different components based on the dependence of the growth rate on conditions. We have also discussed the effects and peculiarities inherent to multicomponent systems, including the precipitation of nanoislands, effect of stresses, chemical reactions, and morphological instabilities. These phenomena can have a significant impact on the growth of multicomponent crystals or MOFs, making it crucial to understand their roles and the opportunities they provide for controlling the growth process in such complex systems. We presented approaches for accounting for new effects in theoretical equations, which enable the modification and extension of the theory to include additional factors. Lastly, we discussed some of the available numerical approaches and software packages, which make it possible to deepen knowledge and simplify the theoretical consideration of the crystal growth process.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Author contributions
The author confirms being the sole contributor of this work and has approved it for publication.
Funding
AR thanks the Russian Science Foundation (grant Nos. 21-72-00032) for the support of this work.
Acknowledgments
The author thanks prof. S.A. Kukushkin for fruitful discussions.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fchem.2023.1189729/full#supplementary-material
References
1
AlexandrovD. V.GalenkoP. K. (2021). A review on the theory of stable dendritic growth. Philosophical Trans. R. Soc. A379 (2205), 20200325. 10.1098/rsta.2020.0325
2
AltafA.HassanS.PejcicB.BaigN.HussainZ.SohailM. (2022). Recent progress in the design, synthesis and applications of chiral metal-organic frameworks. Front. Chem.10, 1014248. 10.3389/fchem.2022.1014248
3
AndersenM.PanosettiC.ReuterK. (2019). A practical guide to surface kinetic Monte Carlo simulations. Front. Chem.7, 202. 10.3389/fchem.2019.00202
4
AndersonM. W.Gebbie-RayetJ. T.HillA. R.FaridaN.AttfieldM. P.CubillasP.et al (2017). Predicting crystal growth via a unified kinetic three-dimensional partition model. Nature544, 456–459. 10.1038/nature21684
5
AnsariM. I. H.QurashiA.NazeeruddinM. K. (2018). Frontiers, opportunities, and challenges in perovskite solar cells: A critical review. Photochem. Rev.35, 1–24. 10.1016/j.jphotochemrev.2017.11.002
6
AnwarJ.ZahnD. (2011). Uncovering molecular processes in crystal nucleation and growth by using molecular simulation. Angew. Chem. Int. Ed.50 (9), 1996–2013. 10.1002/anie.201000463
7
BalesG. S.ZangwillA. (1990). Morphological instability of a terrace edge during step-flow growth. Phys. Rev. B41 (9), 5500–5508. 10.1103/physrevb.41.5500
8
BennemaP. (1984). Spiral growth and surface roughening: Developments since Burton, Cabrera and Frank. J. Cryst. Growth69, 182–197. 10.1016/0022-0248(84)90027-7
9
BurtonW. K.CabreraN.FrankF. C. (1951). The growth of crystals and the equilibrium structure of their surfaces. Philosophical Trans. R. Soc. A243, 299–358. 10.1098/rsta.1951.0006
10
CabreraN.LevineM. M. (1956). XLV. On the dislocation theory of evaporation of crystals. Philos. Mag.1 (4), 450–458. 10.1080/14786435608238124
11
CelebiA. T.JamaliS. H.BardowA.VlugtT. J.MoultosO. A. (2021). Finite-size effects of diffusion coefficients computed from molecular dynamics: A review of what we have learned so far. Mol. Simul.47 (10-11), 831–845. 10.1080/08927022.2020.1810685
12
ChernovA. A. (2004). Notes on interface growth kinetics 50 years after Burton, Cabrera and Frank. J. Cryst. Growth264, 499–518. 10.1016/j.jcrysgro.2003.12.076
13
ChernovA. A. (1961). The spiral growth of crystals. Sov. Phys. Uspekhi4 (1), 116–148. 10.1070/pu1961v004n01abeh003328
14
DongZ.ZhangL.WangS.LuoL. (2020). Direct visualization of dynamic atomistic processes of Cu2O crystal growth through gas-solid reaction. N Energy70, 104527. 10.1016/j.nanoen.2020.104527
15
DubrovskiiV. G. (2014). Nucleation theory and growth of nanostructures. Berlin: Springer, 601.
16
EltsE.GreinerM.BriesenH. (2017). In silico prediction of growth and dissolution rates for organic molecular crystals: A multiscale approach. Crystals7, 288. 10.3390/cryst7100288
17
FreundR.ZarembaO.ArnautsG.AmelootR.SkorupskiG.DincăM.et al (2021). The current status of MOF and COF applications. Angew. Chem. Int. Ed.60 (45), 23975–24001. 10.1002/anie.202106259
18
GibouF.FedkiwR.OsherS. (2018). A review of level-set methods and some recent applications. J. Comput. Phys.353, 82–109. 10.1016/j.jcp.2017.10.006
19
GilmerG. H.GhezR.CabreraN. (1971). An analysis of combined surface and volume diffusion processes in crystal growth. J. Cryst. Growth15, 79–93. 10.1016/0022-0248(71)90027-3
20
GlassC. W.OganovA. R.HansenN. (2006). USPEX—evolutionary crystal structure prediction. Comput. Phys. Commun.175 (11-12), 713–720. 10.1016/j.cpc.2006.07.020
21
GomezH.BuresM.MoureA. (2019). A review on computational modelling of phase-transition problems. Philosophical Trans. R. Soc. A377 (2143), 20180203. 10.1098/rsta.2018.0203
22
HirthJ. P. (1965). “The kinetic and thermodynamic properties of surfaces: Surface energy, surface diffusion, and growth and dissolution or evaporation processes,” in Energetics in metallurgical phenomena (Philadelphia, PA: Gordon and Breach), 2.
23
JagadammaL. K.WangS. (2021). Wide-bandgap halide perovskites for indoor photovoltaics. Front. Chem.9, 632021. 10.3389/fchem.2021.632021
24
JeongH.-C.WilliamsE. D. (1999). Steps on surfaces: Experiment and theory. Surf. Sci. Rep.34, 171–294. 10.1016/S0167-5729(98)00010-7
25
JoswiakM. N.PetersB.DohertyM. F. (2018). In silico crystal growth rate prediction for NaCl from aqueous solution. Cryst. Growth and Des.18 (10), 6302–6306. 10.1021/acs.cgd.8b01184
26
KangS.JeongW.HongC.HwangS.YoonY.HanS. (2022). Accelerated identification of equilibrium structures of multicomponent inorganic crystals using machine learning potentials. npj Comput. Mater.8 (1), 108. 10.1038/s41524-022-00792-w
27
KarthikaS.RadhakrishnanT. K.KalaichelviP. (2016). A review of classical and nonclassical nucleation theories. Cryst. Growth and Des.16 (11), 6663–6681. 10.1021/acs.cgd.6b00794
28
KosolobovS. (2019). Subsurface diffusion in crystals and effect of surface permeability on the atomic step motion. Sci. Rep.9 (1), 13428. 10.1038/s41598-019-49681-1
29
KourR.AryaS.VermaS.SinghA.MahajanP.KhoslaA. (2019). Recent advances and challenges in indium gallium nitride (InxGa1-xN) materials for solid state lighting. ECS J. Solid State Sci. Technol.9 (1), 015011. 10.1149/2.0292001jss
30
KukushkinS. A.OsipovA. V.RedkovA. V. (2014). Morphological stability criterion for a spherical crystallization front in a multicomponent system with chemical reactions. Phys. Solid State56, 2530–2536. 10.1134/s106378341412018x
31
KukushkinS. A.OsipovA. V. (1998). Thin-film condensation processes. Physics-Uspekhi41 (10), 983–1014. 10.1070/pu1998v041n10abeh000461
32
KuvadiaZ. B.DohertyM. F. (2011). Spiral growth model for faceted crystals of non-centrosymmetric organic molecules grown from solution. Cryst. growth and Des.11 (7), 2780–2802. 10.1021/cg101560u
33
KvashninA. G.KvashninD. G.OganovA. R. (2019). Novel unexpected reconstructions of (100) and (111) surfaces of NaCl: Theoretical prediction. Sci. Rep.9 (1), 14267. 10.1038/s41598-019-50548-8
34
MartynecT.KarapanagiotisC.KlappS. H.KowarikS. (2021). Machine learning predictions of surface migration barriers in nucleation and non-equilibrium growth. Commun. Mater.2 (1), 90. 10.1038/s43246-021-00188-1
35
MillerW. (2015). “Simulation of epitaxial growth by means of density functional theory, kinetic Monte Carlo, and phase field methods,” in Handbook of crystal growth (Amsterdam: Elsevier), 521–559.
36
MiuraH.KobayashiR. (2015). Phase-field modeling of step dynamics on growing crystal surface: Direct integration of growth units to step front. Cryst. Growth and Des.15, 2165–2175. 10.1021/cg501806d
37
MiuraH. (2015). Phase-field modeling of step dynamics on growing crystal surface: Step pinning induced by impurities. Cryst. Growth and Des.15, 4142–4148. 10.1021/acs.cgd.5b00762
38
MullinsW. W.SekerkaR. F. (1963). Morphological stability of a particle growing by diffusion or heat flow. J. Appl. Phys.34 (2), 323–329. 10.1063/1.1702607
39
Myers-BeaghtonA. K.VvedenskyD. D. (1991). Generalized Burton-Cabrera-Frank theory for growth and equilibration on stepped surfaces. Phys. Rev. A44 (4), 2457–2468. 10.1103/physreva.44.2457
40
NieY.BartonA. T.AddouR.ZhengY.WalshL. A.EichfeldS. M.et al (2018). Dislocation driven spiral and non-spiral growth in layered chalcogenides. Nanoscale10, 15023–15034. 10.1039/c8nr02280a
41
NovikovI. S.GubaevK.PodryabinkinE. V.ShapeevA. V. (2020). The MLIP package: Moment tensor potentials with MPI and active learning. Mach. Learn. Sci. Technol.2 (2), 025002. 10.1088/2632-2153/abc9fe
42
Pierre-LouisO. (2003). Phase field models for step flow. Phys. Rev. E68, 021604. 10.1103/physreve.68.021604
43
PimpinelliA.VillainJ. (1994). What does an evaporating surface look like?Phys. A Stat. Mech. its Appl.204 (1-4), 521–542. 10.1016/0378-4371(94)90446-4
44
PolitiP.GrenetG.MartyA.PonchetA.VillainJ. (2000). Instabilities in crystal growth by atomic or molecular beams. Phys. Rep.324 (5-6), 271–404. 10.1016/s0370-1573(99)00046-0
45
Ramirez ReinaT.IvanovaS.CentenoM. A.OdriozolaJ. A. (2013). Low-temperature CO oxidation on multicomponent gold-based catalysts. Front. Chem.1, 12. 10.3389/fchem.2013.00012
46
RatschC.GyureM. F.CaflischR. E.GibouF.PetersenM.KangM.et al (2002). Level-set method for island dynamics in epitaxial growth. Phys. Rev. B65 (19), 195403. 10.1103/physrevb.65.195403
47
RedkovA. V.KukushkinS. A. (2020). Development of burton-cabrera-frank theory for the growth of a non-kossel crystal via chemical reaction. Cryst. Growth&Design20 (4), 2590–2601. 10.1021/acs.cgd.9b01721
48
RedkovA. V.KukushkinS. A. (2021). Dynamic interaction of steps and nanoislands during growth of a multicomponent crystal. Cryst. Growth&Design21 (9), 4914–4926. 10.1021/acs.cgd.1c00349
49
RedkovA. V.KukushkinS. A.OsipovA. V. (2020b). Growth of faceted pores in a multi-component crystal by applying mechanical stress. CrystEngComm22 (32), 5280–5288. 10.1039/d0ce00888e
50
RedkovA. V.KukushkinS. A.OsipovA. V. (2020a). Spiral growth of a multicomponent crystal from vapor of its components. J. Cryst. Growth548, 125845. 10.1016/j.jcrysgro.2020.125845
51
RedkovA. V.KukushkinS. A. (2022). Theoretical aspects of the growth of a non-kossel crystal from vapours: The role of advacancies. Faraday Discuss.235, 362–382. 10.1039/d1fd00083g
52
RedkovA. V.OsipovA. V.KukushkinS. A. (2015). Stability of the surface of an elastically strained multicomponent film in a system with chemical reactions. Phys. Solid State57, 2524–2531. 10.1134/s106378341512029x
53
RostM. J.JacobseL.KoperM. T. (2019). The dualism between adatom- and vacancy-based single crystal growth models. Nat. Commun.10, 5233–5237. 10.1038/s41467-019-13188-0
54
SalvalaglioM.PeregoC.GibertiF.MazzottiM.ParrinelloM. (2015). Molecular-dynamics simulations of urea nucleation from aqueous solution. Proc. Natl. Acad. Sci.112, E6–E14. 10.1073/pnas.1421192111
55
SangwalK. (1996). Effects of impurities on crystal growth processes. Prog. Cryst. Growth Charact. Mater.32, 3–43. 10.1016/0960-8974(96)00008-3
56
SchimmelS.SunW.DropkaN. (2022). Artificial intelligence for crystal growth and characterization. Crystals12 (9), 1232. 10.3390/cryst12091232
57
SitnikovS. V.LatyshevA. V.KosolobovS. S. (2017). Advacancy-mediated atomic steps kinetics and two-dimensional negative island nucleation on ultra-flat Si (111) surface. J. Cryst. Growth457, 196–201. 10.1016/j.jcrysgro.2016.05.048
58
SossoG. C.ChenJ.CoxS. J.FitznerM.PedevillaP.ZenA.et al (2016). Crystal nucleation in liquids: Open questions and future challenges in molecular dynamics simulations. Chem. Rev.116 (12), 7078–7116. 10.1021/acs.chemrev.5b00744
59
TanA. K.HamzahN. A.AhmadM. A.NgS. S.HassanZ. (2022). Recent advances and challenges in the MOCVD growth of indium gallium nitride: A brief review. Mater. Sci. Semicond. Process.143, 106545. 10.1016/j.mssp.2022.106545
60
TohidiF.HolaghS. G.ChitsazA. (2022). Thermoelectric generators: A comprehensive review of characteristics and applications. Appl. Therm. Eng.201, 117793. 10.1016/j.applthermaleng.2021.117793
61
TroyanI. A.SemenokD. V.IvanovaA. G.KvashninA. G.ZhouD.SadakovA. V.et al (2022). High-temperature superconductivity in hydrides. Phys. - Uspekhi65 (7), 748–761. 10.3367/ufne.2021.05.039187
62
UwahaM. (2016). Introduction to the BCF theory. Prog. Cryst. Growth Charact. Mater.62, 58–68. 10.1016/j.pcrysgrow.2016.04.002
63
Van Der EerdenJ. P. (1982). The advance velocity of steps under the influence of volume and surface diffusion, by direct and indirect incorporation of growth units. J. Cryst. Growth56, 174–188. 10.1016/0022-0248(82)90027-6
64
Van der HoekB.Van Der EerdenJ. P.BennemaP.SunagawaI. (1982). The influence of stress on spiral growth. J. Cryst. Growth58, 365–380. 10.1016/0022-0248(82)90284-6
65
VekilovP. G. (2007). What determines the rate of growth of crystals from solution?Cryst. Growth&Design7 (12), 2796–2810. 10.1021/cg070427i
66
VvedenskyD. D. (2004). Multiscale modelling of nanostructures. J. Phys. Condens. Matter16, R1537–R1576. 10.1088/0953-8984/16/50/r01
67
WangH. H.HouL.LiY. Z.JiangC. Y.WangY. Y.ZhuZ. (2017). Porous MOF with highly efficient selectivity and chemical conversion for CO2. ACS Appl. Mater. interfaces9, 17969–17976. 10.1021/acsami.7b03835
68
WangP.AnJ.YeZ.CaiW.ZhengX. (2022). Cu-based multicomponent metallic compound materials as electrocatalyst for water splitting. Front. Chem.10, 913874. 10.3389/fchem.2022.913874
69
WangQ.AstrucD. (2019). State of the art and prospects in metal-organic framework (MOF)-based and MOF-derived nanocatalysis. Chem. Rev.120 (2), 1438–1511. 10.1021/acs.chemrev.9b00223
70
WonY. S.LeeJ.KimC. S.ParkS. S. (2009). Computational study of adsorption, diffusion, and dissociation of precursor species on the GaN (0001) surface during GaN MOCVD. Surf. Sci.603 (4), L31–L34. 10.1016/j.susc.2009.01.008
71
Załuska-KoturM.PopovaH.TonchevV. (2021). Step bunches, nanowires and other vicinal “creatures”—ehrlich–schwoebel effect by cellular automata. Crystals11, 1135. 10.3390/cryst11091135
Summary
Keywords
multicomponent, crystal, growth, theory, BCF, Chernov, Gilmer-Ghez-Cabrera, mechanism
Citation
Redkov A (2023) Spiral growth of multicomponent crystals: theoretical aspects. Front. Chem. 11:1189729. doi: 10.3389/fchem.2023.1189729
Received
19 March 2023
Accepted
28 April 2023
Published
12 May 2023
Volume
11 - 2023
Edited by
Xin Huang, Tianjin University, China
Reviewed by
Xiaoxi Yu, China University of Petroleum (Huadong), China
Updates
Copyright
© 2023 Redkov.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Alexey Redkov, avredkov@gmail.com
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.