分形基底上刻蚀模型动力学标度行为的 数值模拟研究

为探讨分形基底结构对生长表面标度行为的影响, 本文采用Kinetic Monte Carlo(KMC)方法模拟了刻蚀模型(etching model)在谢尔宾斯基箭头和蟹状分形基底上刻蚀表面的动力学行为. 研究表明,在两种分形基底上的刻蚀模型都表现出很好的动力学标度行为, 并且满足Family-Vicsek标度规律. 虽然谢尔宾斯基箭头和蟹状分形基底的分形维数相同, 但模拟得到的标度指数却不同, 并且粗糙度指数 α与动力学指数z也不满足在欧几里得基底上成立的标度关系α+z=2. 由此可以看出, 标度指数不仅与基底的分形维数有关, 而且和分形基底的具体结构有关.     

分形基底上受限固-固模型动力学性质的数值模拟研究

为了探讨非完整基底结构对生长表面动力学行为的影响,本文在具有相同分形维数而不同谱维数的谢尔宾斯基箭头和蟹状分形基底上对受限固-固(restricted solid-on-solid,RSOS)模型的生长过程进行了大量的数值模拟研究.通过计算表面宽度和饱和表面极值高度的统计行为对生长表面的动力学行为进行了分析.结果表明,分形基底结构对生长表面的动力学行为具有显著的影响.尽管在两种基底上受限固-固模型的表面宽度均表现出很好的动力学标度行为,仍然满足Family-Vicsek标度规律,但由此计算得到的动力学标度指数并不相同.饱和生长表面的极值高度并不能满足三种常用的极值统计分布,即Weibull,Gumbel和Frechet分布,而是能很好地符合Asym2Sig分布.     

科赫分形基底上受限固-固模型动力学标度行为的数值研究

为分析基底结构对离散生长模型动力学性质的影响, 本文在随机游走指数十分接近而分形维数和谱维数均不相同的科赫格子和科赫曲线分形基底上对受限固-固(restricted solid-on-solid)模型的生长过程进行数值模拟研究. 通过分析表面宽度和饱和表面极值高度的统计行为发现: 随机游走的动力学指数能够对饱和粗化表面的动力学行为起主要贡献. 尽管分形基底具有不同的分形维数和谱维数, 但是在两种分形基底上得到了在误差范围内相同的粗造度指数. 两种分形基底上饱和表面相对生长高度极大(小)值分布分别可以很好的塌缩在一起, 且很好的满足Asym2Sig函数分布.     

含时空关联噪声生长方程的标度奇异性分析

基于动力学重整化群理论研究表面界面生长动力学标度奇异性问题, 得到含时空关联噪声的表面生长方程标度奇异指数的一般结果,并将此方法应用于几种典型的局域生长方程——Kardar-Parisi-Zhang（KPZ）方程、线性生长方程、Lai-Das Sarma-Villain（LDV）方程.结果表明,在长波长极限下局域生长方程的动力学标度奇异性与最相关项、基底维数以及噪声有关,并且若出现标度奇异性,只会是超粗化（super rough）奇异标度行为,而不是内禀（intrinsically）奇异标度行为. 关键词： 标度奇异性 动力学重整化群理论 时空关联噪声     

含关联噪声的空间分数阶随机生长方程的动力学标度行为研究

为探讨含关联噪声的空间分数阶随机生长方程的动力学标度行为,本文利用Riesz分数阶导数和Grümwald-Letnikov分数阶导数定义方法研究了关联噪声驱动下的空间分数阶Edwards-Wilkinson (SFEW)方程在1+1维情况下的数值解,得到了不同噪声关联因子和分数阶数时的生长指数、粗糙度指数、动力学指数等,所求出的临界指数均与标度分析方法的结果相符合.研究表明噪声关联因子和分数阶数均影响到SFEW方程的动力学标度行为,且表现为连续变化的普适类.     

点缺陷对表面生长动力学标度行为的影响

