Analysis of epidemic models with demographics in metapopulation networks

In this paper, two susceptible–infected–susceptible (SIS) epidemic models are presented and analyzed by reaction–diffusion processes with demographics in metapopulation networks. Firstly, an SIS model with constant-inputting is discussed. The model has a disease-free equilibrium, which is locally asymptotically stable when the basic reproduction number is less than unity, otherwise it is unstable. It has an endemic equilibrium, which is globally asymptotically stable. Secondly, in another SIS model, the birth rate is the form of Logistic. Similarly, the stability of disease-free equilibrium and endemic equilibrium is also proved. Finally, numerical simulations are performed to illustrate the analytical results.     

自适应网络中病毒传播的稳定性和分岔行为研究

自适应复杂网络是以节点状态与拓扑结构之间存在反馈回路为特征的网络. 针对自适应网络病毒传播模型, 利用非线性微分动力学系统研究病毒传播行为; 通过分析非线性系统对应雅可比矩阵的特征方程, 研究其平衡点的局部稳定性和分岔行为, 并推导出各种分岔点的计算公式. 研究表明, 当病毒传播阈值小于病毒存在阈值, 即R₀0^c时, 网络中病毒逐渐消除, 系统的无病毒平衡点是局部渐近稳定的; R₀^c0<1时, 网络出现滞后分岔, 产生双稳态现象, 系统存在稳定的无病毒平衡点、较大稳定的地方病平衡点和较小不稳定的地方病平衡点; R₀>1时, 网络中病毒持续存在, 系统唯一的地方病平衡点是局部渐近稳定的. 研究发现, 系统先后出现了鞍结分岔、跨临界分岔、霍普夫分岔等分岔行为. 最后通过数值仿真验证所得结论的正确性. 关键词： 自适应网络 稳定性 分岔 基本再生数     

移动环境下网络病毒传播模型及其稳定性研究

考虑网络节点的随机移动, 基于平均场理论 提出一个移动环境下网络病毒传播的数学模型, 利用微分动力学系统理论研究了病毒传播行为. 研究表明, 当病毒基本再生数R₀ ≤ 1时, 网络中病毒逐渐消除, 系统的无病毒平衡点全局渐进稳定; 当R₀ > 1时, 网络中病毒持续存在, 系统的地方病平衡点全局渐进稳定.通过仿真验证了所得结论的正确性.     

Transition state theory: A generalization to nonequilibrium systems with power-law distributions

Transition state theory (TST) is generalized to nonequilibrium systems with power-law distributions. The stochastic dynamics that gives rise to the power-law distributions for the reaction coordinate and momentum is modeled by Langevin equations and corresponding Fokker-Planck equations. It is considered that a system far away from equilibrium does not have to relax to a thermal equilibrium state with Boltzmann-Gibbs distribution, but asymptotically approaches a nonequilibrium stationary state with a power-law distribution. Thus, we obtain a possible generalization of TST rates to nonequilibrium systems with power-law distributions. Furthermore, we derive the generalized TST rate constants for one-dimensional and n-dimensional Hamiltonian systems away from equilibrium, and obtain a generalized Arrhenius rate for systems with power-law distributions.     

Impact of media coverage on epidemic spreading in complex networks

An SIS network model incorporating the influence of media coverage on transmission rate is formulated and analyzed. We calculate the basic reproduction number _R0$R_0$ by utilizing the local stability of the disease-free equilibrium. Our results show that the disease-free equilibrium is globally asymptotically stable and that the disease dies out if _R0$R_0$ is below 1; otherwise, the disease will persist and converge to a unique positive stationary state. This result may suggest effective control strategies to prevent disease through media coverage and education activities in finite-size scale-free networks. Numerical simulations are also performed to illustrate our results and to give more insights into the dynamical process.     

Global dynamics of a novel multi-group model for computer worms

