Kinetics of Infection-Driven Growth Model with Birth and Death

We propose a two-species infection model, in which an infected aggregate can gain one monomer from a healthy one due to infection when they meet together. Moreover, both the healthy and infected aggregates may lose one monomer because of self-death, but a healthy aggregate can spontaneously yield a new monomer. Consider a simple system in which the birth/death rates are directly proportional to the aggregate size, namely, the birth and death rates of the healthy aggregate of size k are J1 k and J2k while the self-death rate of the infected aggregate of size k is J3k. We then investigate the kinetics of such a system by means of rate equation approach. For the J1 〉 J2 case, the aggregate size distribution of either species approaches the generalized scaling form and the typical size of either species increases wavily at large times. For the J1 = J2 case, the size distribution of healthy aggregates approaches the generalized scaling form while that of infected aggregates satisfies the modified scaling form. For the J1 〈 J2 case, the size distribution of healthy aggregates satisfies the modified scaling form, but that of infected aggregates does not scale.     

Competition Between Self-birth and Catalyzed Death in Aggregation Growth with Catalysis Injection

We propose two irreversible aggregation growth models of aggregates of two distinct species （.4 and B） to study the interactions between virus aggregates and medicine efficacy aggregates in the virus-medicine cooperative evolution system. The A-species aggregates evolve driven by self monomer birth and B-species aggregate-catalyzed monomer death in model I and by self birth, catalyzed death, and self monomer exchange reactions in model II, while the catalyst B-species aggregates are assumed to be injected into the system sustainedly or at a periodic time-dependent rate. The kinetic behaviors of the A-species aggregates are investigated by the rate equation approach based on the mean-field theory with the self birth rate kernel IA（k） = Ik, catalyzed death rate kernel JAB（k） = Jk and self exchange rate kernel KA （k, l） = Kkl. The kinetic behaviors of the A-species aggregates are mainly dominated by the competition between the two effects of the self birth （with the effective rate I） and the catalyzed death （with the effective rate JB0）, while the effects of the self exchanges of the A-species aggregates which appear in an effective rate KAo play important roles in the cases of I 〉 JBo and I = JBo. The evolution behaviors of the total mass M1^A（t） and the total aggregate number MA（t） are obtained, and the aggregate size distribution ak（t） of species A is found to approach a generalized scaling form in the case of I ≥ JBo and a special modified scaling form in the case of I 〈 JB0. The periodical evolution of the B-monomers concentration plays an exponential form of the periodic modulation.     

Competing role of catalysis-coagulation and catalysis-fragmentation in kinetic aggregation behaviours

We propose a kinetic aggregation model where species A aggregates evolve by the catalysis-coagulation and the catalysis-fragmentation,while the catalyst aggregates of the same species B or C perform self-coagulation processes.By means of the generalized Smoluchowski rate equation based on the mean-field assumption,we study the kinetic behaviours of the system with the catalysis-coagulation rate kernel K(i,j;l) ∝ l ν and the catalysis-fragmentation rate kernel F(i,j;l) ∝ l μ,where l is the size of the catalyst aggregate,and ν and μ are two parameters reflecting the dependence of the catalysis reaction on the size of the catalyst aggregate.The relation between the values of parameters ν and μ reflects the competing roles between the two catalysis processes in the kinetic evolution of species A.It is found that the competing roles of the catalysis-coagulation and catalysis-fragmentation in the kinetic aggregation behaviours are not determined simply by the relation between the two parameters ν and μ,but also depend on the values of these two parameters.When ν μ and ν≥ 0,the kinetic evolution of species A is dominated by the catalysis-coagulation and its aggregate size distribution a k(t) obeys the conventional or generalized scaling law;when ν μ and ν≥ 0 or ν 0 but μ≥ 0,the catalysis-fragmentation process may play a dominating role and a k(t) approaches the scale-free form;and in other cases,a balance is established between the two competing processes at large times and a k(t) obeys a modified scaling law.     

Kinetic Behaviors of a Competitive Population and Fitness System in Exchange-Driven Growth

We proposed an aggregation model of two species aggregates of fitness and population to study the interaction between the two species in their exchange-driven processes of the same species by introducing the monomer birth of fitness catalyzed by the population, where the fitness aggregates perform self-death process and the population aggregates perform self-birth process. The kinetic behaviors of the aggregate size distributions of the fitness and population were analyzed by the rate equation approach with their exchange rate kernel K1（k,j） = K1kj and K2（k,j） = K2kj, the fitness aggregate＇s self-death rate kernel J1 （ k ） = J1 k, population aggregate＇s self-birth rate kernel J2（ k ） = J2k and population-catalyzed fitness birth rate kernel I（k,j） = Ikj＇. The kinetic behavior of the fitness was found depending crucially on the parameter v, which reflects the dependence of the population-catalyzed fitness birth rate on the size of the catalyst （population） aggregate. （i） In the v ≤ 0 case, the effect of catalyzed-birth of fitness is rather weak and the exchange-driven aggregation and self-death of the fitness dominate the process, and the fitness aggregate size distribution αk（t） does not have scale form. （ii） When v ≥0, the effect of the population-catalyzed birth of fitness gets strong enough, and the catalyzed-birth and self-death of the fitness aggregates, together with the self-birth of the population aggregates dominate the evolution process of the fitness aggregates. The aggregate size distribution αk （t） approaches a generalized scaling form.     

