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.     

Finite-size scaling properties of the damage distance and dynamical critical exponent for the Ising model

With the damage spreading method, scaling properties of the damage distance on the Ising model with heat bath dynamics are studied numerically. With the parallel flipping scheme, the scaling curves fall on two curves, which depend on the odd or even lattice sizes. The both scaling curves give the consistent dynamical exponent as z = 2.16±0.04 for d = 2 and z = 2.09±0.05 for d = 3, respectively. By shifting one of them, two curves overlap each other perfectly. Meanwhile, all the scaling curves obtained by single-spin flipping processes (with different odd or even lattice sizes) fall on a single curve, from which the consistent dynamical critical exponent with the parallel scheme is obtained z = 2.18±0.02 for d = 2 and z = 2.08±0.04 for d = 3.     

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

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

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

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

三角格点基底上磁性分形团簇形貌演化规律

在扩散限制凝聚模型基础上引入粒子的自旋自由度,将磁耦合系数扩展为随自旋间距离幂次变化的非常数项J/r^a,采用Monte Carlo方法研究在二维三角格点基底上具有幂次相互作用的磁性团簇形貌及其分形维数D_f的演化规律.模拟结果表明,对于较大的幂指数α值,即α≥5时,团簇形貌随耦合参数J的变化较小,其分形维数D_f在1.50~1.70之间;随着α值的减小,团簇形貌随参数J有一明显的演化过程,在模拟范围内,分形维数D_f在1.20~1.90之间.     

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

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

Short-time dynamic behaviour of critical XY systems

Using Monte Carlo methods, the short-time dynamic scaling behaviour of two-dimensional critical XY systems is investigated. Our results for the XY model show that there exists universal scaling behaviour already in the short-time regime, but the values of the dynamic exponent z differ for different initial conditions. For the fully frustrated XY model, power law scaling behaviour is also observed in the short-time regime. However, a violation of the standard scaling relation between the exponents is detected.     

Hyperscaling for polymer rings

The statistics of a long closed self-avoiding walk (SAW) or polymer ring on a d-dimensional lattice obeys hyperscaling. The combination p_NR²_N^d/2μ^−N (where p_N is the number of configurations of an oriented and rooted N-step ring, R²_N a typical average size squared, and μ the SAW effective connectivity constant of the lattice) is equal for N å ∞ to a lattice-dependent constant times a universal amplitude A(d). The latter amplitude is calculated directly from the minimal continous Edwards model to second order in 4 − d. The case of rings at the upper critical dimension d = 4 is also studied. The results are checked against field-theoretical calculations, and former simulations. As a consequence, we show that the universal constant λ appearing to second order in in all critical phenomena amplitude ratios is equal to .     

多孔材料分形扩散模型的Fourier-Bessel级数算法及其应用

将Fick扩散定律的Fourier三角级数算法推广成多孔材料分形扩散模型的Fourier-Bessel级数算法,并把它应用于化学工程中吸附问题涉及的浓度分布与相对吸附量的计算中,取得一些规律性认识.由于分形扩散模型是在Fick扩散定律的基础上增加了表征微观结构的参数d_f和θ,研究多孔材料中的浓度分布与相对吸附量时,与Fick扩散定律的研究结果相比,定性上基本一致,在定量上有差别,d_f和θ对扩散传质过程的影响各有侧重,用它们可更好地描述多孔材料中的扩散过程.     

Transport properties of saturated and unsaturated porous fractal materials

Inviscid, irrotational flow through fractal porous materials is studied. The key parameter is the variation of tortuosity with the filling fraction phi of fluid in the porous material. Altering the filling fraction provides a way of probing the effect of the fractal structure over all its length scales. The variation of tortuosity with phi is found to follow a power law of the form alpha approximately phi (-E) for deterministic and stochastic fractals in two and three dimensions. A phenomenological argument for the scaling of tortuosity alpha with filling fraction phi is presented and is given by alpha approximately phi(D_{w}-2/D_{f}-d_{E}), where D_{f} is the fractal dimension, D_{w} is the random walk dimension, and d_{E} is the Euclidean dimension. Numerically calculated values of the exponents show good agreement with those predicted from the phenomenological argument for both the saturated and the unsaturated model.     

Binding energy of an exciton bound to ionized donors in quantum dots

Binding energies for an exciton (X) trapped in the two-dimensional quantum dot by a positive ion located on the z axis at a distance d from the dot plane are calculated by using the method of few-body physics. This configuration is called a barrier (D⁺,X) center. The dependence of the binding energy of the ground state of the barrier (D⁺,X) center on the dot radius for a few values of the distance d between the fixed positive ion on the z axis and the dot plane is obtained. We find that when d<0.2nm the barrier (D⁺,X) center does not form a bound state.     

含随机裂纹网络孔隙材料渗透率的逾渗模型研究

采用逾渗理论对含随机裂纹网络的孔隙材料渗透性进行研究. 开裂孔隙材料渗透率的影响因素包括裂纹网络的几何特征、孔隙材料本体渗透率以及裂纹开度, 本文使用连续区逾渗理论模型建立了渗透率的标度律. 对于裂纹网络的几何特征, 本文基于连续区逾渗理论并考虑裂纹网络的分形特征提出了有限区域内二维随机裂纹网络的连通度定义; 对随机裂纹网络的几何分析表明, 随机裂纹局部团簇效应会降低裂纹网络的整体连通性, 随机裂纹网络的标度指数并非经典逾渗理论给出的固定值, 而是随着网络的分形维数的减小而增大. 本文在网络连通度和主裂纹团的曲折度的基础上, 提出了开裂孔隙材料渗透率标度律的解析表达, K=K₀(K_m,b)(ρ-ρ_c)^μ, 分别考虑了裂纹网络的几何逾渗特征 (ρ-ρ_c)^μ、孔隙材料渗透率K_m 以及裂纹开度比b; 对有限区域含有随机裂纹网络的孔隙材料渗透过程的有限元模拟表明, K₀ 在裂纹逾渗阈值附近与b呈指数关系, 但当裂纹的局部渗透率与K_m比值高于10⁶ 后, 开度比b对渗透率不再有影响.     

