Evolutionary Game Dynamics in a Fitness-Dependent Wright-FisherProcess with Noise

Evolutionary game dynamics in finite size populations can be described by a fitness-dependent Wright-Fisher process. We consider symmetric 2$\times $2 games in a well-mixed population. In our model, two parameters to describe the level of player's rationality and noise intensity in environment are introduced. In contrast with the fixation probability method that used in a noiseless case, the introducing of the noise intensity parameter makes the process an ergodic Markov process and based on the limit distribution of the process, we can analysis the evolutionary stable strategy (ESS) of the games. We illustrate the effects of the two parameters on the ESS of games using the Prisoner's dilemma games (PDG) and the snowdrift games (SG). We also compare the ESS of our model with that of the replicator dynamics in infinite size populations. The results are determined by simulation experiments.     

Some Analytical Properties of the Model for Stochastic Evolutionary Games in Finite Populations with Non-uniform Interaction Rate

Traditional evolutionary games assume uniform interaction rate, which means that the rate at which individuals meet and interact is independent of their strategies. But in some systems, especially biological systems, the players interact with each other discriminately. Taylor and Nowak (2006) were the first to establish the corresponding non-uniform interaction rate model by allowing the interaction rates to depend on strategies. Their model is based on replicator dynamics which assumes an infinite size population. But in reality, the number of individuals in the population is always finite, and there will be some random interference in the individuals' strategy selection process. Therefore, it is more practical to establish the corresponding stochastic evolutionary model in finite populations. In fact, the analysis of evolutionary games in a finite size population is more difficult. Just as Taylor and Nowak said in the outlook section of their paper, "The analysis of non-uniform interaction rates should be extended to stochastic game dynamics of finite populations." In this paper, we are exactly doing this work. We extend Taylor and Nowak's model from infinite to finite case, especially focusing on the influence of non-uniform connection characteristics on the evolutionary stable state of the system. We model the strategy evolutionary process of the population by a continuous ergodic Markov process. Based on the limit distribution of the process, we can give the evolutionary stable state of the system. We make a complete classification of the symmetric 2×2 games. For each case game, the corresponding limit distribution of the Markov-based process is given when noise intensity is small enough. In contrast with most literatures in evolutionary games using the simulation method, all our results obtained are analytical. Especially, in the dominant-case game, coexistence of the two strategies may become evolutionary stable states in our model. This result can be used to explain the emergence of cooperation in the Prisoner is Dilemma Games to some extent. Some specific examples are given to illustrate our results.     

Some Analytical Properties of the Model for Stochastic Evolutionary Games in Finite Populations with Non-uniform Interaction Rate

Traditional evolutionary games assume uniform interaction rate, which means that the rate at which individuals meet and interact is independent of their strategies. But in some systems, especially biological systems, the players interact with each other discriminately. Taylor and Nowak (2006) were the first to establish the corresponding non-uniform interaction rate model by allowing the interaction rates to depend on strategies. Their model is based on replicator dynamics which assumes an infinite size population. But in reality, the number of individuals in the population is always finite, and there will be some random interference in the individuals' strategy selection process. Therefore, it is more practical to establish the corresponding stochastic evolutionary model in finite populations. In fact, the analysis of evolutionary games in a finite size population is more difficult. Just as Taylor and Nowak said in the outlook section of their paper, "The analysis of non-uniform interaction rates should be extended to stochastic game dynamics of finite populations." In this paper, we are exactly doing this work. We extend Taylor and Nowak's model from infinite to finite case, especially focusing on the influence of non-uniform connection characteristics on the evolutionary stable state of the system. We model the strategy evolutionary process of the population by a continuous ergodic Markov process. Based on the limit distribution of the process, we can give the evolutionary stable state of the system. We make a complete classification of the symmetric 2×2 games. For each case game, the corresponding limit distribution of the Markov-based process is given when noise intensity is small enough. In contrast with most literatures in evolutionary games using the simulation method, all our results obtained are analytical. Especially, in the dominant-case game, coexistence of the two strategies may become evolutionary stable states in our model. This result can be used to explain the emergence of cooperation in the Prisoner is Dilemma Games to some extent. Some specific examples are given to illustrate our results.     

