1
|
Luis EEM, Ferreira SC, de Assis TA. Bifractality in the one-dimensional Wolf-Villain model. Phys Rev E 2024; 110:L012801. [PMID: 39161014 DOI: 10.1103/physreve.110.l012801] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/13/2024] [Accepted: 06/10/2024] [Indexed: 08/21/2024]
Abstract
We introduce a multifractal optimal detrended fluctuation analysis to study the scaling properties of the one-dimensional Wolf-Villain (WV) model for surface growth. This model produces coarsened surface morphologies for long timescales (up to 10^{9} monolayers) and its universality class remains an open problem. Our results for the multifractal exponent τ(q) reveal an effective local roughness exponent consistent with a transient given by the molecular beam epitaxy (MBE) growth regime and Edwards-Wilkinson (EW) universality class for negative and positive q values, respectively. Therefore, although the results corroborate that long-wavelength fluctuations belong to the EW class in the hydrodynamic limit, as conjectured in the recent literature, a bifractal signature of the WV model with an MBE regime at short wavelengths was observed.
Collapse
Affiliation(s)
| | | | - Thiago A de Assis
- Instituto de Física, Universidade Federal Fluminense, Avenida Litorânea s/n, 24210-340, Niterói, RJ, Brazil
- Instituto de Física, Universidade Federal da Bahia, Campus Universitário da Federação, Rua Barão de Jeremoabo s/n, 40170-115, Salvador, BA, Brazil
| |
Collapse
|
2
|
Shin J, Lee IJ. Surface Evolution of Polymer Films Grown by Vapor Deposition: Growth of Local and Global Slopes of Interfaces. Polymers (Basel) 2024; 16:1535. [PMID: 38891479 PMCID: PMC11175125 DOI: 10.3390/polym16111535] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/29/2024] [Revised: 05/20/2024] [Accepted: 05/27/2024] [Indexed: 06/21/2024] Open
Abstract
The kinetic roughening of polymer films grown by vapor deposition polymerization was analyzed using the widely accepted classification framework of "generic scaling ansatz" given for the structure factor. Over the past two decades, this method has played a pivotal role in classifying diverse forms of dynamic scaling and understanding the mechanisms driving interface roughening. The roughness exponents of the polymer films were consistently determined as α=1.25±0.09, αloc=0.73±0.02, and αs=0.99±0.06. However, the inability to unambiguously assign these roughness exponent values to a specific scaling subclass prompts the proposal of a practical alternative. This report illustrates how all potential dynamic scaling can be consistently identified and classified based on the relationship between two temporal scaling exponents measured in real space: the average local slope and the global slope of the interface. The intrinsic anomalous roughening class is conclusively assigned to polymer film growth characterized by anomalous "native (background slope-removed) local height fluctuations". Moreover, the new analysis reveals that interfaces exhibiting anomalous scaling, previously classified as intrinsic anomalous roughening, could potentially belong to the super-rough class, particularly when the spectral roughness exponent αs is equal to 1.
Collapse
Affiliation(s)
| | - I. J. Lee
- Department of Physics, Research Institute of Physics and Chemistry, Jeonbuk National University, Jeonju 54896, Republic of Korea
| |
Collapse
|
3
|
Martynec T, Klapp SHL. Modeling of nonequilibrium surface growth by a limited-mobility model with distributed diffusion length. Phys Rev E 2019; 100:033307. [PMID: 31639962 DOI: 10.1103/physreve.100.033307] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/20/2019] [Indexed: 11/07/2022]
Abstract
Kinetic Monte Carlo (KMC) simulations are a well-established numerical tool to investigate the time-dependent surface morphology in molecular beam epitaxy experiments. In parallel, simplified approaches such as limited mobility (LM) models characterized by a fixed diffusion length have been studied. Here we investigate an extended LM model to gain deeper insight into the role of diffusional processes concerning the growth morphology. Our model is based on the stochastic transition rules of the Das Sarma-Tamborena model but differs from the latter via a variable diffusion length. A first guess for this length can be extracted from the saturation value of the mean-squared displacement calculated from short KMC simulations. Comparing the resulting surface morphologies in the sub- and multilayer growth regime to those obtained from KMC simulations, we find deviations which can be cured by adding fluctuations to the diffusion length. This mimics the stochastic nature of particle diffusion on a substrate, an aspect which is usually neglected in LM models. We propose to add fluctuations to the diffusion length by choosing this quantity for each adsorbed particle from a Gaussian distribution, where the variance of the distribution serves as a fitting parameter. We show that the diffusional fluctuations have a huge impact on cluster properties during submonolayer growth as well as on the surface profile in the high coverage regime. The analysis of the surface morphologies on one- and two-dimensional substrates during sub- and multilayer growth shows that the LM model can produce structures that are indistinguishable to the ones from KMC simulations at arbitrary growth conditions.
Collapse
Affiliation(s)
- Thomas Martynec
- Institut für Theoretische Physik, Technische Universität Berlin, Hardenbergstr. 36, 10623 Berlin, Germany
| | - Sabine H L Klapp
- Institut für Theoretische Physik, Technische Universität Berlin, Hardenbergstr. 36, 10623 Berlin, Germany
| |
Collapse
|
4
|
Pereira AJ, Alves SG, Ferreira SC. Effects of a kinetic barrier on limited-mobility interface growth models. Phys Rev E 2019; 99:042802. [PMID: 31108608 DOI: 10.1103/physreve.99.042802] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/13/2019] [Indexed: 11/07/2022]
Abstract
The role played by a kinetic barrier originated by out-of-plane step edge diffusion, introduced by Leal et al. [J. Phys.: Condens. Matter 23, 292201 (2011)JCOMEL0953-898410.1088/0953-8984/23/29/292201], is investigated in the Wolf-Villain and Das Sarma-Tamborenea models with short-range diffusion. Using large-scale simulations, we observe that this barrier is sufficient to produce growth instability, forming quasiregular mounds in one and two dimensions. The characteristic surface length saturates quickly indicating a uncorrelated growth of the three-dimensional structures, which is also confirmed by a growth exponent β=1/2. The out-of-plane particle current shows a large reduction of the downward flux in the presence of the kinetic barrier enhancing, consequently, the net upward diffusion and the formation of three-dimensional self-assembled structures.
Collapse
Affiliation(s)
- Anderson J Pereira
- Departamento de Física, Universidade Federal de Viçosa, Minas Gerais, 36570-900, Viçosa, Brazil
| | - Sidiney G Alves
- Departamento de Estatística, Física e Matemática, Campus Alto Paraopeba, Universidade Federal de São João Del-Rei, 36420-000, Ouro Branco, MG, Brazil
| | - Silvio C Ferreira
- Departamento de Física, Universidade Federal de Viçosa, Minas Gerais, 36570-900, Viçosa, Brazil.,National Institute of Science and Technology for Complex Systems, 22290-180, Rio de Janeiro, Brazil
| |
Collapse
|
5
|
de Assis TA, Reis FDAA. Smoothening in thin-film deposition on rough substrates. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2015; 92:052405. [PMID: 26651710 DOI: 10.1103/physreve.92.052405] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/25/2015] [Indexed: 06/05/2023]
Abstract
The evolution of the surface roughness W of a thin film deposited on a rough substrate is studied with a model of temperature-activated adatom diffusion, irreversible lateral aggregation, and no step energy barrier, in which the main parameter is the ratio R of diffusion and deposition rates. At sufficiently low temperatures (R≲10), the average number of adatom steps after adsorption is very small, thus W monotonically increases with time t due to an approximately uncorrelated deposition at short times. If the temperature is not very low (R∼10(3) or larger), smoothening occurs at short times and the Villain-Lai-Das Sarma (VLDS) growth equation governs the long time roughening, which is attained after a crossover time t(c) that increases with the correlation length ξ(i) of the substrate. Scaling arguments predict the dependence of t(c) on temperature and on the substrate production time and the scaling relation for the difference between the roughness of films deposited on rough and flat substrates, in good agreement with numerical results. The effect of temperature is not a direct extension of previous results on flat substrates because the short wavelength fluctuations delay the formation of terraces. For this reason, the effective energy obtained from the dependence of t(c) on R is 40% of the energy of activated adatom diffusion. A scaling law for the initial smoothening is proposed as W/W(i)=Ψ(t/t(c1)), with a crossover time t(c1)≡R(-θ)ξ(i)(z), where W(i) is the substrate roughness, θ≈0.4, and z is the VLDS dynamical exponent. It provides good data collapse if W is not very small and is suggested to be tested experimentally.
Collapse
Affiliation(s)
- T A de Assis
- Instituto de Física, Universidade Federal da Bahia, Campus Universitário da Federação, Rua Barão de Jeremoabo s/n, 40170-115, Salvador, BA, Brazil
| | - F D A Aarão Reis
- Instituto de Física, Universidade Federal Fluminense, Avenida Litorânea s/n, 24210-340 Niterói RJ, Brazil
| |
Collapse
|
6
|
Aarão Reis FDA. Normal dynamic scaling in the class of the nonlinear molecular-beam-epitaxy equation. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2013; 88:022128. [PMID: 24032796 DOI: 10.1103/physreve.88.022128] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/30/2013] [Revised: 07/31/2013] [Indexed: 06/02/2023]
Abstract
The scaling of local height fluctuations is studied numerically in lattice growth models of the class of the nonlinear stochastic equation of Villain-Lai-Das Sarma (VLDS) in substrate dimensions d=1 and 2. In d=1, the average local slopes of the conserved restricted solid-on-solid (CRSOS) models converge to a finite value in the long-time limit, with power-law corrections in time whose exponents are close to 0.1. Other VLDS models in d=1, such as that of Das Sarma and Tamborenea, show a divergence of local slopes up to 10(6) monolayers, typical of anomalous roughening, but a comparison of roughness distributions shows that they scale as the linear fourth-order growth equation in those time scales. Normal scaling is also obtained in a modified VLDS equation with instability suppression, in contrast to recent numerical works. In d=2, a CRSOS model and a model with lateral aggregation of diffusing particles show normal scaling of the local slopes, also with small correction exponents. These results consistently show that the VLDS class has normal dynamic scaling in d=1 and 2, in agreement with the theoretical predictions of Phys. Rev. Lett. 94, 166103 (2005), and they show that the apparently anomalous features observed in previous works are effects of large scaling correction terms or crossover effects.
Collapse
Affiliation(s)
- F D A Aarão Reis
- Instituto de Física, Universidade Federal Fluminense, Avenida Litorânea s/n, 24210-340 Niterói RJ, Brazil
| |
Collapse
|
7
|
Xun Z, Tang G, Han K, Xia H, Hao D, Li Y. Asymptotic dynamic scaling behavior of the (1+1)-dimensional Wolf-Villain model. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2012; 85:041126. [PMID: 22680438 DOI: 10.1103/physreve.85.041126] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/15/2012] [Revised: 03/28/2012] [Indexed: 06/01/2023]
Abstract
Extensive kinetic Monte Carlo simulations are presented for the Wolf-Villain model in (1+1) dimensions. Asymptotic dynamic scaling is found for lattice sizes L≥2048. The exponents obtained from our simulations, α=0.50±0.02 and β=0.25±0.02, are in excellent agreement with the exact values α=1/2 and β=1/4 for the one-dimensional Edwards-Wilkinson equation. Our findings explain the widespread discrepancies of previous reports for exponents of the Wolf-Villain model in (1+1) dimensions, and the results are also consistent with the theoretical predictions of López et al. [J. M. López, M. Castro, and R. Gallego, Phys. Rev. Lett. 94, 166103 (2005)].
Collapse
Affiliation(s)
- Zhipeng Xun
- Department of Physics, China University of Mining and Technology, Xuzhou 221116, People's Republic of China.
| | | | | | | | | | | |
Collapse
|
8
|
Kanjanaput W, Limkumnerd S, Chatraphorn P. Growth instability due to lattice-induced topological currents in limited-mobility epitaxial growth models. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2010; 82:041607. [PMID: 21230287 DOI: 10.1103/physreve.82.041607] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/09/2010] [Revised: 06/25/2010] [Indexed: 05/30/2023]
Abstract
The energetically driven Ehrlich-Schwoebel barrier had been generally accepted as the primary cause of the growth instability in the form of quasiregular moundlike structures observed on the surface of thin film grown via molecular-beam epitaxy (MBE) technique. Recently the second mechanism of mound formation was proposed in terms of a topologically induced flux of particles originating from the line tension of the step edges which form the contour lines around a mound. Through large-scale simulations of MBE growth on a variety of crystalline lattice planes using limited-mobility, solid-on-solid models introduced by Wolf-Villain and Das Sarma-Tamborenea in 2+1 dimensions, we show that there exists a topological uphill particle current with strong dependence on specific lattice crystalline structure. Without any energetically induced barriers, our simulations produce spectacular mounds very similar, in some cases, to what have been observed in many recent MBE experiments. On a lattice where these currents cease to exist, the surface appears to be scale invariant, statistically rough as predicted by the conventional continuum growth equation.
Collapse
Affiliation(s)
- Wittawat Kanjanaput
- Department of Physics, Faculty of Science, Chulalongkorn University, Phayathai Road, Patumwan, Bangkok 10330, Thailand.
| | | | | |
Collapse
|
9
|
Pang NN, Tzeng WJ. Extensive studies on linear growth processes with spatiotemporally correlated noise in arbitrary substrate dimensions. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2010; 82:031605. [PMID: 21230084 DOI: 10.1103/physreve.82.031605] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/06/2008] [Revised: 06/25/2010] [Indexed: 05/30/2023]
Abstract
An extensive analytical and numerical study on a class of growth processes with spatiotemporally correlated noise in arbitrary dimension is undertaken. In addition to the conventional investigation on the interface morphology and interfacial widths, we pay special attention to exploring the characteristics of the slope-slope correlation function S(r,t) and the [Q]-th degree residual local interfacial width w[Q](l,t), whose importance has been somewhat overlooked in the literature. Based on the above analysis, we give a plausible theoretical explanation about the various experimental observations of kinetically and thermodynamically unstable surface growth. Furthermore, through explicit examples, we show that the statistical methods of calculating the exponents (including the dynamic exponent z, the global roughness exponent α, and the local roughness exponent α(loc)), based on the scaling of S(r,t) and w[Q](l,t), are very reliable and rarely influenced by the finite time and/or finite-size effects. Another important issue we focus on in this paper is related to numerical calculation. For the specific class of growth processes discussed in this paper, we develop a very efficient and accurate algorithm for numerical calculation of the dynamics of interface configuration, the structure factor, the various correlation functions, the interfacial width and its variants in arbitrary dimensions, even with very large system size and very late time. The proposed systematical algorithm can be easily generalized to other linear processes and some special nonlinear processes.
Collapse
Affiliation(s)
- Ning-Ning Pang
- Department of Physics and Center for Theoretical Sciences, National Taiwan University, Taipei, Taiwan
| | | |
Collapse
|
10
|
Odor G, Liedke B, Heinig KH. Surface pattern formation and scaling described by conserved lattice gases. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2010; 81:051114. [PMID: 20866192 DOI: 10.1103/physreve.81.051114] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/15/2009] [Revised: 02/26/2010] [Indexed: 05/29/2023]
Abstract
We extend our 2+1 -dimensional discrete growth model [Odor, Phys. Rev. E 79, 021125 (2009)] with conserved, local exchange dynamics of octahedra, describing surface diffusion. A roughening process was realized by uphill diffusion and curvature dependence. By mapping the slopes onto particles, two-dimensional nonequilibrium binary lattice model emerges, in which the (smoothing or roughening) surface diffusion can be described by attracting or repelling motion of oriented dimers. The binary representation allows simulations on very large size and time scales. We provide numerical evidence for Mullins-Herring or molecular-beam epitaxy class scaling of the surface width. The competition of inverse Mullins-Herring diffusion with a smoothing deposition, which corresponds to a Kardar-Parisi-Zhang (KPZ) process, generates different patterns: dots or ripples. We analyze numerically the scaling and wavelength growth behavior in these models. In particular, we confirm by large size simulations that the KPZ type of scaling is stable against the addition of this surface diffusion, hence this is the asymptotic behavior of the Kuramoto-Sivashinsky equation as conjectured by field theory in two dimensions, but has been debated numerically. If very strong, normal surface diffusion is added to a KPZ process, we observe smooth surfaces with logarithmic growth, which can describe the mean-field behavior of the strong-coupling KPZ class. We show that ripple coarsening occurs if parallel surface currents are present, otherwise logarithmic behavior emerges.
Collapse
Affiliation(s)
- Géza Odor
- Research Institute for Technical Physics and Materials Science, P.O. Box 49, H-1525 Budapest, Hungary
| | | | | |
Collapse
|
11
|
Aarão Reis FDA. Dynamic scaling in thin-film growth with irreversible step-edge attachment. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2010; 81:041605. [PMID: 20481733 DOI: 10.1103/physreve.81.041605] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/15/2009] [Revised: 03/23/2010] [Indexed: 05/29/2023]
Abstract
We study dynamic scaling in a model with collective diffusion (CD) of isolated atoms in terraces and irreversible aggregation at step edges. Simulations are performed in two-dimensional substrates with several diffusion to deposition ratios R identical with D/F. Data collapse of scaled roughness distributions confirms that this model is in the class of the fourth-order nonlinear growth equation by Villain, Lai, and Das Sarma (VLDS) with negligible finite-size effects, while estimates of scaling exponents show some discrepancies. This result is consistent with the prediction of a recent renormalization group approach and improves previous numerical works on related models. The roughness follows dynamic scaling as W=Lalpha/R1/2f(xi/L), with correlation length xi=(Rt)1/z, where z is the dynamic exponent. We also propose a limited mobility (LM) model where the incident atom executes up to G steps before a new atom is adsorbed, and irreversibly aggregates at step edges. This model is also shown to belong to the VLDS class. The size of the plateaus in the film surface increases as G1/2 and the lateral correlation scales as G1/2t1/z. The time evolution of the roughness reproduces that of the CD model if an equivalent parameter G approximately R2/z is chosen. This suggests the possibility of using LM models with tunable diffusion length to simulate processes with simultaneous diffusion of many atoms. A scaling approach is used to justify exponent values and dynamic relations for both models, including the significant decrease of surface roughness as R or G increases.
Collapse
Affiliation(s)
- F D A Aarão Reis
- Instituto de Física, Universidade Federal Fluminense, Avenida Litorânea s/n, 24210-340 Niterói, RJ, Brazil.
| |
Collapse
|
12
|
Haselwandter CA, Vvedensky DD. Renormalization of stochastic lattice models: epitaxial surfaces. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2008; 77:061129. [PMID: 18643239 DOI: 10.1103/physreve.77.061129] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/28/2007] [Indexed: 05/26/2023]
Abstract
We present the application of a method [C. A. Haselwandter and D. D. Vvedensky, Phys. Rev. E 76, 041115 (2007)] for deriving stochastic partial differential equations from atomistic processes to the morphological evolution of epitaxial surfaces driven by the deposition of new material. Although formally identical to the one-dimensional (1D) systems considered previously, our methodology presents substantial additional technical issues when applied to two-dimensional (2D) surfaces. Once these are addressed, subsequent coarse-graining is accomplished as before by calculating renormalization-group (RG) trajectories from initial conditions determined by the regularized atomistic models. Our applications are to the Edwards-Wilkinson (EW) model [S. F. Edwards and D. R. Wilkinson, Proc. R. Soc. London, Ser. A 381, 17 (1982)], the Wolf-Villain (WV) model [D. E. Wolf and J. Villain, Europhys. Lett. 13, 389 (1990)], and a model with concurrent random deposition and surface diffusion. With our rules for the EW model no appreciable crossover is obtained for either 1D or 2D substrates. For the 1D WV model, discussed previously, our analysis reproduces the crossover sequence known from kinetic Monte Carlo (KMC) simulations, but for the 2D WV model, we find a transition from smooth to unstable growth under repeated coarse-graining. Concurrent surface diffusion does not change this behavior, but can lead to extended transient regimes with kinetic roughening. This provides an explanation of recent experiments on Ge(001) with the intriguing conclusion that the same relaxation mechanism responsible for ordered structures during the early stages of growth also produces an instability at longer times that leads to epitaxial breakdown. The RG trajectories calculated for concurrent random deposition and surface diffusion reproduce the crossover sequences observed with KMC simulations for all values of the model parameters, and asymptotically always approach the fixed point corresponding to the equation proposed by Villain [J. Phys. I 1, 19 (1991)] and by Lai and Das Sarma [Phys. Rev. Lett. 66, 2899 (1991)]. We conclude with a discussion of the application of our methodology to other growth settings.
Collapse
Affiliation(s)
- Christoph A Haselwandter
- Department of Physics, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | | |
Collapse
|
13
|
Haselwandter CA, Vvedensky DD. Renormalization of stochastic lattice models: basic formulation. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2007; 76:041115. [PMID: 17994944 DOI: 10.1103/physreve.76.041115] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/11/2007] [Indexed: 05/25/2023]
Abstract
We describe a general method for the multiscale analysis of stochastic lattice models. Beginning with a lattice Langevin formulation of site fluctuations, we derive stochastic partial differential equations by regularizing the transition rules of the model. Subsequent coarse graining is accomplished by calculating renormalization-group (RG) trajectories from initial conditions determined by the regularized atomistic models. The RG trajectories correspond to hierarchies of continuum equations describing lattice models over expanding length and time scales. These continuum equations retain a quantitative connection over different scales, as well as to the underlying atomistic dynamics. This provides a systematic method for the derivation of continuum equations from the transition rules of lattice models for any length and time scales. As an illustration we consider the one-dimensional (1D) Wolf-Villain (WV) model [Europhys. Lett. 13, 389 (1990)]. The RG analysis of this model, which we develop in detail, is generic and can be applied to a wide range of conservative lattice models. The RG trajectory of the 1D WV model shows a complex crossover sequence of linear and nonlinear stochastic differential equations, which is in excellent agreement with kinetic Monte Carlo simulations of this model. We conclude by discussing possible applications of the multiscale method described here to other nonequilibrium systems.
Collapse
|
14
|
Haselwandter CA, Vvedensky DD. Multiscale theory of fluctuating interfaces: renormalization of atomistic models. PHYSICAL REVIEW LETTERS 2007; 98:046102. [PMID: 17358788 DOI: 10.1103/physrevlett.98.046102] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/29/2005] [Indexed: 05/14/2023]
Abstract
We describe a framework for the multiscale analysis of atomistic surface processes which we apply to a model of homoepitaxial growth with deposition according to the Wolf-Villain model and concurrent surface diffusion. Coarse graining is accomplished by calculating renormalization-group (RG) trajectories from initial conditions determined by the regularized atomistic theory. All of the crossover and asymptotic scaling regimes known from computer simulations are obtained, but we also find that two-dimensional substrates show an intriguing transition from smooth to mounded morphologies along the RG trajectory.
Collapse
|
15
|
Chua ALS, Haselwandter CA, Baggio C, Vvedensky DD. Langevin equations for fluctuating surfaces. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2005; 72:051103. [PMID: 16383589 DOI: 10.1103/physreve.72.051103] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/27/2005] [Indexed: 05/05/2023]
Abstract
Exact Langevin equations are derived for the height fluctuations of surfaces driven by the deposition of material from a molecular beam. We consider two types of model: deposition models, where growth proceeds by the deposition and instantaneous local relaxation of particles, with no subsequent movement, and models with concurrent random deposition and surface diffusion. Starting from a Chapman-Kolmogorov equation the deposition, relaxation, and hopping rules of these models are first expressed as transition rates within a master equation for the joint height probability density function. The Kramers-Moyal-van Kampen expansion of the master equation in terms of an appropriate "largeness" parameter yields, according to a limit theorem due to Kurtz [Stoch. Proc. Appl. 6, 223 (1978)], a Fokker-Planck equation that embodies the statistical properties of the original lattice model. The statistical equivalence of this Fokker-Planck equation, solved in terms of the associated Langevin equation, and solutions of the Chapman-Kolmogorov equation, as determined by kinetic Monte Carlo (KMC) simulations of the lattice transition rules, is demonstrated by comparing the surface roughness and the lateral height correlations obtained from the two formulations for the Edwards-Wilkinson [Proc. R. Soc. London Ser. A 381, 17 (1982)] and Wolf-Villain [Europhys. Lett. 13, 389 (1990)] deposition models, and for a model with random deposition and surface diffusion. In each case, as the largeness parameter is increased, the Langevin equation converges to the surface roughness and lateral height correlations produced by KMC simulations for all times, including the crossover between different scaling regimes. We conclude by examining some of the wider implications of these results, including applications to heteroepitaxial systems and the passage to the continuum limit.
Collapse
Affiliation(s)
- Alvin L-S Chua
- The Blackett Laboratory, Imperial College, London SW7 2BW, United Kingdom
| | | | | | | |
Collapse
|
16
|
Aarão Reis FDA. Numerical study of discrete models in the class of the nonlinear molecular beam epitaxy equation. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2004; 70:031607. [PMID: 15524534 DOI: 10.1103/physreve.70.031607] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/15/2004] [Indexed: 05/24/2023]
Abstract
We study numerically some discrete growth models belonging to the class of the nonlinear molecular beam epitaxy equation, or the Villain-Lai-Das Sarma (VLDS) equation. The conserved restricted solid-on-solid model (CRSOS) with maximum height differences Delta H(max)=1 and Delta H(max)=2 was analyzed in substrate dimensions d=1 and d=2 . The Das Sarma and Tamborenea (DT) model and a competitive model involving random deposition and CRSOS deposition were studied in d=1. For the CRSOS model with Delta H(max)=1, we obtain the more accurate estimates of scaling exponents in d=1:roughness exponent alpha=0.94+/-0.02 and dynamical exponent z=2.88+/-0.04. These estimates are significantly below the values of one-loop renormalization for the VLDS theory, which confirms Janssen's proposal of the existence of higher-order corrections. The roughness exponent in d=2 is very near the one-loop result alpha=2/3, in agreement with previous works. The moments W(n) of orders n=2 , 3, 4 of the height distribution were calculated for all models, and the skewness S triple bond W3/W(3/2)(2) and the kurtosis Q triple bond W4/W(2)2-3 were estimated. At the steady states, the CRSOS models and the competitive model have nearly the same values of S and Q in d=1, which suggests that these amplitude ratios are universal in the VLDS class. The estimates for the DT model are different, possibly due to their typically long crossover to asymptotic values. Results for the CRSOS models in d=2 also suggest that those quantities are universal.
Collapse
Affiliation(s)
- F D A Aarão Reis
- Instituto de Física, Universidade Federal Fluminense, Avenida Litorânea s/n, 24210-340 Niterói RJ, Brazil
| |
Collapse
|
17
|
Castez MF, Salvarezza RC, Solari HG. Probing universality classes in solid-on-solid deposition. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2004; 70:011605. [PMID: 15324063 DOI: 10.1103/physreve.70.011605] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/27/2004] [Indexed: 05/24/2023]
Abstract
We consider several stochastic processes corresponding to the same physical solid-on-solid deposition problem. Simplified models presenting the same (conditional) mean and variance for each process are also introduced as well as generalizations in terms of the deposition of blobs and probabilistic deposition rules. We compare the evolution of the roughness as a function of time for a three-parameter family that includes as limit cases the Family model and the Edwards-Wilkinson equation, showing that in all cases the derived models with the same mean and variance are indistinguishable from the originating models in terms of the evolution of the roughness. Finally, we show that although all the models studied belong to the same universality class, some relevant features such as the final surface roughness are reproduced only for models within a restricted class determined by sharing the same (conditional) mean and variance.
Collapse
Affiliation(s)
- Marcos F Castez
- Instituto de Investigaciones Fisicoquímicas Teóricas y Aplicadas (INIFTA), UNLP, CONICET, Casilla de Correo 16, Sucursal 4, (1900) La Plata, Argentina
| | | | | |
Collapse
|
18
|
Constantin M, Dasgupta C, Chatraphorn PP, Majumdar SN, Sarma SD. Persistence in nonequilibrium surface growth. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2004; 69:061608. [PMID: 15244586 DOI: 10.1103/physreve.69.061608] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/22/2004] [Revised: 03/29/2004] [Indexed: 05/24/2023]
Abstract
Persistence probabilities of the interface height in ( 1+1 ) - and ( 2+1 ) -dimensional atomistic, solid-on-solid, stochastic models of surface growth are studied using kinetic Monte Carlo simulations, with emphasis on models that belong to the molecular beam epitaxy (MBE) universality class. Both the initial transient and the long-time steady-state regimes are investigated. We show that for growth models in the MBE universality class, the nonlinearity of the underlying dynamical equation is clearly reflected in the difference between the measured values of the positive and negative persistence exponents in both transient and steady-state regimes. For the MBE universality class, the positive and negative persistence exponents in the steady-state are found to be theta(S)(+) =0.66+/-0.02 and theta(S)(-) =0.78+/-0.02, respectively, in ( 1+1 ) dimensions, and theta(S)(+) =0.76+/-0.02 and theta(S)(-) =0.85+/-0.02, respectively, in ( 2+1 ) dimensions. The noise reduction technique is applied on some of the ( 1+1 ) -dimensional models in order to obtain accurate values of the persistence exponents. We show analytically that a relation between the steady-state persistence exponent and the dynamic growth exponent, found earlier to be valid for linear models, should be satisfied by the smaller of the two steady-state persistence exponents in the nonlinear models. Our numerical results for the persistence exponents are consistent with this prediction. We also find that the steady-state persistence exponents can be obtained from simulations over times that are much shorter than that required for the interface to reach the steady state. The dependence of the persistence probability on the system size and the sampling time is shown to be described by a simple scaling form.
Collapse
Affiliation(s)
- M Constantin
- Condensed Matter Theory Center, Department of Physics, University of Maryland, College Park, Maryland 20742-4111, USA
| | | | | | | | | |
Collapse
|
19
|
Vvedensky DD. Crossover and universality in the Wolf-Villain model. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2003; 68:010601. [PMID: 12935119 DOI: 10.1103/physreve.68.010601] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/07/2003] [Indexed: 05/24/2023]
Abstract
The transition rules of the Wolf-Villain model for the deposition and instantaneous relaxation of particles on a lattice are expressed as a Langevin equation for the height fluctuations at each site. A coarse-graining transformation of this equation shows directly that this model belongs to the Edwards-Wilkinson universality class, in agreement with kinetic Monte Carlo simulations. The crossover from the Mullins-Herring equation is explained by the transformation under coarse graining of the coefficients in the equation of motion.
Collapse
|
20
|
Chatraphorn PP, Sarma SD. Layer-by-layer epitaxy in limited mobility nonequilibrium models of surface growth. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2002; 66:041601. [PMID: 12443210 DOI: 10.1103/physreve.66.041601] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/04/2002] [Indexed: 05/24/2023]
Abstract
We study, using noise-reduction techniques, layer-by-layer epitaxial growth in limited mobility solid-on-solid nonequilibrium surface growth models, which have been introduced in the context of kinetic surface roughening in ideal molecular beam epitaxy. Multiple hit noise reduction and long surface diffusion length lead to qualitatively similar layer-by-layer epitaxy in (1+1)- and (2+1)-dimensional limited mobility growth simulations. We discuss the dynamic scaling characteristics connecting the transient layer-by-layer growth regime with the asymptotic kinetically rough growth regime.
Collapse
Affiliation(s)
- P Punyindu Chatraphorn
- Condensed Matter Theory Center, Department of Physics, University of Maryland, College Park, MD 20742-4111, USA
| | | |
Collapse
|