Damage-spreading dynamic scaling for the ising model on the sierpinski gasket fractal

We study the relaxation towards equilibrium of the ferromagnetic Ising model on the Sierpinski gasket, which is a fractal lattice. We do this by performing Monte Carlo simulations, based on the heat-bath dynamics, and investigating the time evolution of the Hamming distance between two different configurations of the model. Starting with an initial damage created in all lattice sites, we calculate the average values of two quantities that characterize the relaxation process: the nonlinear damage relaxation time (tau), and the time for all sites to be undamaged at least once (tau(c)). We find that tau diverges, at low temperatures, with a dynamical exponent z which depends linearly on the inverse of temperature, as predicted by a generalized scaling theory developed by Henley. There is a complete breakdown of scaling for tau(c).     

自对耦无序分布随机链Potts模型的临界普适性研究

以蒙特卡罗模拟方法对自对耦分布二维随机链q态Potts模型的短时临界行为进行了数值研究.利用初始非平衡演化阶段存在的普适幂指数和有限体积标度行为，数值模拟了在不同形式随机分布时q＝3和q＝8态Potts模型磁临界指数η和动力学临界指数z.计算结果发现η不依赖于自对偶无序分布的具体形式， 从而以数值方法给出了一个关于淬火掺杂自旋系统的临界普适行为的验证. 关键词： 随机链Potts模型 动力学蒙特卡罗模拟 临界普适性     

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.     

Self-organized criticality in a block lattice model of the brittle crust

An earthquake model is introduced, in which the brittle crust is treated as a two-dimensional system of many blocks divided by faults, and the mechanical behavior of the faults is described by the Burridge-Knopoff stick-slip model. The coherent system naturally evolves into a self-organized critical state. Some universal scaling laws of seismicity, such as the Gutenberg-Richter law with the b value in agreement with the observational result and the fractal feature of fault patterns, are reproduced. Some ambiguity in simple cellular automata models is also solved.     

The Scaling Exponent for 2-Dimensional 6-State Potts Model

Using an optimal Monte Carlo renormalization group method, the scaling exponent for six-state Potts model on 2-dimensional random triangle lattice ie studied. The scaling exponent υ=0.54, is consistent with the expectant value of the scaling theory.     

Electronic shot noise in fractal conductors

By solving a master equation in the Sierpiński lattice and in a planar random-resistor network, we determine the scaling with size L of the shot noise power P due to elastic scattering in a fractal conductor. We find a power-law scaling P proportional, variantL;{d_{f}-2-alpha}, with an exponent depending on the fractal dimension d_{f} and the anomalous diffusion exponent alpha. This is the same scaling as the time-averaged current I[over ], which implies that the Fano factor F=P/2eI[over ] is scale-independent. We obtain a value of F=1/3 for anomalous diffusion that is the same as for normal diffusion, even if there is no smallest length scale below which the normal diffusion equation holds. The fact that F remains fixed at 1/3 as one crosses the percolation threshold in a random-resistor network may explain recent measurements of a doping-independent Fano factor in a graphene flake.     

Dynamics of the formation of two-dimensional ordered structures

The growth of ordered domains in lattice gas models, which occurs after the system is quenched from infinite temperature to a state below the critical temperatureT _c, is studied by Monte Carlo simulation. For a square lattice with repulsion between nearest and next-nearest neighbors, which in equilibrium exhibits fourfold degenerate (2×1) superstructures, the time-dependent energy E(t), domain size L(t), and structure functionS(q, t) are obtained, both for Glauber dynamics (no conservation law) and the case with conserved density (Kawasaki dynamics). At late times the energy excess and halfwidth of the structure factor decrease proportional tot ^–x, whileL(t) t ^x, where the exponent x=1/2 for Glauber dynamics and x1/3 for Kawasaki dynamics. In addition, the structure factor satisfies a scaling lawS(k,t)=t ^2xS(kt^x). The smaller exponent for the conserved density case is traced back to the excess density contained in the walls between ordered domains which must be redistributed during growth. Quenches toT>T _c, T=T_c (where we estimate dynamic critical exponents) andT=0 are also considered. In the latter case, the system becomes frozen in a glasslike domain pattern far from equilibrium when using Kawasaki dynamics. The generalization of our results to other lattices and structures also is briefly discussed.     

New universal short-time scaling behaviour of critical relaxation processes

We study the critical relaxation properties of Model A (purely dissipative relaxation) starting from a macroscopically prepared initial state characterised by non-equilibrium values for order parameter and correlations. Using a renormalisation group approach we observe that even (macroscopically)early stages of the relaxation process display universal behaviour governed by a new, independent initial slip exponent. For large times, the system crosses over to the well-known long-time relaxation behaviour.The new exponent is calculated toO(²) in =4–d, whered is the spatial dimension of the system. The initial slip scaling form of general correlation and response functions as well as the order parameter is derived, exploiting a short-time operator expansion. The leading scaling behaviour is determined by initial states with sharp values of the order parameter. Non-vanishing correlations generate corrections to scaling.     

Velocity intermittency in a buoyancy subrange in a two-dimensional soap film convection experiment