In this paper, we study worm dynamics in computer networks composed of many autonomous systems. A novel multigroup SIQR (susceptible-infected-quarantined-removed) model is proposed for computer worms by explicitly considering anti-virus measures and the network infrastructure. Then, the basic reproduction number of worm R0 is derived and the global dynamics of the model are established. It is shown that if R0 is less than or equal to 1, the disease-free equilibrium is globally asymptotically stable and the worm dies out eventually, whereas, if R0 is greater than 1, one unique endemic equilibrium exists and it is globally asymptotically stable, thus the worm persists in the network. Finally, numerical simulations are given to illustrate the theoretical results.     

On the global asymptotic stability of a system of N generalized chemical rate equations

It is proved that a system of N generalized chemical rate equations, which for N = 2 describes the kinetics of irradiation-produced point defects, has a unique equilibrium point in the part of phase space where all dependent variables are positive. By finding an appropriate Lyapunov function, it is also shown that this equilibrium point is globally asymptotically stable, for all positive initial conditions of the system.     

Critical Behavior of the Kramers Escape Rate in Asymmetric Classical Field Theories

We introduce an asymmetric classical Ginzburg–Landau model in a bounded interval, and study its dynamical behavior when perturbed by weak spatiotemporal noise. The Kramers escape rate from a locally stable state is computed as a function of the interval length. An asymptotically sharp second-order phase transition in activation behavior, with corresponding critical behavior of the rate prefactor, occurs at a critical length ? _c , similar to what is observed in symmetric models. The weak-noise exit time asymptotics, to both leading and subdominant orders, are analyzed at all interval lengthscales. The divergence of the prefactor as the critical length is approached is discussed in terms of a crossover from non-Arrhenius to Arrhenius behavior as noise intensity decreases. More general models without symmetry are observed to display similar behavior, suggesting that the presence of a “phase transition” in escape behavior is a robust and widespread phenomenon.     

New unified formulation of transient phenomena near the instability point on the basis of the fokker-planck equation

A new unified theory of transient phenomena is proposed to treat the passage from an initially unstable state to a final stable state. In the nonlinear intermediate time region, this is reduced to the scaling theory by the present author, and for t→∞ it gives a correct fluctuation asymptotically.     

Non-equilibrium behavior at a liquid-gas critical point

Second-order phase transitions in a non-equilibrium liquid-gas model with reversible mode couplings, i.e., model H for binary-fluid critical dynamics, are studied using dynamic field theory and the renormalization group. The system is driven out of equilibrium either by considering different values for the noise strengths in the Langevin equations describing the evolution of the dynamic variables (effectively placing these at different temperatures), or more generally by allowing for anisotropic noise strengths, i.e., by constraining the dynamics to be at different temperatures in d _|| - and d _⊥-dimensional subspaces, respectively. In the first, isotropic case, we find one infrared-stable and one unstable renormalization group fixed point. At the stable fixed point, detailed balance is dynamically restored, with the two noise strengths becoming asymptotically equal. The ensuing critical behavior is that of the standard equilibrium model H. At the novel unstable fixed point, the temperature ratio for the dynamic variables is renormalized to infinity, resulting in an effective decoupling between the two modes. We compute the critical exponents at this new fixed point to one-loop order. For model H with spatially anisotropic noise, we observe a critical softening only in the d _⊥-dimensional sector in wave vector space with lower noise temperature. The ensuing effective two-temperature model H does not have any stable fixed point in any physical dimension, at least to one-loop order. We obtain formal expressions for the novel critical exponents in a double expansion about the upper critical dimension d _c = 4 - d _|| and with respect to d _|| , i.e., about the equilibrium theory. Received 4 April 2002 Published online 13 August 2002     

The Vlasov equation and the Hamiltonian mean-field model

We show that the quasi-stationary states of homogeneous (zero magnetization) states observed in the N-particle dynamics of the Hamiltonian mean-field (HMF) model are nothing but Vlasov stable homogeneous states. There is an infinity of Vlasov stable homogeneous states corresponding to different initial momentum distributions. Tsallis q-exponentials in momentum, homogeneous in angle, distribution functions are possible, however, they are not special in any respect, among an infinity of others. All Vlasov stable homogeneous states lose their stability because of finite N effects and, after a relaxation time diverging with a power-law of the number of particles, the system converges to the Boltzmann–Gibbs equilibrium.     