A fractal approach to low velocity non-Darcy flow in a low permeability porous medium

In this paper, the mechanism for fluid flow at low velocity in a porous medium is analyzed based on plastic flow of oil in a reservoir and the fractal approach. The analytical expressions for flow rate and velocity of non-Newtonian fluid flow in the low permeability porous medium are derived, and the threshold pressure gradient （TPG） is also obtained. It is notable that the TPG （J） and permeability （K） of the porous medium analytically exhibit the scaling behavior J ~ K-D＇r/（l＋Or）, where DT is the fractal dimension for tortuous capillaries. The fractal characteristics of tortuosity for capillaries should be considered in analysis of non-Darcy flow in a low permeability porous medium. The model predictions of TPG show good agreement with those obtained by the available expression and experimental data. The proposed model may be conducible to a better understanding of the mechanism for nonlinear flow in the low permeability porous medium.     

Invariant Sets and Exact Solutions to Nonlinear Diffusion Equations with x-Dependent Convection and Absorption

In this paper, we introduce a new invariant set Eo={u：ux=f＇（x）F（u）＋ε[g＇（x）-f＇（x）g（x）]F（u）×exp（-∫^u1/F（z）dz）}where f and g are some smooth functions of x, ε is a constant, and F is a smooth function to be determined. The invariant sets and exact sohltions to nonlinear diffusion equation ut = （ D（u）ux）x ＋ Q（x, u）ux ＋ P（x, u）, are discussed. It is shown that there exist several classes of solutions to the equation that belong to the invariant set Eo.     

Coverage Dependence of Growth Mode in Heteroepitaxy of Ni on Cu（100） Surface

A realistic kinetic Monte Carlo （KMC） simulation model with physical parameters is developed, which well reproduces the heteroepitaxial growth of multilayered Ni thin film on Cu（100） surfaces at room temperature. The effects of mass transport between interlayers and edge diffusion of atoms along the islands are included in the simulation model, and the surface roughness and the layer distribution versus total coverage are calculated. Speeially, the simulation model reveals the transition of growth mode with coverage and the difference between the Ni heteroepitaxy on Cu（100） and the Ni homoepitaxy on Ni（100）. Through comparison of KMC simulation with the real scanning tunneling microscopy （STM） experiments, the Ehrlich-Schwoebel （ES） barrier Ees is estimated to be 0.18±0.02 eV for Ni/Cu（100） system while 0.28 eV for Ni/Ni（100）. The simulation also shows that the growth mode depends strongly on the thickness of thin film and the surface temperature, and the critical thickness of growth mode transition is dependent on the growth condition such as surface temperature and deposition flux as well.     

Step instability of the surface during Ino.2Gao.sAs （001） annealing

Anisotropic evolution of the step edges on the compressive-strained In0.2Ga0.8As/GaAs（001） surface has been investigated by scanning tunneling microscopy （STM）. The experiments suggest that step edges are indeed sinuous and protrude somewhere a little way along the [110] direction, which is different from the classical waviness predicted by the theoretical model. We consider that the monatomic step edges undergo a morphological instability induced by the anisotropic diffusion of adatoms on the terrace during annealing, and we improve a kinetic model of step edge based on the classical Burton Cabrer-Frank （BCF） model in order to determine the normal velocity of step enlargement. The results show that the normal velocity is proportional to the arc length of the peninsula, which is consistent with the first result of our kinetic model. Additionally, a significant phenomenon is an excess elongation along the [110] direction at the top of the peninsula with a higher aspect ratio, which is attributed to the restriction of diffusion lengths.     

Critical Behaviors in a Stochastic Local Limited One-Dimensional Rice-Pile Model

A stochastic local limited one-dimensional rice-pile model is numerically investigated. The distributions for avalanche sizes have a clear power-law behavior and it displays a simple finite size scaling. We obtain the avalanche exponents Ts= 1.54±0.10,βs = 2.17±0.10 and TT = 1.80±0.10, βT =1.46 ± 0.10. This self-organized critical model belongs to the same universality class with the Oslo rice-pile model studied by K. Christensen et al. [Phys. Rev. Lett. 77 （1996） 107], a rice-pile model studied by L.A.N. Amaral et al. [Phys. Rev. E 54 （1996） 4512], and a simple deterministic self-organized critical model studied by M.S. Vieira [Phys. Rev. E 61 （2000） 6056].     