In a two-dimensional soap film convection experiment, the velocity fields are found to be strongly intermittent in the buoyancy subrange, which was reported to be nonintermittent in a recent numerical simulation. The structure functions Sq(l)(= ) exhibit self-similar structures and can be described by power laws l(zetaq) for integers 8 > or = q < or = 1. By extending Kolmogorov's refined similarity hypothesis to our system, an analytical form is derived for the scaling exponent zeta(q) = q/2 + (mu/18)(3q - q2). Our measurements yield mu = 0.42, which is significantly greater than 0.2 found in high Reynolds number turbulence in wind tunnels.     

Scaling of particle trajectories on a lattice

The scaling behavior of the closed trajectories of a moving particle generated by randomly placed rotators or mirrors on a square or triangular lattice is studied numerically. On both lattices, for most concentrations of the scatterers the trajectories close exponentially fast. For special critical concentrations infinitely extended trajectories can occur which exhibit a scaling behavior similar to that of the perimeters of percolation clusters.At criticality, in addition to the two critical exponents =15/7 andd _f=7/4 found before, the critical exponent =3/7 appears. This exponent determines structural scaling properties of closed trajectories of finite size when they approach infinity. New scaling behavior was found for the square lattice partially occupied by rotators, indicating a different universality class than that of percolation clusters.Near criticality, in the critical region, two scaling functions were determined numerically:f(x), related to the trajectory length (S) distributionn _s, andh(x), related to the trajectory sizeR _s (gyration radius) distribution, respectively. The scaling functionf(x) is in most cases found to be a symmetric double Gaussian with the same characteristic size exponent =0.433/7 as at criticality, leading to a stretched exponential dependence ofn _S onS, n_Sexp(–S ^6/7). However, for the rotator model on the partially occupied square lattice an alternative scaling function is found, leading to a new exponent =1.6±0.3 and a superexponential dependence ofn _S onS.h(x) is essentially a constant, which depends on the type of lattice and the concentration of the scatterers. The appearance of the same exponent =3/7 at and near a critical point is discussed.     

Footprints of nonextensive Tsallis statistics, selfaffinity and universality in the preparation of the L’Aquila earthquake hidden in a pre-seismic EM emission

Fracture induced physical fields allow a real-time monitoring of damage evolution in materials during mechanical loading. We investigate the preparation of the recently occurred L’Aquila earthquake in terms of a detected precursory electromagnetic anomaly. The precursor is well described by a recently introduced model for earthquake dynamics, which has been rooted in a nonextensive Tsallis framework starting from first principles. The analysis in terms of nonextensivity implies that the well established aspect of self-affine nature of faulting and fracture is hidden into the precursor, namely, the activation of the L’Aquila fault is a reduced self-affine image of the regional seismicity covering many geological faults. The Gutenberg-Richter magnitude-frequency relationship, the best known scaling relation for earthquakes, verifies the results based on nonextensivity. The latter suggests that the activation of the L’Aquila fault is a magnified image of the laboratory seismicity by means of acoustic and electromagnetic emissions. Finally, we present evidence for universality in magnetic storm, earthquake, and electromagnetic precursor occurrence by means of complexity and nonextensivity.     

Critical Behaviors and Finite-Size Scaling of Principal Fluctuation Modes in Complex Systems

Complex systems consisting of N agents can be investigated from the aspect of principal fluctuation modes of agents. From the correlations between agents, an N×N correlation matrix C can be obtained. The principal fluctuation modes are defined by the eigenvectors of C. Near the critical point of a complex system, we anticipate that the principal fluctuation modes have the critical behaviors similar to that of the susceptibity. With the Ising model on a two-dimensional square lattice as an example, the critical behaviors of principal fluctuation modes have been studied. The eigenvalues of the first 9 principal fluctuation modes have been invesitigated. Our Monte Carlo data demonstrate that these eigenvalues of the system with size L and the reduced temperature t follow a finite-size scaling form λ_n(L, t)=L^γ/ν f_n(tL^1/ν), where γ is critical exponent of susceptibility and ν is the critical exponent of the correlation length. Using eigenvalues λ₁, λ₂ and λ₆, we get the finite-size scaling form of the second moment correlation length ξ(L, t)=Lξ(tL^1/ν). It is shown that the second moment correlation length in the two-dimensional square lattice is anisotropic.     

The bond bending model on triangular lattice and square lattice

The bond bending model is studied using the series expansion method on a triangular lattice and on a square lattice. The elastic splay susceptibility χ_SR and the elastic compressional susceptibility χ_el are calculated up to 11th order for the triangular lattice and up to 14th order for the square lattice. The elastic splay crossover exponent, ζ_SP, is found to be ζ_SP ≈ 1.26 ± 0.05 for the triangular lattice and ζ_SP = 1.30 ± 0.04 for the square lattice which is close to the conductivity exponent, ζ_Re, of the resistor network. From the scaling relation ? _B = dv + ζ_SP, we found that the bulk modulus exponent ? _B = 3.93 ± 0.05 for the triangular lattice and ? _B = 3.97 ± 0.04 for the square lattice which is in good agreement with the result ? _B = 3.96 ± 0.04, obtained by Zabolitzky et al. using a transfer matrix technique on a honeycomb lattice.     