Symmetry breaking bifurcations in a circular chain of N coupled logistic maps

A circular chain of N cells with logistic dynamics, coupled together with symmetric nearest neighbor coupling and periodic boundary conditions is investigated. For certain coupling parameters we observe bifurcation of a stable state into two types of period two solutions. By using the symmetry of this Coupled Map Lattice model, we show that the bifurcated system only can have periodic solutions with symmetry group corresponding to certain subgroups of the full symmetry group of the system.     

Positive Commutators in Non-Equilibrium Quantum Statistical Mechanics

The method of positive commutators, developed for zero temperature problems over the last twenty years, has been an essential tool in the spectral analysis of Hamiltonians in quantum mechanics. We extend this method to positive temperatures, i.e. to non-equilibrium quantum statistical mechanics. We use the positive commutator technique to give an alternative proof of a fundamental property of a certain class of large quantum systems, called Return to Equilibrium. This property says that equilibrium states are (asymptotically) stable: if a system is slightly perturbed from its equilibrium state, then it converges back to that equilibrium state as time goes to infinity. Received: 27 December 2000 / Accepted: 21 June 2001     

Ab initio theory of complex electronic ground state of superconductors: II.: Antiadiabatic state—Ground state of superconductors

Based on Q, P-dependent modification of the Born-Oppenheimer approximation (BOA), the ab initio theory of complex electronic ground state of superconductors is presented. As an illustrative example, application of the theory to superconductors of a different character and to the corresponding nonsuperconducting analogues is presented. It is shown that due to electron-phonon (EP) interactions, which drive system from adiabatic into antiadiabatic state, adiabatic translation symmetry is broken and system is stabilized in antiadiabatic state at distorted geometry with respect to adiabatic equilibrium high-symmetry structure. Stabilization effect in the antiadiabatic state is due to strong dependence of the electronic motion on the instantaneous nuclear kinetic energy, i.e. on the effect that is neglected on the adiabatic level within the BOA. At distorted geometry, antiadiabatic ground state is geometrically degenerated with fluxional nuclear configurations in the phonon modes that drive system into this state. It has been shown that until the system remains in antiadiabatic state, nonadiabatic polaron-renormalized phonon interactions are zero in the well-defined k-region of reciprocal lattice. This, along with geometric degeneracy of the antiadiabatic state, enables formation of mobile bipolarons that can move over lattice as supercarriers without dissipation. Moreover, it has been shown that due to EP interactions at transition into antiadiabatic state, k-dependent gap in one-electron spectrum has been opened. Gap opening is associated with shift of the original adiabatic Hartree-Fock orbital energies and with the k-dependent change in density of states of particular band(s) at Fermi level. Corrected one-particle spectrum enables to derive thermodynamic properties in full agreement with corresponding thermodynamic properties of superconductors.Based on the complex ab initio theory, it has been shown that Fröhlich's effective attractive electron-electron interaction term represents correction to electron correlation energy at transition from adiabatic into antiadiabatic state due to EP interactions. It has been shown that increased electron correlation is a consequence of stabilization of the system in superconducting electronic ground state, but not the reason for its formation.     

Relaxation of vacuum energy in <Emphasis Type="Italic">q</Emphasis>-theory

The q-theory formalism aims to describe the thermodynamics and dynamics of the deep quantum vacuum. The thermodynamics leads to an exact cancellation of the quantum-field zero-point-energies in equilibrium, which partly solves the main cosmological constant problem. But, with reversible dynamics, the spatially flat Friedmann–Robertson–Walker universe asymptotically approaches the Minkowski vacuum only if the Big Bang already started out in an initial equilibrium state. Here, we extend q-theory by introducing dissipation from irreversible processes. Neglecting the possible instability of a de-Sitter vacuum, we obtain different scenarios with either a de-Sitter asymptote or collapse to a final singularity. The Minkowski asymptote still requires fine-tuning of the initial conditions. This suggests that, within the q-theory approach, the decay of the de-Sitter vacuum is a necessary condition for the dynamical solution of the cosmological constant problem.     