Evolutionary games on scale-free networks with a preferential selection mechanism

Considering the heterogeneity of individuals’ influence in the real world, we introduce a preferential selection mechanism to evolutionary games (the Prisoner’s Dilemma Game and the Snowdrift Game) on scale-free networks and focus on the cooperative behavior of the system. In every step, each agent chooses an individual from all its neighbors with a probability proportional to k^α indicating the influence of the neighbor, where k is the degree. Simulation results show that the cooperation level has a non-trivial dependence on α. To understand the effect of preferential selection mechanism on the evolution of the system, we investigate the time series of the cooperator frequency in detail. It is found that the cooperator frequency is greatly influenced by the initial strategy of hub nodes when α>0. This observation is confirmed by investigating the system behavior when some hub nodes’ strategies are fixed.     

Exploring Evolutionary Fitness in Biological Systems Using Machine Learning Methods

Here, we propose a computational approach to explore evolutionary fitness in complex biological systems based on empirical data using artificial neural networks. The essence of our approach is the following. We first introduce a ranking order of inherited elements (behavioral strategies or/and life history traits) in considered self-reproducing systems: we use available empirical information on selective advantages of such elements. Next, we introduce evolutionary fitness, which is formally described as a certain function reflecting the introduced ranking order. Then, we approximate fitness in the space of key parameters using a Taylor expansion. To estimate the coefficients in the Taylor expansion, we utilize artificial neural networks: we construct a surface to separate the domains of superior and interior ranking of pair inherited elements in the space of parameters. Finally, we use the obtained approximation of the fitness surface to find the evolutionarily stable (optimal) strategy which maximizes fitness. As an ecologically important study case, we apply our approach to explore the evolutionarily stable diel vertical migration of zooplankton in marine and freshwater ecosystems. Using machine learning we reconstruct the fitness function of herbivorous zooplankton from empirical data and predict the daily trajectory of a dominant species in the northeastern Black Sea.     

Hydrogenic impurity bound polaron in a parabolic quantum dot with arbitrary electron-phonon coupling strength

A variational approach is employed to obtain the ground and the first excited state binding energies of an electron bound to a hydrogenic impurity in a polar semiconductor quantum dot (QD) with symmetric parabolic confinement in both two and three-dimensions. We perform calculations for the entire range of the electron-phonon coupling constant and the Coulomb binding parameter and for arbitrary confinement length. It is found that the binding energy of ground and first excited state is larger in a two-dimension (2D) dot than in a three-dimension (3D) dot and this trend is more pronounced with the increase of the electron-phonon coupling constant for the same value of the Coulomb binding parameter and confinement length. Furthermore, the ground and the first excited state binding energy increases with increasing the Coulomb binding parameter in both 2D and 3D QDs for the same electron-phonon coupling constant.     

理想气体与范氏气体在任意准静态过程中的摩尔热容

导出了理想气体与范氏气体在任意准静态过程中的摩尔热容公式。     

An extension of the direct method and similarity reductions of a generalized Burgers equation with an arbitrary derivative function

In this paper,we extend the well-known direct method proposed by Clarkson and Kruskal for finding similarity reductions of partial differential equations.It follows that some new similarity reductions of the generalized Burgers equation,such as travelling wave reduction,logarithmic reduction,power reduction,rational fractional reduction,etc,are derived,in which some of these cannot be obtained solely by using the direct method.The similarity reductions obtained are interpreted by the nonclassical symmetry Lie group.     

Optical state selection process with optical pumping in a cesium atomic fountain clock