辐照损伤对钨晶格热导率影响的分子动力学研究

钨是最具应用前景的面向等离子体候选材料,但核聚变堆内强烈的辐照环境会使钨的近表面区域产生辐照损伤,进而影响其关键的导热性能.本文构建了包含辐照损伤相关缺陷的晶体钨模型,并采用非平衡分子动力学的方法定量研究了这些缺陷对钨导热性能的影响.结果表明,随中子辐射能量的增加,晶体内部留下的Frenkel缺陷数目增多进而导致钨的晶格热导率降低;间隙原子比空位更易于向晶界偏聚,且钨中的间隙钨原子与空位相比,使晶格热导率下降程度更大.纳米级氦气泡导致晶格热导率的显著降低,气孔率为2.1%时晶格热导率降至完美晶体的约25%.这些不同的缺陷造成不同程度的周围晶格扭曲,增加了声子散射几率,是导致晶格热导率下降的根源.     

Lyapunov exponents and anomalous diffusion of a Lorentz gas with infinite horizon using approximate zeta functions

We compute the Lyapunov exponent, the generalized Lyapunov exponents, and the diffusion constant for a Lorentz gas on a square lattice, thus having infinite horizon. Approximate zeta functions, written in terms of probabilities rather than periodic orbits, are used in order to avoid the convergence problems of cycle expansions. The emphasis is on the relation between the analytic structure of the zeta function, where a branch cut plays an important role, and the asymptotic dynamics of the system. The Lyapunov exponent for the corresponding map agrees with the conjectured limit _map = -2 log(R) + C + O(R) and we derive an approximate value for the constantC in good agreement with numerical simulations. We also find a diverging diffusion constantD(t)logt and a phase transition for the generalized Lyapunov exponents.     

How well do Monte Carlo simulations distinguish between U(1) and SU(2)?: A study of the string tension in U(1) lattice gauge theory

To investigate how a system with a known deconfining phase transition behaves when studied on finite lattices via Monte Carlo simulations, we have made such studies of compact U(1) lattice gauge theory for 8⁴, 10⁴, and 12⁴ lattices. We have concentrated on the mean plaquette energy and the string tension. The string tension does not vanish on a finite lattice, but using finite size scaling arguments the indications are that it does vanish on an infinite lattice, where we predict the critical coupling β_c = 1.008 and the correlation length exponent $\text{ν =}\text{1}\text{3}$. We compare our results to those for SU(2) and find that although there are differences, they are not yet definitive.     

Investigation of the critical behavior of the three-dimensional weakly diluted Potts model

The static critical behavior of the three-dimensional weakly diluted Potts model with the state q = 3 on a simple cubic lattice has been investigated by the Monte Carlo method using the Wolff single-cluster algorithm. It is shown that at the spin concentrations p = 0.9 and 0.8 a second-order phase transition is observed in the three-dimensional weakly diluted Potts model with the state q = 3. On the basis of the finite-size scaling theory, we calculated the static critical exponents of the specific heat α, susceptibility γ, magnetization β, and the correlation-length exponent v.     

Scaling Region in Desynchronous Coupled Chaotic Maps

The largest Lyapunov exponent and the Lyapunov spectrum of a coupled map lattice are studied when the system state is desynchronous chaos. In the large system size limit a scaling region is found in the parameter space where the largest Lyapunov exponent is independent of the system size and the coupling strength. Some scaling relation between the Lyapunov spectrum distributions for different coupling strengths is found when the coupling strengths are taken in the scaling parameter region. The existence of the scaling domain and the scaling relation of Lyapunov spectra there are heuristically explained.     

Scaling Region in Desynchronous Coupled Chaotic Maps

The largest Lyapunov exponent and the Lyapunov spectrum of a coupled map lattice are studied when the system state is desynchronous chaos. In the large system size limit a scaling region is found in the parameter space where the largest Lyapunov exponent is independent of the system size and the coupling strength. Some scaling relation between the Lyapunov spectrum distributions for different coupling strengths is found when the coupling strengths are taken in the scaling parameter region. The existence of the scaling domain and the scaling relation of Lyapunov spectra there are heuristically explained.     

Superconductor-insulator transition of Josephson-junction arrays on a honeycomb lattice in a magnetic field

We study the superconductor to insulator transition at zero temperature in aJosephson-junction array model on a honeycomb lattice with f flux quantum perplaquette. The path integral representation of the model corresponds to a (2 + 1)-dimensional classical model, which isused to investigate the critical behavior by extensive Monte Carlo simulations on largesystem sizes. In contrast to the model on a square lattice, the transition is found to befirst order for f = 1 /3 and continuous for f = 1 / 2 but in a different universality class.The correlation-length critical exponent is estimated from finite-size scaling of vortexcorrelations. The estimated universal conductivity at the transition is approximately fourtimes its value for f =0. The results are compared with experimental observations on ultrathinsuperconducting films with a triangular lattice of nanoholes in a transverse magneticfield.