Trajectory of neutron–neutron–C excited three-body state

The trajectory of the first excited Efimov state is investigated by using a renormalized zero-range three-body model for a system with two bound and one virtual two-body subsystems. The approach is applied to n–n–¹⁸C, where the n–n$n\text{–}n$ virtual energy and the three-body ground state are kept fixed. It is shown that such three-body excited state goes from a bound to a virtual state when the n–¹⁸C binding energy is increased. Results obtained for the n–¹⁹C elastic cross-section at low energies also show dominance of an S-matrix pole corresponding to a bound or virtual Efimov state. It is also presented a brief discussion of these findings in the context of ultracold atom physics with tunable scattering lengths.     

Stability and quasi-equidistant propagation of NLS soliton trains

Using the complex Toda chain we model the asymptotic behavior of the N soliton pulse trains of the nonlinear Schrödinger equation. Stable asymptotic regimes are: (i) asymptotically free propagation of all N solitons; (ii) bound state regime where the N solitons may move quasi-equidistantly (QED); and (iii) various intermediate regimes. Our method allows one to determine analytically the set of initial soliton parameters corresponding to each regime. We list the soliton parameters, which ensure QED propagation of all N solitons since this is important for optical fiber communication.     

Ground state energy of the electron-hole liquid in type-II (GaAs)m/(AlAs)m quantum wells

The ground state energy of quasi-two-dimensional electron-hole liquid (EHL) at zero temperature is calculated for type-II (GaAs)_m/(AlAs)_m (5≤m≤10) quantum wells (QWs). The correlation effects of Coulomb interaction are taken into account by a random phase approximation of Hubbard. Our EHL ground state energy per electron-hole pair is lower than the exciton energy calculated recently for superlattices, so we expected that EHL is more stable state than excitons at high excitation density. It is also demonstrated that the equilibrium density of EHL in type-II GaAs/AlAs QWs is of one order of magnitude larger than that in type-I GaAs/AlAs QWs.     

Delayed feedback control of time-delayed chaotic systems: Analytical approach at Hopf bifurcation

This Letter is concerned with bifurcation and chaos control in scalar delayed differential equations with delay parameter τ. By linear stability analysis, the conditions under which a sequence of Hopf bifurcation occurs at the equilibrium points are obtained. The delayed feedback controller is used to stabilize unstable periodic orbits. To find the controller delay, it is chosen such that the Hopf bifurcation remains unchanged. Also, the controller feedback gain is determined such that the corresponding unstable periodic orbit becomes stable. Numerical simulations are used to verify the analytical results.     

Fluctuation-driven directed transport in the presence of Lévy flights

The role of Lévy flights on fluctuation-driven transport in time independent periodic potentials with broken spatial symmetry is studied. Two complementary approaches are followed. The first one is based on a generalized Langevin model describing overdamped dynamics in a ratchet type external potential driven by Lévy white noise with stability index α in the range 1<α<2. The second approach is based on the space fractional Fokker-Planck equation describing the corresponding probability density function (PDF) of particle displacements. It is observed that, even in the absence of an external tilting force or a bias in the noise, the Lévy flights drive the system out of the thermodynamic equilibrium and generate an up-hill current (i.e., a current in the direction of the steeper side of the asymmetric potential). For small values of the noise intensity there is an optimal value of α yielding the maximum current. The direction and magnitude of the current can be manipulated by changing the Lévy noise asymmetry and the potential asymmetry. For a sharply localized initial condition, the PDF of staying at the minimum of the potential exhibits scaling behavior in time with an exponent bigger than the −1/α exponent corresponding to the force free case.