基于含噪声Kuramoto-Sivashinsky方程, 采用动力学重正化群技术, 研究生长表面存在点缺陷或杂质对表面生长动力学标度行为的影响, 得到了相应的粗糙度指数α 和动力学标度指数z. 所得结果表明, 点缺陷的存在使生长表面粗化, 并缩短达到稳定生长的弛豫时间.     

非局域Sun-Guo-Grant方程的自洽模耦合理论

通过选取具有正确渐近行为的标度函数形式，将自洽的模耦合理论推广应用到对非局域的Sun-Guo-Grant方程的动力学标度性质的研究中.通过分析得到，在强耦合区基底维数d=1,2的情况下，动力学指数z随非局域参数ρ的变化关系.将这一结果与动力学重正化群理论和直接标度分析得到的结果进行了对比. 关键词： 表面粗糙生长动力学 动力学标度 自洽模耦合理论     

Kuramoto-Sivashinsky与Karda-Parisi-Zhang模型中生长界面分形特性研究

在（2＋1）维情况下，利用数值模拟研究了Kuramoto-Sivashinsky (K-S)与Karda-Parisi-Z hang (KPZ)模型所决定的非平衡态界面生长演化过程.结果表明，KPZ与K-S模型都表现出明 显的时间和空间标度特性.相对于KPZ模型而言，K-S模型所对应的表面具有更明显的颗粒特 征，当生长时间较长时，生长界面呈现蜂窝状结构.通过数值相关分析得到了生长界面的粗 糙度指数、生长指数和动态标度指数等参数.从两种模型对应的表面形貌特征和表面参数来 看，在(2+1)维情况下，KPZ与K-S模型所决定的表面具有完全不同的动态标度行为，属于不 同的两类物理模型. 关键词： Kuramoto-Sivashinsky (K-S)模型 Karda-Parisi-Zhang(KPZ)模型 分形 数值模拟     

1+1维抛射沉积模型内部结构动力学行为的数值研究

采用Kinetic Monte Carlo方法对1+1维抛射沉积(BD)模型内部结构的动力学行为进行了大量的数值模拟研究.分别分析了空洞密度和内部界面的动力学行为.研究表明,空洞密度呈高斯型分布,其平均值首先随生长时间快速增长,然后达到一个与基底尺寸无关的饱和值.除表面宽度,还引入了新的极值统计方法来分析该模型内部界面的动力学行为,分析结果显示,1+1维BD模型内部界面的演化满足标准的Family-Vicsek标度规律,并且属Kardar-Parisi-Zhang方程所描述的普适类.最后对表面宽度和极值统计两种理论方法的有限尺寸效应进行了比较.     

1+1维Wolf-Villain模型奇异标度行为的数值模拟

采用Kinetic Monte Carlo(KMC)方法对描述分子束外延生长(MBE)的1+1维Wolf-Villain模型进行大尺寸和长生长时间的数值模拟研究,以消除渡越行为的影响.计算得到整体和局域标度指数.结果显示,在所模拟的空间和时间尺度范围内,1+1维Wolf-Villain模型仍呈现出固有奇异标度行为.这一结论与López等人最近的理论分析结果不一致.     

Influences of Fractal Substrate Structures on Dynamic Scaling Behaviors of Etching Model

The Etching model on various fractal substrates embedded in two dimensions was investigated by means of kinetic Mento Carlo method in order to determine the relationship between dynamic scaling exponents and fractal parameters. The fractal dimensions are from 1.465 to 1.893, and the random walk exponents are from 2.101 to 2.578.It is found that the dynamic behaviors on fractal lattices are more complex than those on integer dimensions. The roughness exponent increases with the increasing of the random walk exponent on the fractal substrates but shows a non-monotonic relation with respect to the fractal dimension. No monotonic change is observed in the growth exponent.     

Critical behavior of a two-species reaction-diffusion problem in 2D