Solvable Catalyzed Birth-Death-Exchange Competition Model of Three Speciese

A competition model of three species in exchange-driven aggregation growth is proposed. In the model, three distinct aggregates grow by exchange of monomers and in parallel, birth of species A is catalyzed by species B and death of species A is catalyzed by species C. The rates for both catalysis processes are proportional to kj^v and ky respectively, where ν（ω） is a parameter reflecting the dependence of the catalysis reaction rate of birth （death） on the catalyst aggregate＇s size. The kinetic evolution behaviors of the three species are investigated by the rate equation approach based on the mean-field theory： The form of the aggregate size distribution of A-species αk（t） is found to be dependent crucially on the two catalysis rate kernel parameters. The results show that （i） in case of ν ≤O, the form of ak （t） mainly depends on the competition between self-exchange of species A and species-C-catalyzed death of species A; （ii） in case of ν 〉 0, the form of αk（t） mainly depends on the competition between species-B-catalyzed birth of species A and species-C-catalyzed death of species A.     

Thermodynamic Performance Characteristics of a Brownian Microscopic Heat Engine Driven by Discrete and Periodic Temperature Field

A Brownian microscopic heat engine with a particle hopping on a one-dimensional lattice driven by a discrete and periodic temperature field in a periodic sawtooth potential is investigated. In order to clarify the underlying physical pictures of the heat engine, the heat flow via the potential energy and the kinetic energy of the particles are considered simultaneously. Based on describing the jumps among the three states, the expressions of the efficiency and power output of the heat engine are derived analytically. The general performance characteristic curves are plotted by numerical calculation. It is found that the power output-efficiency curve is a loop-shaped one, which is similar to one for a real irreversible heat engine. The influence of the ratio of the temperature of the hot and cold reservoirs and the sawtooth potential on the maximum efficiency and power output is analyzed for some given parameters. When the heat flows via the kinetic energy is neglected, the power output-efficiency curve is an open-shaped one, which is similar to one for an endroeversible heat engine.     

Dynamics Behaviors and Scaling in Intermittent Turbulence of a Shell Model

In this paper, the dynamics behaviors on fo-δ parameter surface is investigated for Gledzer-Ohkitani- Yamada model We indicate the type of intermittency chaos transitions is saddle node bifurcation. We plot phase diagram on fo-δ parameter surface, which is divided into periodic, quasi-periodic, and intermittent chaos areas. By means of varying Taylor-microscale Reynolds number, we calculate the extended self-similarity of velocity structure function.     

Effective-Field Theory for Kinetic Ising Model on Honeycomb Lattice

As an analytical method, the effective-field theory （EFT） is used to study the dynamical response of the kinetic Ising model in the presence of a sinusoidal oscillating field. The effective-field equations of motion of the average magnetization are given for the honeycomb lattice （Z = 3）. The Liapunov exponent A is calculated for discussing the stability of the magnetization and it is used to determine the phase boundary. In the field amplitude ho / ZJ-temperature T/ZJ plane, the phase boundary separating the dynamic ordered and the disordered phase has been drawn. In contrast to previous analytical results that predicted a tricritical point separating a dynamic phase boundary line of continuous and discontinuous transitions, we find that the transition is always continuous. There is inconsistency between our results and previous analytical results, because they do not introduce sufficiently strong fluctuations.     

Swelling-collapse transition of self-attracting walks

We study the structural properties of self-attracting walks in d dimensions using scaling arguments and Monte Carlo simulations. We find evidence of a transition analogous to the Theta transition of polymers. Above a critical attractive interaction u(c), the walk collapses and the exponents nu and k, characterizing the scaling with time t of the mean square end-to-end distance approximately t(2nu) and the average number of visited sites approximately t(k), are universal and given by nu=1/(d+1) and k=d/(d+1). Below u(c), the walk swells and the exponents are as with no interaction, i.e., nu=1/2 for all d, k=1/2 for d=1 and k=1 for d>/=2. At u(c), the exponents are found to be in a different universality class.     

全局耦合网络上的两种类集团的聚集-湮没反应动力学