Scaling properties of eden clusters in three and four dimensions

We study the scaling properties of noise reduced Eden clusters in three and four dimensions for variant B in the strip geometry. We find that the width W for large times behaves as a(s)g(L/s^d−1), where L is the width of the strip, s the noise reduction parameter, d the dimension of space, and a(s) a decreasing function of s, g is a scaling function with the property g(u)→1/2 as u→0 and g(u)u^x as u→∞, where χ is the roughness exponent. This scaling result leads to a new way of determining χ. In 3 dimensions, our numerical values for χ support a recent conjecture by Kim and Kosterlitz: χ = 2/(d + 2), and contradict all the former analytical conjectures. In 4 dimensions, we cannot distinguish between the conjectures of Kim and Kosterlitz and the conjecture of Wolf and Kertész, because large crossovers and finite size effects make the measurement of the exponents difficult.     

Vortex phases and dynamics in high-temperature and conventional amorphous superconductors

The vortex phases and dynamics of high-temperature and conventional amorphous superconductors are investigated over a broad frequency range. A second-order vortex-solid melting transition in both YBa₂Cu₃O₇ single crystals with random point defects and amorphous Mo₃Si films is manifested by a set of universal critical exponents (ν≈2/3 and z≈3) from both ac and dc transport as well as ac magnetic susceptibility measurements. In contrast, the melting transition in YBa₂Cu₃O₇ single crystals with correlated columnar defects is found to be consistent with a Bose-glass phase transition characterized by a different set of critical exponents (ν≈1/2 and z′≈4).     

Phase diagram of optimal paths

We show that by choosing appropriate distributions of the randomness the search for optimal paths links diverse problems of disordered media, such as directed percolation, invasion percolation, and directed and nondirected spanning polymers. We also introduce a simple and efficient algorithm, which solves the d-dimensional model numerically in O(N(1+df/d)) steps, where df is the fractal dimension of the path. Using extensive simulations in two dimensions, we identify the phase boundaries of the directed polymer universality class. A new strong-disorder phase occurs where the optimum paths are self-affine with parameter-dependent scaling exponents. Furthermore, the phase diagram contains directed and nondirected percolation as well as the directed random walk models at specific points and lines.     

Simulations of migration, fragmentation and coalescence of non-wetting fluids in porous media

A model based on invasion percolation was used to simulate the migration on a non-wetting fluid through a porous medium filled with an immiscible wetting fluid under the influence of a gradient such as that provided by gravity. The migrating fluid clusters undergo both fragmentation and coalescence. The fragment size distribution obtained from two-dimensional simulations in which the gradient g is slowly increased from 0 can be represented by the scaling form N_s(g)s^-2ƒ(s|g|^-z where z=1+(D−1)ν(ν+1). Here D is the fractal dimensionality of invasion percolation, with trapping, and ν is the ordinary percolation correlation length exponent.     

Oriented polymers: A transfer matrix calculation

Based on transfer matrix techniques and finite-size scaling, we study the oriented polymer (self-avoiding walk) with nearest neighbor interaction. In the repulsive regime, various critical exponents are computed and compared with exact values predicted recently. The polymer is also found to undergo a spiral transition for sufficiently strong attractive interaction. The fractal dimension of the polymer is computed in the repulsive and attractive regimes and at the spiral transition point. The later is found to be different from that at the collapse transition of the ordinary self-avoiding walk.     

Critical exponents of the Ising model on low-dimensional fractal media

The critical behavior of the Ising model on fractal substrates with noninteger Hausdorff dimension d_H<2 and infinite ramification order is studied by means of the short-time critical dynamic scaling approach. Our determinations of the critical temperatures and critical exponents β, γ, and ν are compared to the predictions of the Wilson-Fisher expansion, the Wallace-Zia expansion, the transfer matrix method, and more recent Monte Carlo simulations using finite-size scaling analysis. We also determined the effective dimension (d_ef), which plays the role of the Euclidean dimension in the formulation of the dynamic scaling and in the hyperscaling relationship d_ef=2β/ν+γ/ν. Furthermore, we obtained the dynamic exponent z of the nonequilibrium correlation length and the exponent θ that governs the initial increase of the magnetization. Our results are consistent with the convergence of the lower-critical dimension towards d=1 for fractal substrates and suggest that the Hausdorff dimension may be different from the effective dimension.     

Exact scaling for the mean first-passage time of random walks on a generalized Koch network with a trap

In this paper, we study the scaling for the mean first-passage time (MFPT) of the random walks on a generalized Koch network with a trap. Through the network construction, where the initial state is transformed from a triangle to a polygon, we obtain the exact scaling for the MFPT. We show that the MFPT grows linearly with the number of nodes and the dimensions of the polygon in the large limit of the network order. In addition, we determine the exponents of scaling efficiency characterizing the random walks. Our results are the generalizations of those derived for the Koch network, which shed light on the analysis of random walks over various fractal networks.