We propose and realize a new optical state selection method on a cesium atomic fountain clock by applying a two-laser 3-3' optical pumping configuration to spin polarize atoms. The atoms are prepared in |F=3, m_F=0> clock state with optical pumping directly after being launched up, followed by a pushing beam to push away the atoms remaining in the |F=4> state. With a state selection efficiency exceeding 92%, this optical method can substitute the traditional microwave state selection, and helps to develop a more compact physical package. A Ramsey fringe has been achieved with this optical state selection method, and a contrast of 90% is obtained with a full width half maximum of 0.92 Hz. The short-term frequency stability of 6.8×10^-14 (τ/s)^-1/2 is acquired. In addition, the number of detected atoms is increased by a factor of 1.7 with the optical state selection.     

Bound polaron in a two- and three-dimensional quantum dot for arbitrary electron–phonon coupling strength

A variational approach is presented for calculating the ground-state (GS) binding energies of an electron bound to a Coulomb impurity in a polar semiconductor quantum dot (QD) with parabolic confinement in both two and three dimensions. We perform calculations for the entire range of the electron-phonon coupling constant and the Coulomb binding parameter and for arbitrary confinement length. It is found that the polaronic effect is stronger in a two dimensions (2D) dot than in a three dimensions (3D) dot and this trend is more pronounced with the increase of the coupling constant. Furthermore, the GS binding energy increases with increasing the Coulomb binding parameter in both 2D and 3D QDs for the same electron–phonon coupling constant. The results also indicate that this effect becomes much more pronounced with decreasing dimensionality.     

Zone inhomogeneity with the random asymmetric simple exclusion process in a one-lane system

In this paper we use theoretical analysis and extensive simulations to study zone inhomogeneity with the random asymmetric simple exclusion process (ASEP). In the inhomogeneous zone, the hopping probability is less than 1. Two typical lattice geometries are investigated here. In case A, the lattice includes two equal segments. The hopping probability in the left segment is equal to 1, and in the right segment it is equal to p, which is less than 1. In case B, there are three equal segments in the system; the hopping probabilities in the left and right segments are equal to 1, and in the middle segment it is equal to p, which is less than 1. Through theoretical analysis, we can discover the effect on these systems when p is changed.     

Volatility estimators and the inverse range process in a random volatility random walk and Wiener processes

The purpose of this paper is to study the mean, the variance, the probability distribution and the hazard rate of the inverse range process of an a-priori unknown volatility random walk. Motivation for this process arises when it is necessary to obtain statistics that pertain to a process volatility in addition to the usual variance statistics. As a result, range process statistics are indicated as an additional source of information in the study of processes’ volatility. Examples and applications are considered.     

Combinatorial and approximative analyses in a spatially random division process

For a spatial characteristic, there exist commonly fat-tail frequency distributions of fragment-size and -mass of glass, areas enclosed by city roads, and pore size/volume in random packings. In order to give a new analytical approach for the distributions, we consider a simple model which constructs a fractal-like hierarchical network based on random divisions of rectangles. The stochastic process makes a Markov chain and corresponds to directional random walks with splitting into four particles. We derive a combinatorial analytical form and its continuous approximation for the distribution of rectangle areas, and numerically show a good fitting with the actual distribution in the averaging behavior of the divisions.     

有源等离子体在开放大气环境下的反应扩散过程研究

利用一维模型用数值模拟的方法,研究了低气压开放环境下大气等离子体在存在等离子体源的情况下的反应扩散过程.得到了考虑化学反应、扩散以及漂移共同作用下的大气等离子体的主要成分随等离子体注入流量的变化规律及一定流量情况下主要带电成分随时空的演化规律. 将数值模拟结果与一个近似解析公式相衔接,估计达到稳态时维持一定电子密度所需要的等离子体流量,可据此进一步估计所需功率. 关键词： 大气等离子体 等离子体源 数值模拟 反应扩散过程     

二能级原子与多模光场简并多光子共振相互作用系统中量子保真度的演化特性