利用Monte-Carlo模拟研究了全局耦合网络上扩散限制的不可逆聚集-湮没过程的动力学行为. 在系统中, 同种类集团相遇, 将发生聚集反应; 不同种类的集团相遇, 则发生部分湮没反应. 模拟结果表明:1) 当两种粒子初始浓度相等时, 系统长时间演化后, 集团浓度c(t)和粒子浓度g(t)呈现幂律形式, c(t)～t^{- α}和g(t)～t^-β, 其中幂指数α 和β 满足α=2β 的关系, 且α=2/(2 + q); 集团大小分布随时间的演化满足标度律, a_kt)=k^-τt^-ω\varPhi (k/t^z), 其中τ≈-1.27q, ω≈(3 + 1.27q)/(2 + q), z=α/2=1/(2 + q); 2) 当两种粒子初始浓度不相等时, 系统经长时间演化后, 初始浓度较小的种类完全湮没, 而初始浓度较大的那个种类的集团浓度c_A(t)仍具有幂律形式, c_A(t)～t^-α, 其中α=1/(1+q), 其集团大小分布随时间的演化也满足标度律, 标度指数为τ≈-1.27q, ω≈(2 + 1.27q)/(1 + q)和z=α=1/(1 + q). 模拟结果与已报道的理论分析结果相符得很好.     

The Universality of the Critical Dynamics of the Kinetic Ising Model on the Regular Fractals

The form of the universal scaling law of the critical dynamic exponent, z = D_ƒ + 2/υ, is found on a family of regular fractals by the exact TDRG method. Here, we generate a regular fractal by an anisotropic growing process. Identifying the growing probabilities as the interactions between Ising spins on the fractals, we map the growing probability clouds as a group of the anisotropic Ising Hamiltonians. Applying the RG transformations, we find that the systems of this group of Ising Hamiltonians can be described by two universal static correlation exponents υ₀ = ∞ and υ = 1. So, the growing processes proposed by us capture the essential features in the directed DLA simulations. The studies about their critical dynamic behaviours reveal that unlike the one-dimensional chain the critical dynamics of the kinetic Ising model on the regular fractals is universal. The further discussions show that there is a universal scaling law form of the critical dynamic exponent of the kinetic Ising model, z = D_ƒ + R_max/2υ, on the site models of the regular fractals with R_min = 2. Meanwhile, we discuss Daniel Kandal's correction to the formula of the,critical dynamic exponent in the TDRG method and show that our TDRG calculations are exact.     

Critical exponent crossovers in escape near a bifurcation point

In periodically driven systems, near a bifurcation (critical) point the period-averaged escape rate Wmacr; scales with the field amplitude A as |ln(Wmacr;| proportional, variant (A(c)-A)(xi), where A(c) is a critical amplitude. We find three scaling regions. With increasing field frequency or decreasing |A(c)-A|, the critical exponent xi changes from xi=3/2 for a stationary system to a dynamical value xi=2 and then again to xi=3/2. Monte Carlo simulations agree with the scaling theory.     

Experimental studies of THGEM in different Ar/CO2 mixtures

In this paper, the performance of a type of domestic THGEM （THick Gaseous Electron Multiplier） working in Ar/CO2 mixtures is reported in detail. This kind of single THGEM can provide a gain range from 100 to 1000, which is very suitable for application in neutron detection. In order to study its basic characteristics as a reference for the development of a THGEM based neutron detector, the counting rate plateau, the energy resolution and the gain of the THGEM have been measured in different Ar/CO2 mixtures with a variety of electrical fields. For the Ar/CO2（90%/10%） gas mixture, a wide counting rate plateau is achieved from 720 V to 770 V, with a plateau slope of 2.4%/100 V, and an excellent energy resolution of about 22% is obtained at the 5.9 keV full energy peak of the 55Fe X-ray source.     

Decoherence,Orbit, and Classical Linear Frequency Entropy： Quantum-Classical Correspondence

The decoherence process is analyzed for an open quantum system that is classically chaotic, with a classical linear frequency entropy developed to measure the stability of classical motion. Investigation shows that the decoherence measured by the rate of quantum linear entropy production varies significantly with both the underlying classical orbits and the classical linear frequency entropy. Such correspondence is also supported by the further investigation on the Loschmidt Echo.     

Micro-track structure analysis for 100 MeV Si ions in CR-39 by using atomic force microscopy

To analyze the micro-track structure of heavy ions in a polymer material,parameters including bulk etch rate,track etch rate,etch rate ratio,and track core size were measured.The pieces of CR-39 were exposed to 100 MeV Si ions with normal incidence and were etched in 6.25N NaOH solution at 70 C.Bulk etch rate was read out by a profilemeter after several hours of etching.The other parameters were obtained by using an atomic force microscope(AFM)after a short time of etching.We have measured the second etch pits and minute etch pits to obtain the track growth curve and three dimension track structures to track the core size and etch rate measurements.The local dose of the track core was calculated by theδ-ray theory.In our study,we figure out that the bulk etch rate Vb=(1.58±0.022)μm/h,the track etch rate Vt=(2.90±0.529)μm/h,the etch rate ratio V=1.84±0.031,and the track core radii r≈4.65 nm.In the meantime,we find that the micro-track development violates the traditional track-growth model.For this reason,a scenario is carried out to provide an explanation.