In this work we study the critical behavior of a model that simulates the propagation of an epidemic process over a population. We simulate the model on two-dimensional finite lattices to determine the critical density of the diffusive population. A finite size scaling analysis is employed to determine the order parameter and correlation length critical exponents.     

Scaling of waves in the bak-tang-wiesenfeld sandpile model

We study probability distributions of waves of topplings in the Bak-Tang-Wiesenfeld model on hypercubic lattices for dimensions D>/=2. Waves represent relaxation processes which do not contain multiple toppling events. We investigate bulk and boundary waves by means of their correspondence to spanning trees, and by extensive numerical simulations. While the scaling behavior of avalanches is complex and usually not governed by simple scaling laws, we show that the probability distributions for waves display clear power-law asymptotic behavior in perfect agreement with the analytical predictions. Critical exponents are obtained for the distributions of radius, area, and duration of bulk and boundary waves. Relations between them and fractal dimensions of waves are derived. We confirm that the upper critical dimension D(u) of the model is 4, and calculate logarithmic corrections to the scaling behavior of waves in D=4. In addition, we present analytical estimates for bulk avalanches in dimensions D>/=4 and simulation data for avalanches in D相似文献   

Dynamic scaling behaviors of the restricted-solid-on-solid model on honeycomb and square-octagon lattice substrates

The dynamic scaling behaviors of the restricted-solid-on-solid (RSOS) model on two new types of substrate, which are honeycomb and square-octagon lattice substrates, are studied by means of Kinetic Monte Carlo simulations. The growth exponent β and the roughness exponent α defined, respectively, by the surface width via W ~ t ^β and the saturated width via W _sat ~ L ^α , L being the system size, were obtained by a power-counting analysis. Our simulation results show that the Family-Vicsek scaling is still satisfied. However, the structures of the substrates indeed affect the dynamic behavior of the growth model. The values of the roughness exponents fall between regular and fractal lattices. Deeper analysis show that the coordination number of the substrates play an crucial role.     

Critical behavior of a cubic-lattice 3D Ising model for systems with quenched disorder

A Monte Carlo method is applied to simulate the static critical behavior of a cubic-lattice 3D Ising model for systems with quenched disorder. Numerical results are presented for the spin concentrations of p = 1.0, 0.95, 0.9, 0.8, 0.6 on L × L × L lattices with L = 20–60 under periodic boundary conditions. The critical temperature is determined by the Binder cumulant method. A finite-size scaling technique is used to calculate the static critical exponents α, β, γ, and ν (for specific heat, susceptibility, magnetization, and correlation length, respectively) in the range of p under study. Universality classes of critical behavior are discussed for three-dimensional diluted systems.     

Computer simulation of the critical behavior of 3D disordered ising model

The critical behavior of the disordered ferromagnetic Ising model is studied numerically by the Monte Carlo method in a wide range of variation of concentration of nonmagnetic impurity atoms. The temperature dependences of correlation length and magnetic susceptibility are determined for samples with various spin concentrations and various linear sizes. The finite-size scaling technique is used for obtaining scaling functions for these quantities, which exhibit a universal behavior in the critical region; the critical temperatures and static critical exponents are also determined using scaling corrections. On the basis of variation of the scaling functions and values of critical exponents upon a change in the concentration, the conclusion is drawn concerning the existence of two universal classes of the critical behavior of the diluted Ising model with different characteristics for weakly and strongly disordered systems.     

Percolation on patchwise heterogeneous surfaces under equilibrium conditions

The site percolation problem on square lattices whose sites are grouped in two types of energetically different patches is studied. Several lattices formed by collections of either randomly or orderly localized and no overlapped patches of different sizes are generated. The system is characterized by two parameters, namely, the size of each patch, l, and the energy difference between the two kind of sites, ΔE. Particles are adsorbed at equilibrium on the lattice. The critical coverage is determined by means of Monte Carlo simulations and finite-size scaling analysis. The percolative behavior of the system as a function of the parameters characterizing the heterogeneity of the energetic surface topography is presented and discussed.