利用全量子理论和数值计算方法研究了多模相干态光场与单个二能级原子通过任意N_j度简并N~∑光子共振相互作用系统中量子保真度的时间演化特性,给出了三模场与原子相互作用过程中光场和原子保真度的数值计算结果,详细讨论了初始平均光子数、原子分布角、原子偶极相位角、光场激发角以及原子简并度等对量子保真度的影响.数值计算结果表明:以上诸多因素对量子保真度影响的结果均导致其发生振荡性变化.光场和原子保真度随着初始光场增强而急剧减小,说明初始光强敏感地影响着保真度的大小;量子保真度的变化快慢程度强烈地依赖于原子简并度及场一原子的耦合系数;原子分布角、光场激发角不同程度地对量子保真度的大小和频率有所影响;而原子偶极相位角的变化对场和原子的量子保真度几乎没有影响.根据这些特性,通过某些条件的约束可以适当控制保真度变化的快慢及其大小.     

Reactive process and hysteresis effect in magnetron with magnetized hollow cathode enhanced target

The hysteresis effect in the reactive process was investigated in the magnetron with a magnetized hollow cathode enhanced target (HoCET) in which the target is coupled with the hollow cathode magnetized by the magnetic field of the magnetron. The process, where both the magnetron and hollow cathode plasmas are combined, is compared to the magnetron sputtering. The hysteresis curve in the magnetized HoCET magnetron, recording the titanium emission intensity versus varying content of nitrogen in the gas mixture exhibits a local maximum on the increasing part of the curve. The hysteresis curve is shifted to lower contents of nitrogen than the hysteresis curve for the magnetron. It is concluded that more efficient utilization of the reactive gas takes place in this device.     

Pattern selection in a predation model with self and cross diffusion

In this paper, we present the amplitude equations for the excited modes in a cross-diffusive predator--prey model with zero-flux boundary conditions. From these equations, the stability of patterns towards uniform and inhomogenous perturbations is determined. Furthermore, we present novel numerical evidence of six typical turing patterns, and find that the model dynamics exhibits complex pattern replications: for μ₁<μ≤μ₂, the steady state is the only stable solution of the model; for μ₂<μ≤μ₄, by increasing the control parameter μ, the sequence H_π-hexagons → H₀-hexagon-stripe mixtures rightarrow stripes → H_π-hexagon-stripe mixtures → H₀-hexagons is observed; for μ>μ₄, the stripe pattern emerges. This may enrich the pattern formation in the cross-diffusive predator--prey model.     

Analytical and numerical investigations of displaced thermal state evolutions in a laser process

We investigate how displaced thermal states(DTSs) evolve in a laser channel. Remarkably, the initial DTS, an example of a mixed state, still remains mixed and thermal. At long times, they finally decay to a highly classical thermal field only related to the laser parameters κ and g. The normal ordering product of density operator of the DTS in the laser channel leads to obtaining the analytical time-evolution expressions of the photon number, Wigner function, and von Neumann entropy. Also, some interesting results are presented via numerically investigating these explicit time-dependent expressions.     

充模过程中熔接痕的改进光滑粒子动力学方法模拟与预测

提出了一种黏弹性流体的改进光滑粒子动力学(SPH)方法以试探性地模拟和预测黏弹性FENE-P熔体充模过程中熔接痕的形态演化. 首先基于SPH方法建立了聚合物流动的宏微观耦合模型, 同时提出了黏弹性流体的改进SPH离散格式. 随后, 通过模拟一些基准算例验证了改进的SPH方法模拟聚合物宏微观耦合问题的有效性及收敛性, 以及所提出的黏弹性温度模型的有效性. 最后, 模拟了环型腔内的充模过程, 试探性地展示了充模过程中微观分子的变形过程. 同时采用顺序热流道技术模拟了多浇口C形腔内的充模过程, 并与其他数值结果做比较. 数值结果表明: 对于大制件多浇口充模过程, 顺序热流道技术能够改善甚至消除充模过程中的熔接痕.