Evolutionary Algorithms and Dynamic Programming
Doerr, Benjamin; Eremeev, Anton; Neumann, Frank; Theile, Madeleine; Thyssen, Christian
2013-01-01
Recently, it has been proven that evolutionary algorithms produce good results for a wide range of combinatorial optimization problems. Some of the considered problems are tackled by evolutionary algorithms that use a representation which enables them to construct solutions in a dynamic programming fashion. We take a general approach and relate the construction of such algorithms to the development of algorithms using dynamic programming techniques. Thereby, we give general guidelines on how ...
Neuronal boost to evolutionary dynamics
de Vladar, Harold P.; Szathmáry, Eörs
2015-01-01
Standard evolutionary dynamics is limited by the constraints of the genetic system. A central message of evolutionary neurodynamics is that evolutionary dynamics in the brain can happen in a neuronal niche in real time, despite the fact that neurons do not reproduce. We show that Hebbian learning and structural synaptic plasticity broaden the capacity for informational replication and guided variability provided a neuronally plausible mechanism of replication is in place. The synergy between learning and selection is more efficient than the equivalent search by mutation selection. We also consider asymmetric landscapes and show that the learning weights become correlated with the fitness gradient. That is, the neuronal complexes learn the local properties of the fitness landscape, resulting in the generation of variability directed towards the direction of fitness increase, as if mutations in a genetic pool were drawn such that they would increase reproductive success. Evolution might thus be more efficient within evolved brains than among organisms out in the wild. PMID:26640653
Modeling tumor evolutionary dynamics
Beatriz eStransky
2013-02-01
Full Text Available Tumorigenesis can be seen as an evolutionary process, in which the transformation of a normal cell into a tumor cell involves a number of limiting genetic and epigenetic events, occurring in a series of discrete stages. However, not all mutations in a cell are directly involved in cancer development and it is likely that most of them (passenger mutations do not contribute in any way to tumorigenesis. Moreover, the process of tumor evolution is punctuated by selection of advantageous (driver mutations and clonal expansions. Regarding these driver mutations, it is uncertain how many limiting events are required and / or sufficient to promote a tumorigenic process or what are the values associated with the adaptive advantage of different driver mutations. In spite of the availability of high-quality cancer data, several assumptions about the mechanistic process of cancer initiation and development remain largely untested, both mathematically and statistically. Here we review the development of mathematical/computational models where some assumptions were tested and discuss the impact of these models to the field of tumor biology.
Evolutionary dynamics of mammalian karyotypes
Carlo Alberto Redi
2012-12-01
Full Text Available This special volume of Cytogenetic and Genome Research (edited by Roscoe Stanyon, University of Florence and Alexander Graphodatsky, Siberian division of the Russian Academy of Sciences is dedicated to the fascinating long search of the forces behind the evolutionary dynamics of mammalian karyotypes, revealed after the hypotonic miracle of the 1950s....
Organisations’ evolutionary dynamics: a group dynamics approach
Germán Eduardo Vargas
2010-04-01
Full Text Available Colombian entrepreneurs’ straggling, reactionary and inertial orientation has been inconsistently lustified by the availability of internal and leveraged resources, a concept intensifying deficient technological capacity. Company activity (seen as being a socioeconomic unit has been integrally orientated within an evolutionary framework by company identity and cohesion as well as adaptation and evolutionary mechanisms. The present document uses a group dynamics’ model to illustrate how knowledge-based strategic orientation and integration for innovation have become an imperative for development, from slight leverage, distinguishing between two evolutionary company forms: traditional economic (inertial, as they introduce sporadic incremental improvements and modern companies (dynamic and radical innovators. Revealing conclusions obtained from such model may be used for intervening in and modernising company activity.
Dynamics of a Simple Evolutionary Process
Stauffer, Dietrich; Newman, M. E. J.
We study the simple evolutionary process in which we repeatedly find the least fit agent in a population of agents and give it a new fitness, which is chosen independently at random from a specified distribution. We show that many of the average properties of this process can be calculated exactly using analytic methods. In particular, we find the distribution of fitnesses at arbitrary time, and the distribution of the lengths of runs of hits on the same agent, the latter being found to follow a power law with exponent -1, similar to the distribution of times between evolutionary events in the Bak-Sneppen model and models based on the so-called record dynamics. We confirm our analytic results with extensive numerical simulations.
Evolutionary dynamics of group fairness.
Santos, Fernando P; Santos, Francisco C; Paiva, Ana; Pacheco, Jorge M
2015-08-01
The emergence and impact of fairness is commonly studied in the context of 2-person games, notably the Ultimatum Game. Often, however, humans face problems of collective action involving more than two individuals where fairness is known to play a very important role, and whose dynamics cannot be inferred from what is known from 2-person games. Here, we propose a generalization of the Ultimatum Game for an arbitrary number of players--the Multiplayer Ultimatum Game. Proposals are made to a group of responders who must individually reject or accept the proposal. If the total number of individual acceptances stands below a given threshold, the offer will be rejected; otherwise, the offer will be accepted, and equally shared by all responders. We investigate the evolution of fairness in populations of individuals by means of evolutionary game theory, providing both analytical insights and results from numerical simulations. We show how imposing stringent consensuses significantly increases the value of the proposals, leading to fairer outcomes and more tolerant players. Furthermore, we show how stochastic effects--such as imitation errors and/or errors when assessing the fitness of others--may further enhance the overall success in reaching fair collective action. PMID:25936348
Stochastic evolutionary dynamics of direct reciprocity
Lorens A. Imhof; Nowak, Martin A.
2009-01-01
Evolutionary game theory is the study of frequency-dependent selection. The success of an individual depends on the frequencies of strategies that are used in the population. We propose a new model for studying evolutionary dynamics in games with a continuous strategy space. The population size is finite. All members of the population use the same strategy. A mutant strategy is chosen from some distribution over the strategy space. The fixation probability of the mutant strategy in the reside...
On evolutionary ray-projection dynamics
Joosten, Reinoud; Roorda, Berend
2011-01-01
We introduce the ray-projection dynamics in evolutionary game theory by employing a ray projection of the relative fitness (vector) function, i.e., a projection unto the unit simplex along a ray through the origin. Ray-projection dynamics are weakly compatible in the terminology of Friedman (Econome
Evolutionary Dynamics of Nationalism and Migration
André Barreira da Silva Rocha
2012-01-01
I present a dynamic evolutionary game model to address the relation between nationalism against immigrants and assimilation of the latter into the host country culture. I assume a country composed of two different large polymorphic populations, one of native citizens and the other of immigrants. A native citizen may behave nationalistically or may welcome immigrants. Immigrants may have an interest in learning the host country language or not. Evolution is modelled using replicator dynamics (...
On the Dynamic Foundation of Evolutionary Stability in Continuous Models
Oechssler, Jörg; Riedel, Frank
2000-01-01
We show in this paper that none of the existing static evolutionary stability concepts (ESS, CSS, uninvadability, NIS) is sufficient to guarantee dynamic stability in the weak topology with respect to standard evolutionary dynamics if the strategy space is continuous. We propose a new concept, evolutionary robustness, which is stronger than the previous concepts. Evolutionary robustness ensures dynamic stability for replicator dynamics in doubly symmetric games.
Multiscale structure in eco-evolutionary dynamics
Stacey, Blake C.
In a complex system, the individual components are neither so tightly coupled or correlated that they can all be treated as a single unit, nor so uncorrelated that they can be approximated as independent entities. Instead, patterns of interdependency lead to structure at multiple scales of organization. Evolution excels at producing such complex structures. In turn, the existence of these complex interrelationships within a biological system affects the evolutionary dynamics of that system. I present a mathematical formalism for multiscale structure, grounded in information theory, which makes these intuitions quantitative, and I show how dynamics defined in terms of population genetics or evolutionary game theory can lead to multiscale organization. For complex systems, "more is different," and I address this from several perspectives. Spatial host--consumer models demonstrate the importance of the structures which can arise due to dynamical pattern formation. Evolutionary game theory reveals the novel effects which can result from multiplayer games, nonlinear payoffs and ecological stochasticity. Replicator dynamics in an environment with mesoscale structure relates to generalized conditionalization rules in probability theory. The idea of natural selection "acting at multiple levels" has been mathematized in a variety of ways, not all of which are equivalent. We will face down the confusion, using the experience developed over the course of this thesis to clarify the situation.
A Simple General Model of Evolutionary Dynamics
Thurner, Stefan
Evolution is a process in which some variations that emerge within a population (of, e.g., biological species or industrial goods) get selected, survive, and proliferate, whereas others vanish. Survival probability, proliferation, or production rates are associated with the "fitness" of a particular variation. We argue that the notion of fitness is an a posteriori concept in the sense that one can assign higher fitness to species or goods that survive but one can generally not derive or predict fitness per se. Whereas proliferation rates can be measured, fitness landscapes, that is, the inter-dependence of proliferation rates, cannot. For this reason we think that in a physical theory of evolution such notions should be avoided. Here we review a recent quantitative formulation of evolutionary dynamics that provides a framework for the co-evolution of species and their fitness landscapes (Thurner et al., 2010, Physica A 389, 747; Thurner et al., 2010, New J. Phys. 12, 075029; Klimek et al., 2009, Phys. Rev. E 82, 011901 (2010). The corresponding model leads to a generic evolutionary dynamics characterized by phases of relative stability in terms of diversity, followed by phases of massive restructuring. These dynamical modes can be interpreted as punctuated equilibria in biology, or Schumpeterian business cycles (Schumpeter, 1939, Business Cycles, McGraw-Hill, London) in economics. We show that phase transitions that separate phases of high and low diversity can be approximated surprisingly well by mean-field methods. We demonstrate that the mathematical framework is suited to understand systemic properties of evolutionary systems, such as their proneness to collapse, or their potential for diversification. The framework suggests that evolutionary processes are naturally linked to self-organized criticality and to properties of production matrices, such as their eigenvalue spectra. Even though the model is phrased in general terms it is also practical in the sense
Molecular Dynamics Calculations
1996-01-01
The development of thermodynamics and statistical mechanics is very important in the history of physics, and it underlines the difficulty in dealing with systems involving many bodies, even if those bodies are identical. Macroscopic systems of atoms typically contain so many particles that it would be virtually impossible to follow the behavior of all of the particles involved. Therefore, the behavior of a complete system can only be described or predicted in statistical ways. Under a grant to the NASA Lewis Research Center, scientists at the Case Western Reserve University have been examining the use of modern computing techniques that may be able to investigate and find the behavior of complete systems that have a large number of particles by tracking each particle individually. This is the study of molecular dynamics. In contrast to Monte Carlo techniques, which incorporate uncertainty from the outset, molecular dynamics calculations are fully deterministic. Although it is still impossible to track, even on high-speed computers, each particle in a system of a trillion trillion particles, it has been found that such systems can be well simulated by calculating the trajectories of a few thousand particles. Modern computers and efficient computing strategies have been used to calculate the behavior of a few physical systems and are now being employed to study important problems such as supersonic flows in the laboratory and in space. In particular, an animated video (available in mpeg format--4.4 MB) was produced by Dr. M.J. Woo, now a National Research Council fellow at Lewis, and the G-VIS laboratory at Lewis. This video shows the behavior of supersonic shocks produced by pistons in enclosed cylinders by following exactly the behavior of thousands of particles. The major assumptions made were that the particles involved were hard spheres and that all collisions with the walls and with other particles were fully elastic. The animated video was voted one of two
Calculation of Triggering Angle of Thyristor Rectifiers with Evolutionary Algorithms
Fahri VATANSEVER
2015-05-01
Full Text Available Rectifier circuits have important role in electrical energy systems. Especially in thyristor rectifiers which have capability of generating multiple voltage level, determining/calculating the optimal trigger angle and applying trigger signal at this angles are among the main process. In this study, desired level of average output voltage according to trigger angles is obtained both classically (mathematically solving of equations and using evolutionary algorithms which are genetic algorithms and differential evolutions. In this way, a software can be used in educational purposes which can calculate optimal trigger angles using both mathematically and heuristically, show results and many properties/parameters of circuit graphically and numerically is developed. Analysis/simulations performed with the designed software indicates that evolutionary algorithms can be used in this field effectively and efficiently.
Eco-evolutionary dynamics in a changing world.
Hanski, Ilkka
2012-02-01
Fast evolutionary changes are common in natural populations, though episodes of rapid evolution do not generally last for long and are typically associated with changing environments. During such periods, evolutionary dynamics may influence ecological population dynamics and vice versa. This review is concerned with spatial eco-evolutionary dynamics with a focus on the occurrence of species in marginal habitats and on metapopulations inhabiting heterogeneous environments. Dispersal and gene flow are key processes in both cases, linking demographic and evolutionary dynamics to each other, facilitating but also constraining the expansion of the current niche and the geographical range of species and determining the spatial scale and pattern of adaptation in heterogeneous environments. An eco-evolutionary metapopulation model helps explain the contrasting responses of species to habitat loss and fragmentation. Eco-evolutionary dynamics may facilitate the persistence of species in changing environments, but typically the evolutionary response only partially compensates for the negative ecological consequences of adverse environmental changes. PMID:22335524
Transition matrix model for evolutionary game dynamics
Ermentrout, G. Bard; Griffin, Christopher; Belmonte, Andrew
2016-03-01
We study an evolutionary game model based on a transition matrix approach, in which the total change in the proportion of a population playing a given strategy is summed directly over contributions from all other strategies. This general approach combines aspects of the traditional replicator model, such as preserving unpopulated strategies, with mutation-type dynamics, which allow for nonzero switching to unpopulated strategies, in terms of a single transition function. Under certain conditions, this model yields an endemic population playing non-Nash-equilibrium strategies. In addition, a Hopf bifurcation with a limit cycle may occur in the generalized rock-scissors-paper game, unlike the replicator equation. Nonetheless, many of the Folk Theorem results are shown to hold for this model.
Evolutionary dynamics with fluctuating population sizes and strong mutualism
Chotibut, Thiparat; Nelson, David R.
2015-08-01
Game theory ideas provide a useful framework for studying evolutionary dynamics in a well-mixed environment. This approach, however, typically enforces a strictly fixed overall population size, deemphasizing natural growth processes. We study a competitive Lotka-Volterra model, with number fluctuations, that accounts for natural population growth and encompasses interaction scenarios typical of evolutionary games. We show that, in an appropriate limit, the model describes standard evolutionary games with both genetic drift and overall population size fluctuations. However, there are also regimes where a varying population size can strongly influence the evolutionary dynamics. We focus on the strong mutualism scenario and demonstrate that standard evolutionary game theory fails to describe our simulation results. We then analytically and numerically determine fixation probabilities as well as mean fixation times using matched asymptotic expansions, taking into account the population size degree of freedom. These results elucidate the interplay between population dynamics and evolutionary dynamics in well-mixed systems.
Evolutionary dynamics with fluctuating population sizes and strong mutualism.
Chotibut, Thiparat; Nelson, David R
2015-08-01
Game theory ideas provide a useful framework for studying evolutionary dynamics in a well-mixed environment. This approach, however, typically enforces a strictly fixed overall population size, deemphasizing natural growth processes. We study a competitive Lotka-Volterra model, with number fluctuations, that accounts for natural population growth and encompasses interaction scenarios typical of evolutionary games. We show that, in an appropriate limit, the model describes standard evolutionary games with both genetic drift and overall population size fluctuations. However, there are also regimes where a varying population size can strongly influence the evolutionary dynamics. We focus on the strong mutualism scenario and demonstrate that standard evolutionary game theory fails to describe our simulation results. We then analytically and numerically determine fixation probabilities as well as mean fixation times using matched asymptotic expansions, taking into account the population size degree of freedom. These results elucidate the interplay between population dynamics and evolutionary dynamics in well-mixed systems. PMID:26382443
Evolutionary dynamics of nationalism and migration
Barreira da Silva Rocha, André
2013-08-01
I present a dynamic evolutionary game model to address the relation between nationalism against immigrants and assimilation of the latter into the host country culture. I assume a country composed of two different large polymorphic populations, one of native citizens and the other of immigrants. A native citizen may behave nationalistically or may welcome immigrants. Immigrants may have an interest in learning the host country language or not. Evolution is modeled using replicator dynamics (RD). I also account for the presence of an enclave of immigrants in the host country. In the RD, the latter represents the immigrants’ own population effect, which contribution to fitness is controlled using a parameter ρ, 0≤ρ≤1, that represents the enclave size. In line with the empirical literature on migration, the existence of an enclave of immigrants makes assimilation less likely to occur. For large values of ρ, complete assimilation may not occur even if immigrants and natives share very close cultures and norms. Government policy regarding nationalism is modeled both exogenously and endogenously. A single or multiple asymptotically stable states exist for all cases studied but one in which the dynamics is similar to that found in the predator-prey model of Lotka-Volterra for competing species.
无
2011-01-01
Compared with ellipse cavity, the spoke cavity has many advantages, especially for the low and medium beam energy. It will be used in the superconductor accelerator popular in the future. Based on the spoke cavity, we design and calculate an accelerator
Quantifying Slow Evolutionary Dynamics in RNA Fitness Landscapes
Sulc, P; Wagner, A.; Martin, O. C.
2009-01-01
We re-examine the evolutionary dynamics of RNA secondary structures under directional selection towards an optimum RNA structure. We find that the punctuated equilibria lead to a very slow approach to the optimum, following on average an inverse power of the evolutionary time. In addition, our study of the trajectories shows that the out-of-equilibrium effects due to the evolutionary process are very weak. In particular, the distribution of genotypes is close to that arising during equilibriu...
Dynamical characteristics of software trustworthiness and their evolutionary complexity
ZHENG ZhiMing; MA ShiLong; LI Wei; WEI Wei; JIANG Xin; ZHANG ZhanLi; GUO BingHui
2009-01-01
Developing trusted $oftwares has become an important trend and a natural choice In the development of software technology and applications, and software trustworthiness modeling has become a prerequisite and necessary means. To discuss and explain the basic scientific problems in software trustworthiness and to establish theoretical foundations for software trustworthiness measurement, combining the Ideas of dynamical system study, this paper studies evolutionary laws of software trustworthiness and the dynamical mechanism under the effect of various internal and external factors, and proposes dynamical models for software trustworthiness, thus, software trustworthiness can be considered as the statistical characteristics of behaviors of software systems in the dynamical and open environment. By analyzing two simple examples, the paper explains the relationship between the limit evolutionary behaviors of software trustworthiness attributes and dynamical system characteristics, and interprets the dynamical characteristics of software trustworthiness and their evolutionary complexity.
Eco-evolutionary spatial dynamics in the Glanville fritillary butterfly.
Hanski, Ilkka A
2011-08-30
Demographic population dynamics, gene flow, and local adaptation may influence each other and lead to coupling of ecological and evolutionary dynamics, especially in species inhabiting fragmented heterogeneous environments. Here, I review long-term research on eco-evolutionary spatial dynamics in the Glanville fritillary butterfly inhabiting a large network of approximately 4,000 meadows in Finland. The metapopulation persists in a balance between frequent local extinctions and recolonizations. The genetic spatial structure as defined by neutral markers is much more coarse-grained than the demographic spatial structure determined by the fragmented habitat, yet small-scale spatial structure has important consequences for the dynamics. I discuss three examples of eco-evolutionary spatial dynamics. (i) Extinction-colonization metapopulation dynamics influence allele frequency changes in the phosphoglucose isomerase (Pgi) gene, which leads to strong associations between genetic variation in Pgi and dispersal, recolonization, and local population dynamics. (ii) Inbreeding in local populations increases their risk for extinction, whereas reciprocal effects between inbreeding, population size, and emigration represent likely eco-evolutionary feedbacks. (iii) Genetically determined female oviposition preference for two host plant species exhibits a cline paralleling a gradient in host plant relative abundances, and host plant preference of dispersing females in relation to the host plant composition of habitat patches influences immigration (gene flow) and recolonization (founder events). Eco-evolutionary spatial dynamics in heterogeneous environments may not lead to directional evolutionary changes unless the environment itself changes, but eco-evolutionary dynamics may contribute to the maintenance of genetic variation attributable to fluctuating selection in space and time. PMID:21788506
An evolutionary dynamics model adapted to eusocial insects.
Louise van Oudenhove
Full Text Available This study aims to better understand the evolutionary processes allowing species coexistence in eusocial insect communities. We develop a mathematical model that applies adaptive dynamics theory to the evolutionary dynamics of eusocial insects, focusing on the colony as the unit of selection. The model links long-term evolutionary processes to ecological interactions among colonies and seasonal worker production within the colony. Colony population dynamics is defined by both worker production and colony reproduction. Random mutations occur in strategies, and mutant colonies enter the community. The interactions of colonies at the ecological timescale drive the evolution of strategies at the evolutionary timescale by natural selection. This model is used to study two specific traits in ants: worker body size and the degree of collective foraging. For both traits, trade-offs in competitive ability and other fitness components allows to determine conditions in which selection becomes disruptive. Our results illustrate that asymmetric competition underpins diversity in ant communities.
Quantifying slow evolutionary dynamics in RNA fitness landscapes.
Sulc, Petr; Wagner, Andreas; Martin, Olivier C
2010-12-01
We re-examine the evolutionary dynamics of RNA secondary structures under directional selection towards an optimum RNA structure. We find that the punctuated equilibria lead to a very slow approach to the optimum, following on average an inverse power of the evolutionary time. In addition, our study of the trajectories shows that the out-of-equilibrium effects due to the evolutionary process are very weak. In particular, the distribution of genotypes is close to that arising during equilibrium stabilizing selection. As a consequence, the evolutionary dynamics leave almost no measurable out-of-equilibrium trace, only the transition genotypes (close to the border between different periods of stasis) have atypical mutational properties. PMID:21121025
The evolutionary rate dynamically tracks changes in HIV-1 epidemics
Maljkovic-berry, Irina [Los Alamos National Laboratory; Athreya, Gayathri [Los Alamos National Laboratory; Daniels, Marcus [Los Alamos National Laboratory; Bruno, William [Los Alamos National Laboratory; Korber, Bette [Los Alamos National Laboratory; Kuiken, Carla [Los Alamos National Laboratory; Ribeiro, Ruy M [Los Alamos National Laboratory
2009-01-01
Large-sequence datasets provide an opportunity to investigate the dynamics of pathogen epidemics. Thus, a fast method to estimate the evolutionary rate from large and numerous phylogenetic trees becomes necessary. Based on minimizing tip height variances, we optimize the root in a given phylogenetic tree to estimate the most homogenous evolutionary rate between samples from at least two different time points. Simulations showed that the method had no bias in the estimation of evolutionary rates and that it was robust to tree rooting and topological errors. We show that the evolutionary rates of HIV-1 subtype B and C epidemics have changed over time, with the rate of evolution inversely correlated to the rate of virus spread. For subtype B, the evolutionary rate slowed down and tracked the start of the HAART era in 1996. Subtype C in Ethiopia showed an increase in the evolutionary rate when the prevalence increase markedly slowed down in 1995. Thus, we show that the evolutionary rate of HIV-1 on the population level dynamically tracks epidemic events.
Evolutionary Dynamics of Chronic Myeloid Leukemia
Dingli, David; Traulsen, Arne; Lenaerts, Tom; Pacheco, Jorge M.
2010-01-01
Cancer is an evolutionary process that arises due to mutations and expands through the selection of clones with higher reproductive success that will outcompete their peers. Most tumors require many mutations to explain the cancer phenotype, making it difficult to identify the gene(s) that confer the reproductive fitness to the clone. Moreover, the impact of any oncogene is context dependent: it can increase the fitness of particular stages of cell differentiation but not other stages. In add...
HIV evolutionary dynamics within and among hosts.
Lemey, Philippe; Rambaut, Andrew; Pybus, Oliver G
2006-01-01
The HIV evolutionary processes continuously unfold, leaving a measurable footprint in viral gene sequences. A variety of statistical models and inference techniques have been developed to reconstruct the HIV evolutionary history and to investigate the population genetic processes that shape viral diversity. Remarkably different population genetic forces are at work within and among hosts. Population-level HIV phylogenies are mainly shaped by selectively neutral epidemiologic processes, implying that genealogy-based population genetic inference can be useful to study the HIV epidemic history. Such evolutionary analyses have shed light on the origins of HIV, and on the epidemic spread of viral variants in different geographic locations and in different populations. The HIV genealogies reconstructed from within-host sequences indicate the action of selection pressure. In addition, recombination has a significant impact on HIV genetic diversity. Accurately quantifying both the adaptation rate and the population recombination rate of HIV will contribute to a better understanding of immune escape and drug resistance. Characterizing the impact of HIV transmission on viral genetic diversity will be a key factor in reconciling the different population genetic processes within and among hosts. PMID:17078483
Spatial effect on stochastic dynamics of bistable evolutionary games
We consider the lifetimes of metastable states in bistable evolutionary games (coordination games), and examine how they are affected by spatial structure. A semiclassical approximation based on a path integral method is applied to stochastic evolutionary game dynamics with and without spatial structure, and the lifetimes of the metastable states are evaluated. It is shown that the population dependence of the lifetimes is qualitatively different in these two models. Our result indicates that spatial structure can accelerate the transitions between metastable states. (paper)
Evolutionary dynamics of group interactions on structured populations: A review
Perc, Matjaz; Szolnoki, Attila; Floría, Luis M; Moreno, Yamir; 10.1098/rsif.2012.0997
2013-01-01
Interactions among living organisms, from bacteria colonies to human societies, are inherently more complex than interactions among particles and nonliving matter. Group interactions are a particularly important and widespread class, representative of which is the public goods game. In addition, methods of statistical physics have proven valuable for studying pattern formation, equilibrium selection, and self-organisation in evolutionary games. Here we review recent advances in the study of evolutionary dynamics of group interactions on structured populations, including lattices, complex networks and coevolutionary models. We also compare these results with those obtained on well-mixed populations. The review particularly highlights that the study of the dynamics of group interactions, like several other important equilibrium and non-equilibrium dynamical processes in biological, economical and social sciences, benefits from the synergy between statistical physics, network science and evolutionary game theory...
Saakian, David B
2009-01-01
The parallel mutation-selection evolutionary dynamics, in which mutation and replication are independent events, is solved exactly in the case that the Malthusian fitnesses associated to the genomes are described by the Random Energy Model (REM) and by a ferromagnetic version of the REM. The solution method uses the mapping of the evolutionary dynamics into a quantum Ising chain in a transverse field and the Suzuki-Trotter formalism to calculate the transition probabilities between configurations at different times. We find that in the case of the REM landscape the dynamics can exhibit three distinct regimes: pure diffusion or stasis for short times, depending on the fitness of the initial configuration, and a spin-glass regime for large times. The dynamic transition between these dynamical regimes is marked by discontinuities in the mean-fitness as well as in the overlap with the initial reference sequence. The relaxation to equilibrium is described by an inverse time decay. In the ferromagnetic REM, we find...
Emergence of structured communities through evolutionary dynamics.
Shtilerman, Elad; Kessler, David A; Shnerb, Nadav M
2015-10-21
Species-rich communities, in which many competing species coexist in a single trophic level, are quite frequent in nature, but pose a formidable theoretical challenge. In particular, it is known that complex competitive systems become unstable and unfeasible when the number of species is large. Recently, many studies have attributed the stability of natural communities to the structure of the interspecific interaction network, yet the nature of such structures and the underlying mechanisms responsible for them remain open questions. Here we introduce an evolutionary model, based on the generic Lotka-Volterra competitive framework, from which a stable, structured, diverse community emerges spontaneously. The modular structure of the competition matrix reflects the phylogeny of the community, in agreement with the hierarchial taxonomic classification. Closely related species tend to have stronger niche overlap and weaker fitness differences, as opposed to pairs of species from different modules. The competitive-relatedness hypothesis and the idea of emergent neutrality are discussed in the context of this evolutionary model. PMID:26231415
Dynamic Ising model: reconstruction of evolutionary trees
An evolutionary tree is a cascade of bifurcations starting from a single common root, generating a growing set of daughter species as time goes by. ‘Species’ here is a general denomination for biological species, spoken languages or any other entity which evolves through heredity. From the N currently alive species within a clade, distances are measured through pairwise comparisons made by geneticists, linguists, etc. The larger is such a distance that, for a pair of species, the older is their last common ancestor. The aim is to reconstruct the previously unknown bifurcations, i.e. the whole clade, from knowledge of the N(N − 1)/2 quoted distances, which are taken for granted. A mechanical method is presented and its applicability is discussed. (paper)
Eco-evolutionary dynamics of dispersal in spatially heterogeneous environments.
Hanski, Ilkka; Mononen, Tommi
2011-10-01
Ecology Letters (2011) 14: 1025-1034 ABSTRACT: Evolutionary changes in natural populations are often so fast that the evolutionary dynamics may influence ecological population dynamics and vice versa. Here we construct an eco-evolutionary model for dispersal by combining a stochastic patch occupancy metapopulation model with a model for changes in the frequency of fast-dispersing individuals in local populations. We test the model using data on allelic variation in the gene phosphoglucose isomerase (Pgi), which is strongly associated with dispersal rate in the Glanville fritillary butterfly. Population-specific measures of immigration and extinction rates and the frequency of fast-dispersing individuals among the immigrants explained 40% of spatial variation in Pgi allele frequency among 97 local populations. The model clarifies the roles of founder events and gene flow in dispersal evolution and resolves a controversy in the literature about the consequences of habitat loss and fragmentation on the evolution of dispersal. PMID:21794053
Bridging Developmental Systems Theory and Evolutionary Psychology Using Dynamic Optimization
Frankenhuis, Willem E.; Panchanathan, Karthik; Clark Barrett, H.
2013-01-01
Interactions between evolutionary psychologists and developmental systems theorists have been largely antagonistic. This is unfortunate because potential synergies between the two approaches remain unexplored. This article presents a method that may help to bridge the divide, and that has proven fruitful in biology: dynamic optimization. Dynamic…
A stochastic evolutionary model for survival dynamics
Fenner, Trevor; Loizou, George
2014-01-01
The recent interest in human dynamics has led researchers to investigate the stochastic processes that explain human behaviour in different contexts. Here we propose a generative model to capture the essential dynamics of survival analysis, traditionally employed in clinical trials and reliability analysis in engineering. In our model, the only implicit assumption made is that the longer an actor has been in the system, the more likely it is to have failed. We derive a power-law distribution for the process and provide preliminary empirical evidence for the validity of the model from two well-known survival analysis data sets.
Evolutionary dynamics of a smoothed war of attrition game.
Iyer, Swami; Killingback, Timothy
2016-05-01
In evolutionary game theory the War of Attrition game is intended to model animal contests which are decided by non-aggressive behavior, such as the length of time that a participant will persist in the contest. The classical War of Attrition game assumes that no errors are made in the implementation of an animal׳s strategy. However, it is inevitable in reality that such errors must sometimes occur. Here we introduce an extension of the classical War of Attrition game which includes the effect of errors in the implementation of an individual׳s strategy. This extension of the classical game has the important feature that the payoff is continuous, and as a consequence admits evolutionary behavior that is fundamentally different from that possible in the original game. We study the evolutionary dynamics of this new game in well-mixed populations both analytically using adaptive dynamics and through individual-based simulations, and show that there are a variety of possible outcomes, including simple monomorphic or dimorphic configurations which are evolutionarily stable and cannot occur in the classical War of Attrition game. In addition, we study the evolutionary dynamics of this extended game in a variety of spatially and socially structured populations, as represented by different complex network topologies, and show that similar outcomes can also occur in these situations. PMID:26903203
Temporal and evolutionary dynamics of two-component signaling pathways.
Salazar, Michael E; Laub, Michael T
2015-04-01
Bacteria sense and respond to numerous environmental signals through two-component signaling pathways. Typically, a given stimulus will activate a sensor histidine kinase to autophosphorylate and then phosphotransfer to a cognate response regulator, which can mount an appropriate response. Although these signaling pathways often appear to be simple switches, they can also orchestrate surprisingly sophisticated and complex responses. These temporal dynamics arise from several key regulatory features, including the bifunctionality of histidine kinases as well as positive and negative feedback loops. Two-component signaling pathways are also dynamic on evolutionary time-scales, expanding dramatically in many species through gene duplication and divergence. Here, we review recent work probing the temporal and evolutionary dynamics of two-component signaling systems. PMID:25589045
Evolutionary Dynamics of the World Wide Web
Bernardo A. Huberman; Adamic, Lada A.
1999-01-01
We present a theory for the growth dynamics of the World Wide Web that takes into account the wide range of stochastic growth rates in the number of pages per site, as well as the fact that new sites are created at different times. This leads to the prediction of a universal power law in the distribution of the number of pages per site which we confirm experimentally by analyzing data from large crawls made by the search engines Alexa and Infoseek. The existence of this power law not only imp...
Reconstructing the Nonlinear Dynamical Systems by Evolutionary Computation Techniques
LIU Minzhong; KANG Lishan
2006-01-01
We introduce a new dynamical evolutionary algorithm(DEA) based on the theory of statistical mechanics and investigate the reconstruction problem for the nonlinear dynamical systems using observation data. The convergence of the algorithm is discussed. We make the numerical experiments and test our model using the two famous chaotic systems (mainly the Lorenz and Chen systems ). The results show the relatively accurate reconstruction of these chaotic systems based on observational data can be obtained. Therefore we may conclude that there are broad prospects using our method to model the nonlinear dynamical systems.
Adversarial Scheduling in Evolutionary Game Dynamics
Istrate, Gabriel; Ravi, S S
2008-01-01
Consider a system in which players at nodes of an underlying graph G repeatedly play Prisoner's Dilemma against their neighbors. The players adapt their strategies based on the past behavior of their opponents by applying the so-called win-stay lose-shift strategy. This dynamics has been studied in (Kittock 94), (Dyer et al. 2002), (Mossel and Roch, 2006). With random scheduling, starting from any initial configuration with high probability the system reaches the unique fixed point in which all players cooperate. This paper investigates the validity of this result under various classes of adversarial schedulers. Our results can be sumarized as follows: 1. An adversarial scheduler that can select both participants to the game can preclude the system from reaching the unique fixed point on most graph topologies. 2. A nonadaptive scheduler that is only allowed to choose one of the participants is no more powerful than a random scheduler. With this restriction even an adaptive scheduler is not significantly more ...
Evolutionary dynamics of collective index insurance.
Pacheco, Jorge M; Santos, Francisco C; Levin, Simon A
2016-03-01
Index-based insurances offer promising opportunities for climate-risk investments in developing countries. Indeed, contracts conditional on, e.g., weather or livestock indexes can be cheaper to set up than conventional indemnity-based insurances, while offering a safety net to vulnerable households, allowing them to eventually escape poverty traps. Moreover, transaction costs by insurance companies may be additionally reduced if contracts, instead of arranged with single households, are endorsed by collectives of households that bear the responsibility of managing the division of the insurance coverage by its members whenever the index is surpassed, allowing for additional flexibility in what concerns risk-sharing and also allowing insurance companies to avoid the costs associated with moral hazard. Here we resort to a population dynamics framework to investigate under which conditions household collectives may find collective index insurances attractive, when compared with individual index insurances. We assume risk sharing among the participants of each collective, and model collective action in terms of an N-person threshold game. Compared to less affordable individual index insurances, we show how collective index insurances lead to a coordination problem in which the adoption of index insurances may become the optimal decision, spreading index insurance coverage to the entire population. We further investigate the role of risk-averse and risk-prone behaviors, as well as the role of partial correlation between insurance coverage and actual loss of crops, and in which way these affect the original coordination thresholds. PMID:26486802
Evolutionary history of Pacific salmon in dynamic environments
Waples, Robin S.; George R. Pess; Beechie, Tim
2008-01-01
Contemporary evolution of Pacific salmon (Oncorhynchus spp.) is best viewed in the context of the evolutionary history of the species and the dynamic ecosystems they inhabit. Speciation was complete by the late Miocene, leaving c. six million years for intraspecific diversification. Following the most recent glacial maximum, large areas became available for recolonization. Current intraspecific diversity is thus the product of recent evolution overlaid onto divergent historical lineages forge...
Evolutionary dynamics of bacteria in a human host environment
Lei YANG; Jelsbak, Lars; Marvig, Rasmus Lykke; Damkiær, Søren; Workman, Christopher T; Rau, Martin Holm; Hansen, Susse Kirkelund; Folkesson, Anders; Johansen, Helle Krogh; Ciofu, Oana; Høiby, Niels; Morten O. A. Sommer; Molin, Søren
2011-01-01
Laboratory evolution experiments have led to important findings relating organism adaptation and genomic evolution. However, continuous monitoring of long-term evolution has been lacking for natural systems, limiting our understanding of these processes in situ. Here we characterize the evolutionary dynamics of a lineage of a clinically important opportunistic bacterial pathogen, Pseudomonas aeruginosa, as it adapts to the airways of several individual cystic fibrosis patients over 200,000 ba...
Evolutionary dynamics in a simple model of self-assembly
Johnston, Iain G.; Ahnert, Sebastian E.; Doye, Jonathan P. K.; Louis, Ard A.
2011-06-01
We investigate the evolutionary dynamics of an idealized model for the robust self-assembly of two-dimensional structures called polyominoes. The model includes rules that encode interactions between sets of square tiles that drive the self-assembly process. The relationship between the model’s rule set and its resulting self-assembled structure can be viewed as a genotype-phenotype map and incorporated into a genetic algorithm. The rule sets evolve under selection for specified target structures. The corresponding complex fitness landscape generates rich evolutionary dynamics as a function of parameters such as the population size, search space size, mutation rate, and method of recombination. Furthermore, these systems are simple enough that in some cases the associated model genome space can be completely characterized, shedding light on how the evolutionary dynamics depends on the detailed structure of the fitness landscape. Finally, we apply the model to study the emergence of the preference for dihedral over cyclic symmetry observed for homomeric protein tetramers.
An Evolutionary Dynamic Clustering based Colour Image Segmentation
Amiya Halder, Nilvra Pathak
2011-02-01
Full Text Available We have presented a novel Dynamic Colour Image Segmentation (DCISSystem for colour image. In this paper, we have proposed an efficient colourimage segmentation algorithm based on evolutionary approach i.e. dynamic GAbased clustering (GADCIS. The proposed technique automatically determinesthe optimum number of clusters for colour images. The optimal number ofclusters is obtained by using cluster validity criterion with the help of Gaussiandistribution. The advantage of this method is that no a priori knowledge isrequired to segment the color image. The proposed algorithm is evaluated onwell known natural images and its performance is compared to other clusteringtechniques. Experimental results show the performance of the proposedalgorithm producing comparable segmentation results.
DeLong, John P; Gibert, Jean P
2016-02-01
Heritable trait variation is a central and necessary ingredient of evolution. Trait variation also directly affects ecological processes, generating a clear link between evolutionary and ecological dynamics. Despite the changes in variation that occur through selection, drift, mutation, and recombination, current eco-evolutionary models usually fail to track how variation changes through time. Moreover, eco-evolutionary models assume fitness functions for each trait and each ecological context, which often do not have empirical validation. We introduce a new type of model, Gillespie eco-evolutionary models (GEMs), that resolves these concerns by tracking distributions of traits through time as eco-evolutionary dynamics progress. This is done by allowing change to be driven by the direct fitness consequences of model parameters within the context of the underlying ecological model, without having to assume a particular fitness function. GEMs work by adding a trait distribution component to the standard Gillespie algorithm - an approach that models stochastic systems in nature that are typically approximated through ordinary differential equations. We illustrate GEMs with the Rosenzweig-MacArthur consumer-resource model. We show not only how heritable trait variation fuels trait evolution and influences eco-evolutionary dynamics, but also how the erosion of variation through time may hinder eco-evolutionary dynamics in the long run. GEMs can be developed for any parameter in any ordinary differential equation model and, furthermore, can enable modeling of multiple interacting traits at the same time. We expect GEMs will open the door to a new direction in eco-evolutionary and evolutionary modeling by removing long-standing modeling barriers, simplifying the link between traits, fitness, and dynamics, and expanding eco-evolutionary treatment of a greater diversity of ecological interactions. These factors make GEMs much more than a modeling advance, but an important
Iterative Dynamic Diversity Evolutionary Algorithm for Constrained Optimization
GAO Wei-Shang; SHAO Cheng
2014-01-01
Evolutionary algorithms (EAs) were shown to be effective for complex constrained optimization problems. However, inflexible exploration in general EAs would lead to losing the global optimum nearby the ill-convergence regions. In this paper, we propose an iterative dynamic diversity evolutionary algorithm (IDDEA) with contractive subregions guiding exploitation through local extrema to the global optimum in suitable steps. In IDDEA, a novel optimum estimation strategy with multi-agents evolving diversely is suggested to eﬃciently compute dominance trend and establish a subregion. In addition, a subregion converging iteration is designed to redistrict a smaller subregion in current subregion for next iteration, which is based on a special dominance estimation scheme. Meanwhile, an infimum penalty function is embedded into IDDEA to judge agents and penalize adaptively the unfeasible agents with the lowest fitness of feasible agents. Furthermore, several engineering design optimization problems taken from the specialized literature are successfully solved by the present algorithm with high reliable solutions.
Dynamics Calculation of Travel Wave Tube
无
2011-01-01
During the dynamics calculating of the travel tube, we must obtain the field map in the tube. The field map can be affected by not only the beam loading, but also the attenuation coefficient. The calculation of the attenuation coefficient
Evolutionary network games: Imitation and Best-Response dynamics
Cimini, Giulio; Sánchez, Angel
2014-01-01
We consider games of strategic substitutes and strategic complements on networks. We introduce two different evolutionary dynamics in order to refine their multiplicity of equilibria, that can be related to alternative informational contexts. We find that for the best-shot game, taken as a model for substitutes, a replicator-like dynamics does not lead to Nash equilibria, whereas it leads to unique equilibria (full cooperation or full defection, depending on the initial condition and the game parameter) for complements, represented by a coordination game. On the other hand, when the dynamics becomes more cognitively demanding in the form of a best response evolution, predictions are always Nash equilibria (at least when individuals are fully rational): For the best-shot game we find equilibria with a definite value of the fraction of contributors, whereas for the coordination game symmetric equilibria arise only for low or high initial fractions of cooperators. We also consider extensions of the natural incom...
Evolutionary dynamics of time-resolved social interactions
Cardillo, Alessio; Nicosia, Vincenzo; Sinatra, Roberta; Gómez-Gardeñes, Jesús; Latora, Vito
2013-01-01
Cooperation among unrelated individuals is frequently observed in social groups when their members join efforts and resources to obtain a shared benefit which is unachievable by singles. However, understanding why cooperation arises despite the natural tendency of individuals towards selfish behaviors is still an open problem and represents one of the most fascinating challenges in volutionary dynamics. Very recently, the structural characterization of the networks upon which social interactions take place has shed some light on the mechanisms by which cooperative behaviours emerge and eventually overcome the individual temptation to defect. In particular, it has been found that the heterogeneity in the number of social ties and the presence of tightly-knit communities lead to a significant increase of cooperation as compared with the unstructured and homogeneous connection patterns considered in classical evolutionary dynamics. Here we investigate the role of social ties dynamics for the emergence of coopera...
Quantifying the Determinants of Evolutionary Dynamics Leading to Drug Resistance.
Guillaume Chevereau
Full Text Available The emergence of drug resistant pathogens is a serious public health problem. It is a long-standing goal to predict rates of resistance evolution and design optimal treatment strategies accordingly. To this end, it is crucial to reveal the underlying causes of drug-specific differences in the evolutionary dynamics leading to resistance. However, it remains largely unknown why the rates of resistance evolution via spontaneous mutations and the diversity of mutational paths vary substantially between drugs. Here we comprehensively quantify the distribution of fitness effects (DFE of mutations, a key determinant of evolutionary dynamics, in the presence of eight antibiotics representing the main modes of action. Using precise high-throughput fitness measurements for genome-wide Escherichia coli gene deletion strains, we find that the width of the DFE varies dramatically between antibiotics and, contrary to conventional wisdom, for some drugs the DFE width is lower than in the absence of stress. We show that this previously underappreciated divergence in DFE width among antibiotics is largely caused by their distinct drug-specific dose-response characteristics. Unlike the DFE, the magnitude of the changes in tolerated drug concentration resulting from genome-wide mutations is similar for most drugs but exceptionally small for the antibiotic nitrofurantoin, i.e., mutations generally have considerably smaller resistance effects for nitrofurantoin than for other drugs. A population genetics model predicts that resistance evolution for drugs with this property is severely limited and confined to reproducible mutational paths. We tested this prediction in laboratory evolution experiments using the "morbidostat", a device for evolving bacteria in well-controlled drug environments. Nitrofurantoin resistance indeed evolved extremely slowly via reproducible mutations-an almost paradoxical behavior since this drug causes DNA damage and increases the mutation
A quantitative evolutionary theory of adaptive behavior dynamics.
McDowell, J J
2013-10-01
The idea that behavior is selected by its consequences in a process analogous to organic evolution has been discussed for over 100 years. A recently proposed theory instantiates this idea by means of a genetic algorithm that operates on a population of potential behaviors. Behaviors in the population are represented by numbers in decimal integer (phenotypic) and binary bit string (genotypic) forms. One behavior from the population is emitted at random each time tick, after which a new population of potential behaviors is constructed by recombining parent behavior bit strings. If the emitted behavior produced a benefit to the organism, then parents are chosen on the basis of their phenotypic similarity to the emitted behavior; otherwise, they are chosen at random. After parent behavior recombination, the population is subjected to a small amount of mutation by flipping random bits in the population's bit strings. The behavior generated by this process of selection, reproduction, and mutation reaches equilibrium states that conform to every empirically valid equation of matching theory, exactly and without systematic error. These equations are known to describe the behavior of many vertebrate species, including humans, in a variety of experimental, naturalistic, natural, and social environments. The evolutionary theory also generates instantaneous dynamics and patterns of preference change in constantly changing environments that are consistent with the dynamics of live-organism behavior. These findings support the assertion that the world of behavior we observe and measure is generated by evolutionary dynamics. PMID:24219847
Evolutionary dynamics of N-person snowdrift game
Sui, Xiukai; Cong, Rui; Li, Kun; Wang, Long
2015-12-01
In this letter, we investigate the evolutionary dynamics of N-person snowdrift game in both well-mixed and structured populations. For well-mixed populations, we construct a double-threshold model considering both the necessary and the minimum cost players should pay for completing the task. We have explored the influences of these thresholds on both equilibrium points in infinite populations and the fixation probabilities in finite populations. Results present complicated behaviors that show characteristics of both stag-hunt game and snowdrift game. For structured populations, we use pair approximation and diffusion approximation to derive the critical benefit-to-cost ratio in favor of cooperation.
The Dynamic Evolutionary History of Pancrustacean Eyes and Opsins.
Henze, Miriam J; Oakley, Todd H
2015-11-01
Pancrustacea (Hexapoda plus Crustacea) display an enormous diversity of eye designs, including multiple types of compound eyes and single-chambered eyes, often with color vision and/or polarization vision. Although the eyes of some pancrustaceans are well-studied, there is still much to learn about the evolutionary paths to this amazing visual diversity. Here, we examine the evolutionary history of eyes and opsins across the principle groups of Pancrustacea. First, we review the distribution of lateral and median eyes, which are found in all major pancrustacean clades (Oligostraca, Multicrustacea, and Allotriocarida). At the same time, each of those three clades has taxa that lack lateral and/or median eyes. We then compile data on the expression of visual r-opsins (rhabdomeric opsins) in lateral and median eyes across Pancrustacea and find no evidence for ancient opsin clades expressed in only one type of eye. Instead, opsin clades with eye-specific expression are products of recent gene duplications, indicating a dynamic past, during which opsins often changed expression from one type of eye to another. We also investigate the evolutionary history of peropsins and r-opsins, which are both known to be expressed in eyes of arthropods. By searching published transcriptomes, we discover for the first time crustacean peropsins and suggest that previously reported odonate opsins may also be peropsins. Finally, from analyzing a reconciled, phylogenetic tree of arthropod r-opsins, we infer that the ancestral pancrustacean had four visual opsin genes, which we call LW2, MW1, MW2, and SW. These are the progenitors of opsin clades that later were variously duplicated or lost during pancrustacean evolution. Together, our results reveal a particularly dynamic history, with losses of eyes, duplication and loss of opsin genes, and changes in opsin expression between types of eyes. PMID:26319405
Bidirectional Dynamic Diversity Evolutionary Algorithm for Constrained Optimization
Weishang Gao
2013-01-01
Full Text Available Evolutionary algorithms (EAs were shown to be effective for complex constrained optimization problems. However, inflexible exploration-exploitation and improper penalty in EAs with penalty function would lead to losing the global optimum nearby or on the constrained boundary. To determine an appropriate penalty coefficient is also difficult in most studies. In this paper, we propose a bidirectional dynamic diversity evolutionary algorithm (Bi-DDEA with multiagents guiding exploration-exploitation through local extrema to the global optimum in suitable steps. In Bi-DDEA potential advantage is detected by three kinds of agents. The scale and the density of agents will change dynamically according to the emerging of potential optimal area, which play an important role of flexible exploration-exploitation. Meanwhile, a novel double optimum estimation strategy with objective fitness and penalty fitness is suggested to compute, respectively, the dominance trend of agents in feasible region and forbidden region. This bidirectional evolving with multiagents can not only effectively avoid the problem of determining penalty coefficient but also quickly converge to the global optimum nearby or on the constrained boundary. By examining the rapidity and veracity of Bi-DDEA across benchmark functions, the proposed method is shown to be effective.
Lattice Dynamics Calculation in MGB2
In Present report, We have introduced a new theoretical results for MgB2 by using home design programme Lattice Dynamics. we have calculated partial and total density of states (PDOS, TDOS), infrared and Raman spectrums and specific heat capacity. Dispersion curves in different symmetry points are calculated and found that there is agreement with other calculations. Also we have tried to investigate the Boron Isotope effect on the calculated properties
Evolutionary dynamics of plants and animals: a comparative approach
Valentine, J. W.; Tiffney, B. H.; Sepkoski, J. J. Jr; Sepkoski JJ, J. r. (Principal Investigator)
1991-01-01
Patterns of longevity and rate of appearance of taxa in the fossil record indicate a different evolutionary dynamic between land plants and marine invertebrates. Among marine invertebrates, rates of taxonomic turnover declined through the Phanerozoic, with increasingly extinction-resistant, long-lived, clades coming to dominate. Among terrestrial vascular plants, rates of turnover increased through the Phanerozoic, with short-lived, extinction-prone clades coming to dominate from the Devonian to the present. Terrestrial vertebrates appear to approximate the marine invertebrate pattern more closely than the plant record. We identify two features which individually or jointly may have influenced this distinction. First, land plants continuously invaded stressful environments during their evolution, while marine invertebrates and terrestrial vertebrates did not. Second, the relative structural simplicity and indeterminate mode of plant growth vs. the relative structural complexity and determinate mode of animal growth may have influenced the timing of major clade origin in the two groups.
The mathematical law of evolutionary information dynamics and an observer's evolution regularities
Lerner, Vladimir S
2011-01-01
An interactive stochastics, evaluated by an entropy functional (EF) of a random field and informational process' path functional (IPF), allows us modeling the evolutionary information processes and revealing regularities of evolution dynamics. Conventional Shannon's information measure evaluates a sequence of the process' static events for each information state and do not reveal hidden dynamic connections between these events. The paper formulates the mathematical forms of the information regularities, based on a minimax variation principle (VP) for IPF, applied to the evolution's both random microprocesses and dynamic macroprocesses. The paper shows that the VP single form of the mathematical law leads to the following evolutionary regularities: -creation of the order from stochastics through the evolutionary macrodynamics, described by a gradient of dynamic potential, evolutionary speed and the evolutionary conditions of a fitness and diversity; -the evolutionary hierarchy with growing information values a...
A New, Efficient Stellar Evolution Code for Calculating Complete Evolutionary Tracks
Kovetz, Attay; Prialnik, Dina
2008-01-01
We present a new stellar evolution code and a set of results, demonstrating its capability at calculating full evolutionary tracks for a wide range of masses and metallicities. The code is fast and efficient, and is capable of following through all evolutionary phases, without interruption or human intervention. It is meant to be used also in the context of modeling the evolution of dense stellar systems, for performing live calculations for both normal star models and merger-products. The code is based on a fully implicit, adaptive-grid numerical scheme that solves simultaneously for structure, mesh and chemical composition. Full details are given for the treatment of convection, equation of state, opacity, nuclear reactions and mass loss. Results of evolutionary calculations are shown for a solar model that matches the characteristics of the present sun to an accuracy of better than 1%; a $1 \\Msun$ model for a wide range of metallicities; a series of models of stellar populations I and II, for the mass rang...
Ribosome dynamics and the evolutionary history of ribosomes
Fox, George E.; Paci, Maxim; Tran, Quyen; Petrov, Anton S.; Williams, Loren D.
2015-09-01
The ribosome is a dynamic nanomachine responsible for coded protein synthesis. Its major subsystems were essentially in place at the time of the last universal common ancestor (LUCA). Ribosome evolutionary history thus potentially provides a window into the pre- LUCA world. This history begins with the origins of the peptidyl transferase center where the actual peptide is synthesized and then continues over an extended timeframe as additional functional centers including the GTPase center are added. The large ribosomal RNAs (rRNAs) have grown over time by an accretion process and a model exists that proposes a relative age of each accreted element. We have compared atomic resolution ribosome structures before and after EF-G bound GTP hydrolysis and thereby identified the location of 23 pivot points in the large rRNAs that facilitate ribosome dynamics. Pivots in small subunit helices h28 and h44 appear to be especially central to the process and according to the accretion model significantly older than the other helices containing pivots. Overall, the results suggest that ribosomal dynamics occurred in two phases. In the first phase, an inherently mobile h28/h44 combination provided the flexibility needed to create a dynamic ribosome that was essentially a Brownian machine. This addition likely made coded peptide synthesis possible by facilitating movement of a primitive mRNA. During the second phase, addition of pivoting elements and the creation of a factor binding site allowed the regulation of the inherent motion created by h28/h44. All of these events likely occurred before LUCA.
Evolutionary dynamics of general group interactions in structured populations
Li, Aming; Broom, Mark; Du, Jinming; Wang, Long
2016-02-01
The evolution of populations is influenced by many factors, and the simple classical models have been developed in a number of important ways. Both population structure and multiplayer interactions have been shown to significantly affect the evolution of important properties, such as the level of cooperation or of aggressive behavior. Here we combine these two key factors and develop the evolutionary dynamics of general group interactions in structured populations represented by regular graphs. The traditional linear and threshold public goods games are adopted as models to address the dynamics. We show that for linear group interactions, population structure can favor the evolution of cooperation compared to the well-mixed case, and we see that the more neighbors there are, the harder it is for cooperators to persist in structured populations. We further show that threshold group interactions could lead to the emergence of cooperation even in well-mixed populations. Here population structure sometimes inhibits cooperation for the threshold public goods game, where depending on the benefit to cost ratio, the outcomes are bistability or a monomorphic population of defectors or cooperators. Our results suggest, counterintuitively, that structured populations are not always beneficial for the evolution of cooperation for nonlinear group interactions.
Modeling Evolutionary Dynamics of Lurking in Social Networks
Javarone, Marco Alberto; Tagarelli, Andrea
2016-01-01
Lurking is a complex user-behavioral phenomenon that occurs in all large-scale online communities and social networks. It generally refers to the behavior characterizing users that benefit from the information produced by others in the community without actively contributing back to the production of social content. The amount and evolution of lurkers may strongly affect an online social environment, therefore understanding the lurking dynamics and identifying strategies to curb this trend are relevant problems. In this regard, we introduce the Lurker Game, i.e., a model for analyzing the transitions from a lurking to a non-lurking (i.e., active) user role, and vice versa, in terms of evolutionary game theory. We evaluate the proposed Lurker Game by arranging agents on complex networks and analyzing the system evolution, seeking relations between the network topology and the final equilibrium of the game. Results suggest that the Lurker Game is suitable to model the lurking dynamics, showing how the adoption ...
Nuclear Research Center IRT reactor dynamics calculation
The main features of the code DIRT, for dynamical calculations are described in the paper. With the results obtained by the program, an analysis of the dynamic behaviour of the Research Reactor IRT of the Nuclear Research Center (CIN) is performed. Different transitories were considered such as variation of the system reactivity, coolant inlet temperature variation and also variations of the coolant velocity through the reactor core. 3 refs
Evolutionary game dynamics of controlled and automatic decision-making
Toupo, Danielle F. P.; Strogatz, Steven H.; Cohen, Jonathan D.; Rand, David G.
2015-07-01
We integrate dual-process theories of human cognition with evolutionary game theory to study the evolution of automatic and controlled decision-making processes. We introduce a model in which agents who make decisions using either automatic or controlled processing compete with each other for survival. Agents using automatic processing act quickly and so are more likely to acquire resources, but agents using controlled processing are better planners and so make more effective use of the resources they have. Using the replicator equation, we characterize the conditions under which automatic or controlled agents dominate, when coexistence is possible and when bistability occurs. We then extend the replicator equation to consider feedback between the state of the population and the environment. Under conditions in which having a greater proportion of controlled agents either enriches the environment or enhances the competitive advantage of automatic agents, we find that limit cycles can occur, leading to persistent oscillations in the population dynamics. Critically, however, these limit cycles only emerge when feedback occurs on a sufficiently long time scale. Our results shed light on the connection between evolution and human cognition and suggest necessary conditions for the rise and fall of rationality.
Computational complexity of ecological and evolutionary spatial dynamics.
Ibsen-Jensen, Rasmus; Chatterjee, Krishnendu; Nowak, Martin A
2015-12-22
There are deep, yet largely unexplored, connections between computer science and biology. Both disciplines examine how information proliferates in time and space. Central results in computer science describe the complexity of algorithms that solve certain classes of problems. An algorithm is deemed efficient if it can solve a problem in polynomial time, which means the running time of the algorithm is a polynomial function of the length of the input. There are classes of harder problems for which the fastest possible algorithm requires exponential time. Another criterion is the space requirement of the algorithm. There is a crucial distinction between algorithms that can find a solution, verify a solution, or list several distinct solutions in given time and space. The complexity hierarchy that is generated in this way is the foundation of theoretical computer science. Precise complexity results can be notoriously difficult. The famous question whether polynomial time equals nondeterministic polynomial time (i.e., P = NP) is one of the hardest open problems in computer science and all of mathematics. Here, we consider simple processes of ecological and evolutionary spatial dynamics. The basic question is: What is the probability that a new invader (or a new mutant) will take over a resident population? We derive precise complexity results for a variety of scenarios. We therefore show that some fundamental questions in this area cannot be answered by simple equations (assuming that P is not equal to NP). PMID:26644569
Behavioral variability in an evolutionary theory of behavior dynamics.
Popa, Andrei; McDowell, J J
2016-03-01
McDowell's evolutionary theory of behavior dynamics (McDowell, 2004) instantiates populations of behaviors (abstractly represented by integers) that evolve under the selection pressure of the environment in the form of positive reinforcement. Each generation gives rise to the next via low-level Darwinian processes of selection, recombination, and mutation. The emergent patterns can be analyzed and compared to those produced by biological organisms. The purpose of this project was to explore the effects of high mutation rates on behavioral variability in environments that arranged different reinforcer rates and magnitudes. Behavioral variability increased with the rate of mutation. High reinforcer rates and magnitudes reduced these effects; low reinforcer rates and magnitudes augmented them. These results are in agreement with live-organism research on behavioral variability. Various combinations of mutation rates, reinforcer rates, and reinforcer magnitudes produced similar high-level outcomes (equifinality). These findings suggest that the independent variables that describe an experimental condition interact; that is, they do not influence behavior independently. These conclusions have implications for the interpretation of high levels of variability, mathematical undermatching, and the matching theory. The last part of the discussion centers on a potential biological counterpart for the rate of mutation, namely spontaneous fluctuations in the brain's default mode network. PMID:27002687
Evolutionary Voluntary Prisoner’s Dilemma Game under Deterministic and Stochastic Dynamics
Qian Yu
2015-03-01
Full Text Available The voluntary prisoner’s dilemma (VPD game has sparked interest from various fields since it was proposed as an effective mechanism to incentivize cooperative behavior. Current studies show that the inherent cyclic dominance of the strategies of the VPD game results in periodic oscillations in population. This paper investigated the influence of the level of individual rationality and the size of a population on the evolutionary dynamics of the VPD game. Different deterministic dynamics, such as the replicator dynamic, the Smith dynamic, the Brown-von Neumann-Nash (BNN dynamic and the best response (BR dynamic, for the evolutionary VPD game were modeled and simulated. The stochastic evolutionary dynamics based on quasi birth and death (QBD process was proposed for the evolutionary VPD game and compared with deterministic dynamics. The results indicated that with the increase of the loners’ fixed payoff, the loner is more likely to remain in the stable state of a VPD game under any of the dynamics mentioned above. However, the different speeds of motion under the dynamics in the cycle dominance proved to be diverse under different evolutionary dynamics and also highly sensitive to the rationality of individuals in a population. Furthermore, in QBD stochastic dynamics, the size of the population has a remarkable effect on the possibility distribution. When the population size increases, the limited distribution of the QBD process will be in accordance with the results in the deterministic dynamics.
Modeling evolutionary dynamics of epigenetic mutations in hierarchically organized tumors.
Andrea Sottoriva
2011-05-01
Full Text Available The cancer stem cell (CSC concept is a highly debated topic in cancer research. While experimental evidence in favor of the cancer stem cell theory is apparently abundant, the results are often criticized as being difficult to interpret. An important reason for this is that most experimental data that support this model rely on transplantation studies. In this study we use a novel cellular Potts model to elucidate the dynamics of established malignancies that are driven by a small subset of CSCs. Our results demonstrate that epigenetic mutations that occur during mitosis display highly altered dynamics in CSC-driven malignancies compared to a classical, non-hierarchical model of growth. In particular, the heterogeneity observed in CSC-driven tumors is considerably higher. We speculate that this feature could be used in combination with epigenetic (methylation sequencing studies of human malignancies to prove or refute the CSC hypothesis in established tumors without the need for transplantation. Moreover our tumor growth simulations indicate that CSC-driven tumors display evolutionary features that can be considered beneficial during tumor progression. Besides an increased heterogeneity they also exhibit properties that allow the escape of clones from local fitness peaks. This leads to more aggressive phenotypes in the long run and makes the neoplasm more adaptable to stringent selective forces such as cancer treatment. Indeed when therapy is applied the clone landscape of the regrown tumor is more aggressive with respect to the primary tumor, whereas the classical model demonstrated similar patterns before and after therapy. Understanding these often counter-intuitive fundamental properties of (non-hierarchically organized malignancies is a crucial step in validating the CSC concept as well as providing insight into the therapeutical consequences of this model.
Modeling viral evolutionary dynamics after telaprevir-based treatment.
Eric L Haseltine
2014-08-01
Full Text Available For patients infected with hepatitis C virus (HCV, the combination of the direct-acting antiviral agent telaprevir, pegylated-interferon alfa (Peg-IFN, and ribavirin (RBV significantly increases the chances of sustained virologic response (SVR over treatment with Peg-IFN and RBV alone. If patients do not achieve SVR with telaprevir-based treatment, their viral population is often significantly enriched with telaprevir-resistant variants at the end of treatment. We sought to quantify the evolutionary dynamics of these post-treatment resistant variant populations. Previous estimates of these dynamics were limited by analyzing only population sequence data (20% sensitivity, qualitative resistance information from 388 patients enrolled in Phase 3 clinical studies. Here we add clonal sequence analysis (5% sensitivity, quantitative for a subset of these patients. We developed a computational model which integrates both the qualitative and quantitative sequence data, and which forms a framework for future analyses of drug resistance. The model was qualified by showing that deep-sequence data (1% sensitivity from a subset of these patients are consistent with model predictions. When determining the median time for viral populations to revert to 20% resistance in these patients, the model predicts 8.3 (95% CI: 7.6, 8.4 months versus 10.7 (9.9, 12.8 months estimated using solely population sequence data for genotype 1a, and 1.0 (0.0, 1.4 months versus 0.9 (0.0, 2.7 months for genotype 1b. For each individual patient, the time to revert to 20% resistance predicted by the model was typically comparable to or faster than that estimated using solely population sequence data. Furthermore, the model predicts a median of 11.0 and 2.1 months after treatment failure for viral populations to revert to 99% wild-type in patients with HCV genotypes 1a or 1b, respectively. Our modeling approach provides a framework for projecting accurate, quantitative assessment of
Evolutionary aspects of gift-giving dynamics among Norwegian students
2007-01-01
Data from questionnaires filled out by 336 students during a nine day period in January and February 1999 at the University of Oslo, Norway was analyzed to find patterns in gift-giving behavior corresponding to predicted evolutionary biological and evolutionary psychological hypotheses. Gifts given and received, people given to and received from, monetary value of gifts given and estimated monetary value of gifts received were tallied. We tested the effects of four main factors: kin and non-k...
Bias in Dynamic Monte Carlo Alpha Calculations
Sweezy, Jeremy Ed [Los Alamos National Lab. (LANL), Los Alamos, NM (United States); Nolen, Steven Douglas [Los Alamos National Lab. (LANL), Los Alamos, NM (United States); Adams, Terry R. [Los Alamos National Lab. (LANL), Los Alamos, NM (United States); Trahan, Travis John [Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
2015-02-06
A 1/N bias in the estimate of the neutron time-constant (commonly denoted as α) has been seen in dynamic neutronic calculations performed with MCATK. In this paper we show that the bias is most likely caused by taking the logarithm of a stochastic quantity. We also investigate the known bias due to the particle population control method used in MCATK. We conclude that this bias due to the particle population control method is negligible compared to other sources of bias.
The co-evolutionary dynamics of directed network of spin market agents
Horvath, D.; Z. Kuscsik; Gmitra, M.
2005-01-01
The spin market model [S. Bornholdt, Int.J.Mod.Phys. C 12 (2001) 667] is extended into co-evolutionary version, where strategies of interacting and competitive traders are represented by local and global couplings between the nodes of dynamic directed stochastic network. The co-evolutionary principles are applied in the frame of Bak - Sneppen self-organized dynamics [P. Bak, K. Sneppen, Phys. Rev. Letter 71 (1993) 4083] that includes the processes of selection and extinction actuated by the l...
Eco-evolutionary dynamics in a coevolving host-virus system.
Frickel, Jens; Sieber, Michael; Becks, Lutz
2016-04-01
Eco-evolutionary dynamics have been shown to be important for understanding population and community stability and their adaptive potential. However, coevolution in the framework of eco-evolutionary theory has not been addressed directly. Combining experiments with an algal host and its viral parasite, and mathematical model analyses we show eco-evolutionary dynamics in antagonistic coevolving populations. The interaction between antagonists initially resulted in arms race dynamics (ARD) with selective sweeps, causing oscillating host-virus population dynamics. However, ARD ended and populations stabilised after the evolution of a general resistant host, whereas a trade-off between host resistance and growth then maintained host diversity over time (trade-off driven dynamics). Most importantly, our study shows that the interaction between ecology and evolution had important consequences for the predictability of the mode and tempo of adaptive change and for the stability and adaptive potential of populations. PMID:26898162
EVOLUTIONARY DYNAMIC MODEL OF POPULATION WITH NICHE CONSTRUCTION AND ITS APPLICATION RESEARCH
无
2006-01-01
Based on the theories and approaches in biomechanics, the mechanism and pattern of niche construction were discussed systematically. Through establishing the spatial pattern of niche and its measuring-fitness formula, and the dynamic system models of single- and two-population with niche construction, including corresponding theoretical analysis and numerical simulation on their evolutionary dynamics of population and the mechanism of competitive coexistence, the co-evolutionary relationship between organisms and their environments was revealed. The results indicate that population dynamics is governed by positive feedback between primary ecological factors and resource content.Niche construction generates an evolutionary effect in system by influencing the fitness of population. A threshold effect exists in single population dynamic system. In dynamic system of two competitive populations, niche construction can lead to alternative competitive consequences, which may be a potential mechanism to explain the competitive coexistence of species.
Design games: A conceptual framework for dynamic evolutionary design
Sönmez, N.O.; Erdem, A.
2014-01-01
Most evolutionary computation (EC) applications in design fields either assume simplified, static, performance-oriented procedures for design or focus on well-defined sub-problems, to be able to impose problem-solving and optimization schemes on design tasks, which render known EC techniques directl
Despotism, Democracy, and the Evolutionary Dynamics of Leadership and Followership
Van Vugt, Mark
2009-01-01
Responds to comments made by George B. Graen and Stephen J. Guastello on the current author's article Leadership, followership, and evolution: Some lessons from the past by Van Vugt, Hogan, and Kaiser. In the original article my co-authors and I proposed a new way of thinking about leadership, informed by evolutionary (neo-Darwinian) theory. In…
Electrostatic Energy Calculations for Molecular Dynamics
Love, M J; Comment, Henri J.F. Jansen; Love, Michael J.
1995-01-01
The evaluation of Coulomb forces is a difficult task. The summations that are involved converge only conditionally and care has to be taken in selecting the appropriate procedure to define the limits. The Ewald method is a standard method for obtaining Coulomb forces, but this method is rather slow, since it depends on the square of the number of atoms in a unit cell. In this paper we have adapted the plane-wise summation method for the evaluation of Coulomb forces. The use of this method allows for larger computational cells in molecular dynamics calculations.
Stochastic evolutionary dynamics of minimum-effort coordination games
Li, Kun; Wang, Long
2016-01-01
The minimum-effort coordination game, having potentially important implications in both evolutionary biology and sociology, draws recently more attention for the fact that human behavior in this social dilemma is often inconsistent with the predictions of classic game theory. In the framework of classic game theory, any common effort level is a strict and trembling hand perfect Nash equilibrium, so that no desideratum is provided for selecting among them. Behavior experiments, however, show that the effort levels employed by subjects are inversely related to the effort costs. Here, we combine coalescence theory and evolutionary game theory to investigate this game in finite populations. Both analytic results and individual-based simulations show that effort costs play a key role in the evolution of contribution levels, which is in good agreement with those observed experimentally. Besides well-mixed populations, set structured populations, where the population structure itself is a consequence of the evolutio...
Eco-evolutionary metapopulation dynamics and the spatial scale of adaptation.
Hanski, Ilkka; Mononen, Tommi; Ovaskainen, Otso
2011-01-01
We construct a model that combines extinction-colonization dynamics with the dynamics of local adaptation in a network of habitat patches of dissimilar qualities. We derive a deterministic approximation for the stochastic model that allows the calculation of patch-specific incidences of occupancy and levels of adaptation at steady state. Depending on (i) the strength of local selection, (ii) the amount of genetic variance, (iii) the demographic cost of maladaptation, (iv) the spatial scale of gene flow, and (v) the amount of habitat heterogeneity, the model predicts adaptation at different spatial scales. Local adaptation is predicted when there is much genetic variance and strong selection, while network-level adaptation occurs when the demographic cost of maladaptation is low. For little genetic variance and high cost of maladaptation, the model predicts network-level habitat specialization in species with long-range migration but an intermediate scale of adaptation (mosaic specialization) in species with short-range migration. In fragmented landscapes, the evolutionary dynamics of adaptation may both decrease and enhance metapopulation viability in comparison with no evolution. The model can be applied to real patch networks with given sizes, qualities, and spatial positions of habitat patches. PMID:21090992
The stability concept of evolutionary game theory a dynamic approach
1992-01-01
These Notes grew from my research in evolutionary biology, specifically on the theory of evolutionarily stable strategies (ESS theory), over the past ten years. Personally, evolutionary game theory has given me the opportunity to transfer my enthusiasm for abstract mathematics to more practical pursuits. I was fortunate to have entered this field in its infancy when many biologists recognized its potential but were not prepared to grant it general acceptance. This is no longer the case. ESS theory is now a rapidly expanding (in both applied and theoretical directions) force that no evolutionary biologist can afford to ignore. Perhaps, to continue the life-cycle metaphor, ESS theory is now in its late adolescence and displays much of the optimism and exuberance of this exciting age. There are dangers in writing a text about a theory at this stage of development. A comprehensive treatment would involve too many loose ends for the reader to appreciate the central message. On the other hand, the current central m...
An experimental investigation of evolutionary dynamics in the Rock-Paper-Scissors game
Moshe Hoffman; Sigrid Suetens; Uri Gneezy; Nowak, Martin A.
2015-01-01
Game theory describes social behaviors in humans and other biological organisms. By far, the most powerful tool available to game theorists is the concept of a Nash Equilibrium (NE), which is motivated by perfect rationality. NE specifies a strategy for everyone, such that no one would benefit by deviating unilaterally from his/her strategy. Another powerful tool available to game theorists are evolutionary dynamics (ED). Motivated by evolutionary and learning processes, ED specify changes in...
Female polymorphism, frequency dependence, and rapid evolutionary dynamics in natural populations
Svensson, Erik; Abbott, Jessica; Härdling, Roger
2005-01-01
Rapid evolutionary change over a few generations has been documented in natural populations. Such changes are observed as organisms invade new environments, and they are often triggered by changed interspecific interactions, such as differences in predation regimes. However, in spite of increased recognition of antagonistic male-female mating interactions, there is very limited evidence that such intraspecific interactions could cause rapid evolutionary dynamics in nature. This is because eco...
Improving the accuracy of dynamic mass calculation
Oleksandr F. Dashchenko
2015-06-01
Full Text Available With the acceleration of goods transporting, cargo accounting plays an important role in today's global and complex environment. Weight is the most reliable indicator of the materials control. Unlike many other variables that can be measured indirectly, the weight can be measured directly and accurately. Using strain-gauge transducers, weight value can be obtained within a few milliseconds; such values correspond to the momentary load, which acts on the sensor. Determination of the weight of moving transport is only possible by appropriate processing of the sensor signal. The aim of the research is to develop a methodology for weighing freight rolling stock, which increases the accuracy of the measurement of dynamic mass, in particular wagon that moves. Apart from time-series methods, preliminary filtration for improving the accuracy of calculation is used. The results of the simulation are presented.
Stochastic evolutionary dynamics in minimum-effort coordination games
Li, Kun; Cong, Rui; Wang, Long
2016-08-01
The minimum-effort coordination game draws recently more attention for the fact that human behavior in this social dilemma is often inconsistent with the predictions of classical game theory. Here, we combine evolutionary game theory and coalescence theory to investigate this game in finite populations. Both analytic results and individual-based simulations show that effort costs play a key role in the evolution of contribution levels, which is in good agreement with those observed experimentally. Besides well-mixed populations, set structured populations have also been taken into consideration. Therein we find that large number of sets and moderate migration rate greatly promote effort levels, especially for high effort costs.
Eco-evolutionary dynamics, coding structure and the information threshold
Hogeweg Paulien
2010-11-01
Full Text Available Abstract Background The amount of information that can be maintained in an evolutionary system of replicators is limited by genome length, the number of errors during replication (mutation rate and various external factors that influence the selection pressure. To date, this phenomenon, known as the information threshold, has been studied (both genotypically and phenotypically in a constant environment and with respect to maintenance (as opposed to accumulation of information. Here we take a broader perspective on this problem by studying the accumulation of information in an ecosystem, given an evolvable coding structure. Moreover, our setup allows for individual based as well as ecosystem based solutions. That is, all functions can be performed by individual replicators, or complementing functions can be performed by different replicators. In this setup, where both the ecosystem and the individual genomes can evolve their structure, we study how populations cope with high mutation rates and accordingly how the information threshold might be alleviated. Results We observe that the first response to increased mutation rates is a change in coding structure. At moderate mutation rates evolution leads to longer genomes with a higher diversity than at high mutation rates. Thus, counter-intuitively, at higher mutation rates diversity is reduced and the efficacy of the evolutionary process is decreased. Therefore, moderate mutation rates allow for more degrees of freedom in exploring genotype space during the evolutionary trajectory, facilitating the emergence of solutions. When an individual based solution cannot be attained due to high mutation rates, spatial structuring of the ecosystem can accommodate the evolution of ecosystem based solutions. Conclusions We conclude that the evolutionary freedom (eg. the number of genotypes that can be reached by evolution is increasingly restricted by higher mutation rates. In the case of such severe mutation
Geomedium as a nonlinear dynamic system. An evolutionary concept of earthquake development
Makarov, Pavel V., E-mail: pvm@ispms.tsc.ru [National Research Tomsk State University, Tomsk, 634050, Russia and Institute of Strength Physics and Materials Science SB RAS, Tomsk, 634055 (Russian Federation)
2014-11-14
An evolutionary approach to earthquake development is proposed. A medium under loading is treated as a multiscale nonlinear dynamic system. Its failure involves a number of stages typical of any dynamic system: dynamic chaos, self-organized criticality, and global stability loss in the final stage of its evolution. In the latter stage, the system evolves in a blow-up mode accompanied by catastrophic superfast movements of the elements of this geomedium.
Shahamatnia Ehsan
2016-01-01
Full Text Available Developing specialized software tools is essential to support studies of solar activity evolution. With new space missions such as Solar Dynamics Observatory (SDO, solar images are being produced in unprecedented volumes. To capitalize on that huge data availability, the scientific community needs a new generation of software tools for automatic and efficient data processing. In this paper a prototype of a modular framework for solar feature detection, characterization, and tracking is presented. To develop an efficient system capable of automatic solar feature tracking and measuring, a hybrid approach combining specialized image processing, evolutionary optimization, and soft computing algorithms is being followed. The specialized hybrid algorithm for tracking solar features allows automatic feature tracking while gathering characterization details about the tracked features. The hybrid algorithm takes advantages of the snake model, a specialized image processing algorithm widely used in applications such as boundary delineation, image segmentation, and object tracking. Further, it exploits the flexibility and efficiency of Particle Swarm Optimization (PSO, a stochastic population based optimization algorithm. PSO has been used successfully in a wide range of applications including combinatorial optimization, control, clustering, robotics, scheduling, and image processing and video analysis applications. The proposed tool, denoted PSO-Snake model, was already successfully tested in other works for tracking sunspots and coronal bright points. In this work, we discuss the application of the PSO-Snake algorithm for calculating the sidereal rotational angular velocity of the solar corona. To validate the results we compare them with published manual results performed by an expert.
Shahamatnia, Ehsan; Dorotovič, Ivan; Fonseca, Jose M.; Ribeiro, Rita A.
2016-03-01
Developing specialized software tools is essential to support studies of solar activity evolution. With new space missions such as Solar Dynamics Observatory (SDO), solar images are being produced in unprecedented volumes. To capitalize on that huge data availability, the scientific community needs a new generation of software tools for automatic and efficient data processing. In this paper a prototype of a modular framework for solar feature detection, characterization, and tracking is presented. To develop an efficient system capable of automatic solar feature tracking and measuring, a hybrid approach combining specialized image processing, evolutionary optimization, and soft computing algorithms is being followed. The specialized hybrid algorithm for tracking solar features allows automatic feature tracking while gathering characterization details about the tracked features. The hybrid algorithm takes advantages of the snake model, a specialized image processing algorithm widely used in applications such as boundary delineation, image segmentation, and object tracking. Further, it exploits the flexibility and efficiency of Particle Swarm Optimization (PSO), a stochastic population based optimization algorithm. PSO has been used successfully in a wide range of applications including combinatorial optimization, control, clustering, robotics, scheduling, and image processing and video analysis applications. The proposed tool, denoted PSO-Snake model, was already successfully tested in other works for tracking sunspots and coronal bright points. In this work, we discuss the application of the PSO-Snake algorithm for calculating the sidereal rotational angular velocity of the solar corona. To validate the results we compare them with published manual results performed by an expert.
The Ecology and Evolutionary Dynamics of Meiotic Drive.
Lindholm, Anna K; Dyer, Kelly A; Firman, Renée C; Fishman, Lila; Forstmeier, Wolfgang; Holman, Luke; Johannesson, Hanna; Knief, Ulrich; Kokko, Hanna; Larracuente, Amanda M; Manser, Andri; Montchamp-Moreau, Catherine; Petrosyan, Varos G; Pomiankowski, Andrew; Presgraves, Daven C; Safronova, Larisa D; Sutter, Andreas; Unckless, Robert L; Verspoor, Rudi L; Wedell, Nina; Wilkinson, Gerald S; Price, Tom A R
2016-04-01
Meiotic drivers are genetic variants that selfishly manipulate the production of gametes to increase their own rate of transmission, often to the detriment of the rest of the genome and the individual that carries them. This genomic conflict potentially occurs whenever a diploid organism produces a haploid stage, and can have profound evolutionary impacts on gametogenesis, fertility, individual behaviour, mating system, population survival, and reproductive isolation. Multiple research teams are developing artificial drive systems for pest control, utilising the transmission advantage of drive to alter or exterminate target species. Here, we review current knowledge of how natural drive systems function, how drivers spread through natural populations, and the factors that limit their invasion. PMID:26920473
Alvaro Sanchez
Full Text Available The evolutionary spread of cheater strategies can destabilize populations engaging in social cooperative behaviors, thus demonstrating that evolutionary changes can have profound implications for population dynamics. At the same time, the relative fitness of cooperative traits often depends upon population density, thus leading to the potential for bi-directional coupling between population density and the evolution of a cooperative trait. Despite the potential importance of these eco-evolutionary feedback loops in social species, they have not yet been demonstrated experimentally and their ecological implications are poorly understood. Here, we demonstrate the presence of a strong feedback loop between population dynamics and the evolutionary dynamics of a social microbial gene, SUC2, in laboratory yeast populations whose cooperative growth is mediated by the SUC2 gene. We directly visualize eco-evolutionary trajectories of hundreds of populations over 50-100 generations, allowing us to characterize the phase space describing the interplay of evolution and ecology in this system. Small populations collapse despite continual evolution towards increased cooperative allele frequencies; large populations with a sufficient number of cooperators "spiral" to a stable state of coexistence between cooperator and cheater strategies. The presence of cheaters does not significantly affect the equilibrium population density, but it does reduce the resilience of the population as well as its ability to adapt to a rapidly deteriorating environment. Our results demonstrate the potential ecological importance of coupling between evolutionary dynamics and the population dynamics of cooperatively growing organisms, particularly in microbes. Our study suggests that this interaction may need to be considered in order to explain intraspecific variability in cooperative behaviors, and also that this feedback between evolution and ecology can critically affect the
Evolutionary Fates and Dynamic Functionalization of Young Duplicate Genes in Arabidopsis Genomes.
Wang, Jun; Tao, Feng; Marowsky, Nicholas C; Fan, Chuanzhu
2016-09-01
Gene duplication is a primary means to generate genomic novelties, playing an essential role in speciation and adaptation. Particularly in plants, a high abundance of duplicate genes has been maintained for significantly long periods of evolutionary time. To address the manner in which young duplicate genes were derived primarily from small-scale gene duplication and preserved in plant genomes and to determine the underlying driving mechanisms, we generated transcriptomes to produce the expression profiles of five tissues in Arabidopsis thaliana and the closely related species Arabidopsis lyrata and Capsella rubella Based on the quantitative analysis metrics, we investigated the evolutionary processes of young duplicate genes in Arabidopsis. We determined that conservation, neofunctionalization, and specialization are three main evolutionary processes for Arabidopsis young duplicate genes. We explicitly demonstrated the dynamic functionalization of duplicate genes along the evolutionary time scale. Upon origination, duplicates tend to maintain their ancestral functions; but as they survive longer, they might be likely to develop distinct and novel functions. The temporal evolutionary processes and functionalization of plant duplicate genes are associated with their ancestral functions, dynamic DNA methylation levels, and histone modification abundances. Furthermore, duplicate genes tend to be initially expressed in pollen and then to gain more interaction partners over time. Altogether, our study provides novel insights into the dynamic retention processes of young duplicate genes in plant genomes. PMID:27485883
Wang, Baokui; Wang, Long
2014-01-01
The world in which we are living is a huge network of networks and should be described by interdependent networks. The interdependence between networks significantly affects the evolutionary dynamics of cooperation on them. Meanwhile, due to the diversity and complexity of social and biological systems, players on different networks may not interact with each other by the same way, which should be described by multiple models in evolutionary game theory, such as the Prisoner's Dilemma and Snowdrift Game. We therefore study the evolutionary dynamics of cooperation on two interdependent networks playing different games respectively. We clearly evidence that, with the increment of network interdependence, the evolution of cooperation is dramatically promoted on the network playing Prisoner's Dilemma. The cooperation level of the network playing Snowdrift Game reduces correspondingly, although it is almost invisible. In particular, there exists an optimal intermediate region of network interdependence maximizing ...
System dynamics of behaviour-evolutionary mix-game models
Gou Cheng-Ling; Gao Jie-Ping; Chen Fang
2010-01-01
In real financial markets there are two kinds of traders:one is fundamentalist,and the other is a trend-follower.The mix-game model is proposed to mimic such phenomena.In a mix-game model there are two groups of agents:Group 1 plays the majority game and Group 2 plays the minority game.In this paper,we investigate such a case that some traders in real financial markets could change their investment behaviours by assigning the evolutionary abilities to agents:if the winning rates of agents are smaller than a threshold,they will join the other group;and agents will repeat such an evolution at certain time intervals.Through the simulations,we obtain the following findings:(i) the volatilities of systems increase with the increase of the number of agents in Group 1 and the times of behavioural changes of all agents;(ii) the performances of agents in both groups and the stabilities of systems become better if all agents take more time to observe their new investment behaviours;(iii) there are two-phase zones of market and non-market and two-phase zones of evolution and non-evolution;(iv) parameter configurations located within the cross areas between the zones of markets and the zones of evolution are suited for simulating the financial markets.
System dynamics of behaviour-evolutionary mix-game models
Gou, Cheng-Ling; Gao, Jie-Ping; Chen, Fang
2010-11-01
In real financial markets there are two kinds of traders: one is fundamentalist, and the other is a trend-follower. The mix-game model is proposed to mimic such phenomena. In a mix-game model there are two groups of agents: Group 1 plays the majority game and Group 2 plays the minority game. In this paper, we investigate such a case that some traders in real financial markets could change their investment behaviours by assigning the evolutionary abilities to agents: if the winning rates of agents are smaller than a threshold, they will join the other group; and agents will repeat such an evolution at certain time intervals. Through the simulations, we obtain the following findings: (i) the volatilities of systems increase with the increase of the number of agents in Group 1 and the times of behavioural changes of all agents; (ii) the performances of agents in both groups and the stabilities of systems become better if all agents take more time to observe their new investment behaviours; (iii) there are two-phase zones of market and non-market and two-phase zones of evolution and non-evolution; (iv) parameter configurations located within the cross areas between the zones of markets and the zones of evolution are suited for simulating the financial markets.
Non-Payoff Monotonic Dynamics in an Evolutionary Game of Courtship
Chacoma, Andrés; Zanette, Damián H
2015-01-01
We propose an evolutionary coordination game to formalize a simplified model of the evolution of strategies during human courtship. The dynamics, derived from the consideration of experimental observations on human social behavior driven by self-esteem, turns out to be non-payoff monotonic. This property gives rise to nontrivial evolution in the players' strategies, which we study both numerically and analytically.
A stochastic evolutionary model for capturing human dynamics
Fenner, Trevor; Loizou, George
2015-01-01
The recent interest in human dynamics has led researchers to investigate the stochastic processes that explain human behaviour in various contexts. Here we propose a generative model to capture the dynamics of survival analysis, traditionally employed in clinical trials and reliability analysis in engineering. We derive a general solution for the model in the form of a product, and then a continuous approximation to the solution via the renewal equation describing age-structured population dynamics. This enables us to model a wide rage of survival distributions, according to the choice of the mortality distribution. We provide empirical evidence for the validity of the model from a longitudinal data set of popular search engine queries over 114 months, showing that the survival function of these queries is closely matched by the solution for our model with power-law mortality.
Dynamic evolutionary community detection algorithms based on the modularity matrix
Motivated by the relationship of the dynamic behaviors and network structure, in this paper, we present two efficient dynamic community detection algorithms. The phases of the nodes in the network can evolve according to our proposed differential equations. In each iteration, the phases of the nodes are controlled by several parameters. It is found that the phases of the nodes are ultimately clustered into several communities after a short period of evolution. They can be adopted to detect the communities successfully. The second differential equation can dynamically adjust several parameters, so it can obtain satisfactory detection results. Simulations on some test networks have verified the efficiency of the presented algorithms. (interdisciplinary physics and related areas of science and technology)
Dynamic Waste Management (DWM):Towards an evolutionary decision-making approach
Rojo, Gabriel; Glaus, Mathias; Hausler, Robert; Laforest, Valérie; Bourgois, Jacques
2013-01-01
To guarantee sustainable and dynamic waste management, the dynamic waste management approach (DWM) suggests an evolutionary new approach that maintains a constant flow towards the most favourable waste treatment processes (facilities) within a system. To that end, DWM is based on the law of conservation of energy, which allows the balancing of a network, while considering the constraints of incoming (h1) and outgoing (h2) loads, as well as the distribution network (ΔH) characteristics. The de...
The Tangled Nature Model of evolutionary dynamics reconsidered
Andersen, Christian Walther; Sibani, Paolo
2016-01-01
elephant. To bring out the structural and dynamical effects of trait inheritance , we introduce and numerically analyze a family of TNM models where a positive integer $K$ parametrises correlations between the interactions of an agent and those of its mutated offspring. For $K=1$ a single point mutation...
Achieving Reliability in Master-worker Computing via Evolutionary Dynamics
Christoforou, Evgenia; Fernández Anta, Antonio; Georgiou, Chryssis; Mosteiro, Miguel A.; Sánchez, Angel
2012-01-01
This work considers Internet-based task computations in which a master process assigns tasks, over the Internet, to rational workers and collect their responses. The objective is for the master to obtain the correct task outcomes. For this purpose we formulate and study the dynamics of evolution of Internet-based master-worker computations through reinforcement learning.
Bandy, D. K.; Hall, J. R.; Denker, M. E.
2015-07-01
We show that the role of the Lyapunov exponents can be extended beyond the customary local instability, such as limit cycle behavior, to include its use as an evolutionary predictor of the dynamics of a laser with injected signal (LIS). Numerical studies of LIS reveal that as a function of the input-signal strength the evolution of two nonzero Lyapunov exponents (generally equal) distinctively predicts the evolutionary trend of the fundamental frequency of the laser output signal (an important dynamic characteristic of the LIS) even with the presence of some noise. This globally predictive behavior of the Lyapunov exponents includes also the dynamic behavior of the individual coexisting attractors. Different coexisting attractors of LIS and configurations of Lyapunov exponents for both individual attractors and the global system are reported. Two LIS case studies are considered: (I) a high-gain system with a rich history of nonlinear behavior but not experimentally accessible, and (II) a low-gain system that has complex dynamics and is experimentally accessible for Class B lasers. Universality arguments support the thesis that these different configurations and the extended role of the Lyapunov exponents as an evolutionary predictor of the dynamics will be observed in other nonlinear, dynamic dissipative systems as well.
Pesole Graziano
2009-09-01
Full Text Available Abstract Background The conservation of sequences between related genomes has long been recognised as an indication of functional significance and recognition of sequence homology is one of the principal approaches used in the annotation of newly sequenced genomes. In the context of recent findings that the number non-coding transcripts in higher organisms is likely to be much higher than previously imagined, discrimination between conserved coding and non-coding sequences is a topic of considerable interest. Additionally, it should be considered desirable to discriminate between coding and non-coding conserved sequences without recourse to the use of sequence similarity searches of protein databases as such approaches exclude the identification of novel conserved proteins without characterized homologs and may be influenced by the presence in databases of sequences which are erroneously annotated as coding. Results Here we present a machine learning-based approach for the discrimination of conserved coding sequences. Our method calculates various statistics related to the evolutionary dynamics of two aligned sequences. These features are considered by a Support Vector Machine which designates the alignment coding or non-coding with an associated probability score. Conclusion We show that our approach is both sensitive and accurate with respect to comparable methods and illustrate several situations in which it may be applied, including the identification of conserved coding regions in genome sequences and the discrimination of coding from non-coding cDNA sequences.
Temporal and Evolutionary Dynamics of Two-Component Signaling Pathways
Salazar, Michael E.; Laub, Michael T.
2015-01-01
Bacteria sense and respond to numerous environmental signals through two-component signaling pathways. Typically, a given stimulus will activate a sensor histidine kinase to autophosphorylate and then phosphotransfer to a cognate response regulator, which can mount an appropriate response. Although these signaling pathways often appear to be simple switches, they can also orchestrate surprisingly sophisticated and complex responses. These temporal dynamics arise from several key regulatory fe...
Modeling evolutionary dynamics of epigenetic mutations in hierarchically organized tumors.
Andrea Sottoriva; Louis Vermeulen; Simon Tavaré
2011-01-01
The cancer stem cell (CSC) concept is a highly debated topic in cancer research. While experimental evidence in favor of the cancer stem cell theory is apparently abundant, the results are often criticized as being difficult to interpret. An important reason for this is that most experimental data that support this model rely on transplantation studies. In this study we use a novel cellular Potts model to elucidate the dynamics of established malignancies that are driven by a small subset of ...
Modeling Evolutionary Dynamics of Epigenetic Mutations in Hierarchically Organized Tumors
Sottoriva, A.; Vermeulen, L; Tavare, S
2011-01-01
The cancer stem cell (CSC) concept is a highly debated topic in cancer research. While experimental evidence in favor of the cancer stem cell theory is apparently abundant, the results are often criticized as being difficult to interpret. An important reason for this is that most experimental data that support this model rely on transplantation studies. In this study we use a novel cellular Potts model to elucidate the dynamics of established malignancies that are driven by a small subset of ...
Eco-evolutionary feedback promotes Red Queen dynamics and selects for sex in predator populations.
Haafke, Julia; Abou Chakra, Maria; Becks, Lutz
2016-03-01
Although numerous hypotheses exist to explain the overwhelming presence of sexual reproduction across the tree of life, we still cannot explain its prevalence when considering all inherent costs involved. The Red Queen hypothesis states that sex is maintained because it can create novel genotypes with a selective advantage. This occurs when the interactions between species induce frequent environmental change. Here, we investigate whether coevolution and eco-evolutionary feedback dynamics in a predator-prey system allows for indirect selection and maintenance of sexual reproduction in the predator. Combining models and chemostat experiments of a rotifer-algae system we show a continuous feedback between population and trait change along with recurrent shifts from selection by predation and competition for a limited resource. We found that a high propensity for sex was indirectly selected and was maintained in rotifer populations within environments containing these eco-evolutionary dynamics; whereas within environments under constant conditions, predators evolved rapidly to lower levels of sex. Thus, our results indicate that the influence of eco-evolutionary feedback dynamics on the overall evolutionary change has been underestimated. PMID:26899793
Calculation of dynamics of silver sorption by zirconium phosphate
The experimental and calculation data on kinetics and dynamics of silver ions sorption by zirconium phosphate are considered in present article. The evaluation of sorption process in dynamic conditions is carried out.
Dynamic calculations of pressurized water reactor internals
A mathematical model is briefly described for the calculation of oscillations in the WWER-440 reactor internals. The model was developed for improved safety of the type of reactors. It allows calculating vibrations resistance of reactor components, mainly during accidents, such as loss of coolant accidents. Some results are given of the calculation of forces acting in the rupture of the reactor inlet and outlet pipes. (Z.M.)
How the Magnitude of Prey Genetic Variation Alters Predator-Prey Eco-Evolutionary Dynamics.
Cortez, Michael H
2016-09-01
Evolution can alter the stability and dynamics of ecological communities; for example, prey evolution can drive cyclic dynamics in predator-prey systems that are not possible in the absence of evolution. However, it is unclear how the magnitude of additive genetic variation in the evolving species mediates those effects. In this study, I explore how the magnitude of prey additive genetic variation determines what effects prey evolution has on the dynamics and stability of predator-prey systems. I use linear stability analysis to decompose the stability of a general eco-evolutionary predator-prey model into components representing the stabilities of the ecological and evolutionary subsystems as well as the interactions between those subsystems. My results show that with low genetic variation, the cyclic dynamics and stability of the system are determined by the ecological subsystem. With increased genetic variation, disruptive selection always destabilizes stable communities, stabilizing selection can stabilize or destabilize communities, and prey evolution can alter predator-prey phase lags. Stability changes occur approximately when the magnitude of genetic variation balances the (in)stabilities of the ecological and evolutionary subsystems. I discuss the connections between my stability results and prior results from the theory of adaptive dynamics. PMID:27501090
Friman, Ville-Petri; Dupont, Alessandra; Bass, David; Murrell, David J; Bell, Thomas
2016-06-01
Community dynamics are often studied in subsets of pairwise interactions. Scaling pairwise interactions back to the community level is, however, problematic because one given interaction might not reflect ecological and evolutionary outcomes of other functionally similar species interactions or capture the emergent eco-evolutionary dynamics arising only in more complex communities. Here we studied this experimentally by exposing Pseudomonas fluorescens SBW25 prey bacterium to four different protist predators (Tetrahymena pyriformis, Tetrahymena vorax, Chilomonas paramecium and Acanthamoeba polyphaga) in all possible single-predator, two-predator and four-predator communities for hundreds of prey generations covering both ecological and evolutionary timescales. We found that only T. pyriformis selected for prey defence in single-predator communities. Although T. pyriformis selection was constrained in the presence of the intraguild predator, T. vorax, T. pyriformis selection led to evolution of specialised prey defence strategies in the presence of C. paramecium or A. polyphaga. At the ecological level, adapted prey populations were phenotypically more diverse, less stable and less productive compared with non-adapted prey populations. These results suggest that predator community composition affects the relative importance of ecological and evolutionary processes and can crucially determine when rapid evolution has the potential to change ecological properties of microbial communities. PMID:26684728
Evolutionary dynamics on stochastic evolving networks for multiple-strategy games
Wu, Bin; Zhou, Da; Wang, Long
2011-10-01
Evolutionary game theory on dynamical networks has received much attention. Most of the work has been focused on 2×2 games such as prisoner's dilemma and snowdrift, with general n×n games seldom addressed. In particular, analytical methods are still lacking. Here we generalize the stochastic linking dynamics proposed by Wu, Zhou, Fu, Luo, Wang, and Traulsen [PLoS ONEBSYMBO1932-620310.1371/journal.pone.0011187 5, e11187 (2010)] to n×n games. We analytically obtain that the fast linking dynamics results in the replicator dynamics with a rescaled payoff matrix. In the rescaled matrix, intuitively, each entry is the product of the original entry and the average duration time of the corresponding link. This result is shown to be robust to a wide class of imitation processes. As applications, we show both analytically and numerically that the biodiversity, modeled as the stability of a zero-sum rock-paper-scissors game, cannot be altered by the fast linking dynamics. In addition, we show that the fast linking dynamics can stabilize tit-for-tat as an evolutionary stable strategy in the repeated prisoner's dilemma game provided the interaction between the identical strategies happens sufficiently often. Our method paves the way for an analytical study of the multiple-strategy coevolutionary dynamics.
Evolutionary calculations of carbon dredge-up in helium envelope white dwarfs
Macdonald, James; Hernanz, Margarita; José, Jordi
1998-01-01
We investigate the evolution of cooling helium atmosphere white dwarfs using a full evolutionary code, specifically developed for following the effects of element diffusion and gravitational settling on white dwarf cooling. The major difference between this work and previous work is that we use more recent opacity data from the OPAL project. Since, in general, these opacities are higher than those available ten years ago, at a given effective temperature, convection zones go deeper than in mo...
Epidemiological and Evolutionary Dynamics of Influenza B Viruses in Malaysia, 2012-2014
Xiang Yong Oong; Kim Tien Ng; Tommy Tsan-Yuk Lam; Yong Kek Pang; Kok Gan Chan; Nik Sherina Hanafi; Adeeba Kamarulzaman; Kok Keng Tee
2015-01-01
Epidemiological and evolutionary dynamics of influenza B Victoria and Yamagata lineages remained poorly understood in the tropical Southeast Asia region, despite causing seasonal outbreaks worldwide. From 2012-2014, nasopharyngeal swab samples collected from outpatients experiencing acute upper respiratory tract infection symptoms in Kuala Lumpur, Malaysia, were screened for influenza viruses using a multiplex RT-PCR assay. Among 2,010/3,935 (51.1%) patients infected with at least one respira...
Mutation-selection dynamics and error threshold in an evolutionary model for Turing Machines
Musso, Fabio; Feverati, Giovanni
2011-01-01
We investigate the mutation-selection dynamics for an evolutionary computation model based on Turing Machines that we introduced in a previous article. The use of Turing Machines allows for very simple mechanisms of code growth and code activation/inactivation through point mutations. To any value of the point mutation probability corresponds a maximum amount of active code that can be maintained by selection and the Turing machines that reach it are said to be at the error threshold. Simulat...
The Role of Clonal Interference in the Evolutionary Dynamics of Plasmid-Host Adaptation
Hughes, Julie M; Lohman, Brian K.; Deckert, Gail E.; Nichols, Eric P.; Settles, Matt; Abdo, Zaid; Eva M. Top
2012-01-01
ABSTRACT Promiscuous plasmids replicate in a wide range of bacteria and therefore play a key role in the dissemination of various host-beneficial traits, including antibiotic resistance. Despite the medical relevance, little is known about the evolutionary dynamics through which drug resistance plasmids adapt to new hosts and thereby persist in the absence of antibiotics. We previously showed that the incompatibility group P-1 (IncP-1) minireplicon pMS0506 drastically improved its stability i...
Evolutionary supervision of a dynamical neural network allows learning with on-going weights
Meunier, David; Paugam-Moisy, Hélène
2005-01-01
Recent electrophysiological data show that synaptic weights are highly influenced by electrical activities displayed by neurons. Weights are not stable as assumed in classical neural network models. What is the nature of engrams, if not stored in synaptic weights? Adopting the theory of dynamical systems, which allows an implicit form of memory, we propose a new framework for learning, where synaptic weights are continuously adapted. Evolutionary computation has been applied to a population o...
Senerchia, Natacha; Wicker, Thomas; Felber, François; Parisod, Christian
2013-01-01
Transposable elements (TEs) represent a major fraction of plant genomes and drive their evolution. An improved understanding of genome evolution requires the dynamics of a large number of TE families to be considered. We put forward an approach bypassing the required step of a complete reference genome to assess the evolutionary trajectories of high copy number TE families from genome snapshot with high-throughput sequencing. Low coverage sequencing of the complex genomes of Aegilops cylindri...
Bando, Kathy J.
2005-01-01
Dwarf eelgrass (Zostera japonica) and smooth cordgrass (Spartina alterniflora) are ecologically important invaders of intertidal mudflats in the eastern Pacific. S. alterniflora and Z. japonica invasions alter estuarine nutrient dynamics, cause sediment and infaunal community changes, and modify waterfowl foraging habitats. Ecological and evolutionary mechanisms of the invasion success of Z. japonica and S. alterniflora were addressed with a combination of experimental and observational field...
Macroecology meets macroevolution: evolutionary niche dynamics in the seaweed Halimeda
Verbruggen, H.; Tyberghein, L.; Pauly, K.; C. Vlaeminck; VAN NIEUWENHUYZE, K; Kooistra, W.H.C.F.; Leliaert, F.; De Clerck, O
2011-01-01
Aim Because of their broad distribution in geographical and ecological dimensions, seaweeds (marine macroalgae) offer great potential as models for marine biogeographical inquiry and exploration of the interface between macroecology and macroevolution. This study aims to characterize evolutionary niche dynamics in the common green seaweed genus Halimeda, use the observed insights to gain understanding of the biogeographical history of the genus and predict habitats that can be targeted for t...
Fournier-Level, Alexandre; Perry, Emily O; Wang, Jonathan A; Braun, Peter T; Migneault, Andrew; Cooper, Martha D; Metcalf, C Jessica E; Schmitt, Johanna
2016-05-17
Predicting whether and how populations will adapt to rapid climate change is a critical goal for evolutionary biology. To examine the genetic basis of fitness and predict adaptive evolution in novel climates with seasonal variation, we grew a diverse panel of the annual plant Arabidopsis thaliana (multiparent advanced generation intercross lines) in controlled conditions simulating four climates: a present-day reference climate, an increased-temperature climate, a winter-warming only climate, and a poleward-migration climate with increased photoperiod amplitude. In each climate, four successive seasonal cohorts experienced dynamic daily temperature and photoperiod variation over a year. We measured 12 traits and developed a genomic prediction model for fitness evolution in each seasonal environment. This model was used to simulate evolutionary trajectories of the base population over 50 y in each climate, as well as 100-y scenarios of gradual climate change following adaptation to a reference climate. Patterns of plastic and evolutionary fitness response varied across seasons and climates. The increased-temperature climate promoted genetic divergence of subpopulations across seasons, whereas in the winter-warming and poleward-migration climates, seasonal genetic differentiation was reduced. In silico "resurrection experiments" showed limited evolutionary rescue compared with the plastic response of fitness to seasonal climate change. The genetic basis of adaptation and, consequently, the dynamics of evolutionary change differed qualitatively among scenarios. Populations with fewer founding genotypes and populations with genetic diversity reduced by prior selection adapted less well to novel conditions, demonstrating that adaptation to rapid climate change requires the maintenance of sufficient standing variation. PMID:27140640
Evolutionary dynamics of imatinib-treated leukemic cells by stochastic approach
Pizzolato, Nicola; Valenti, Davide; Adorno, Dominique; Spagnolo, Bernardo
2009-09-01
The evolutionary dynamics of a system of cancerous cells in a model of chronic myeloid leukemia (CML) is investigated by a statistical approach. Cancer progression is explored by applying a Monte Carlo method to simulate the stochastic behavior of cell reproduction and death in a population of blood cells which can experience genetic mutations. In CML front line therapy is represented by the tyrosine kinase inhibitor imatinib which strongly affects the reproduction of leukemic cells only. In this work, we analyze the effects of a targeted therapy on the evolutionary dynamics of normal, first-mutant and cancerous cell populations. Several scenarios of the evolutionary dynamics of imatinib-treated leukemic cells are described as a consequence of the efficacy of the different modelled therapies. We show how the patient response to the therapy changes when a high value of the mutation rate from healthy to cancerous cells is present. Our results are in agreement with clinical observations. Unfortunately, development of resistance to imatinib is observed in a fraction of patients, whose blood cells are characterized by an increasing number of genetic alterations. We find that the occurrence of resistance to the therapy can be related to a progressive increase of deleterious mutations.
The amazing evolutionary dynamics of non-linear optical systems with feedback
Yaroslavsky, Leonid
2013-09-01
Optical systems with feedback are, generally, non-linear dynamic systems. As such, they exhibit evolutionary behavior. In the paper we present results of experimental investigation of evolutionary dynamics of several models of such systems. The models are modifications of the famous mathematical "Game of Life". The modifications are two-fold: "Game of Life" rules are made stochastic and mutual influence of cells is made spatially non-uniform. A number of new phenomena in the evolutionary dynamics of the models are revealed: - "Ordering of chaos". Formation, from seed patterns, of stable maze-like patterns with chaotic "dislocations" that resemble natural patterns, such as skin patterns of some animals and fishes, see shell, fingerprints, magnetic domain patterns and alike, which one can frequently find in the nature. These patterns and their fragments exhibit a remarkable capability of unlimited growth. - "Self-controlled growth" of chaotic "live" formations into "communities" bounded, depending on the model, by a square, hexagon or octagon, until they reach a certain critical size, after which the growth stops. - "Eternal life in a bounded space" of "communities" after reaching a certain size and shape. - "Coherent shrinkage" of "mature", after reaching a certain size, "communities" into one of stable or oscillating patterns preserving in this process isomorphism of their bounding shapes until the very end.
Ito, Takashi; Pilat, Marcin L; Suzuki, Reiji; Arita, Takaya
2016-01-01
Recent studies have reported that population dynamics and evolutionary dynamics, occurring at different time scales, can be affected by each other. Our purpose is to explore the interaction between population and evolutionary dynamics using an artificial life approach based on a 3D physically simulated environment in the context of predator-prey and morphology-behavior coevolution. The morphologies and behaviors of virtual prey creatures are evolved using a genetic algorithm based on the predation interactions between predators and prey. Both population sizes are also changed, depending on the fitness. We observe two types of cyclic behaviors, corresponding to short-term and long-term dynamics. The former can be interpreted as a simple population dynamics of Lotka-Volterra type. It is shown that the latter cycle is based on the interaction between the changes in the prey strategy against predators and the long-term change in both population sizes, resulting partly from a tradeoff between their defensive success and the cost of defense. PMID:26934093
Dynamical collective calculation of supernova neutrino signals.
Gava, Jérôme; Kneller, James; Volpe, Cristina; McLaughlin, G C
2009-08-14
We present the first calculations with three flavors of collective and shock wave effects for neutrino propagation in core-collapse supernovae using hydrodynamical density profiles and the S matrix formalism. We explore the interplay between the neutrino-neutrino interaction and the effects of multiple resonances upon the time signal of positrons in supernova observatories. A specific signature is found for the inverted hierarchy and a large third neutrino mixing angle and we predict, in this case, a dearth of lower energy positrons in Cherenkov detectors midway through the neutrino signal and the simultaneous revelation of valuable information about the original fluxes. We show that this feature is also observable with current generation neutrino detectors at the level of several sigmas. PMID:19792628
Analysis of Ant Colony Optimization and Population-Based Evolutionary Algorithms on Dynamic Problems
Lissovoi, Andrei
exist more complex oscillations that cannot be tracked with a polynomial-size colony. MMAS and (μ+1) EA on Maze We analyse the behaviour of a (μ + 1) EA with genotype diversity on a dynamic fitness function Maze, extended to a finite-alphabet search space. We prove that the (μ + 1) EA is able to track...... the dynamic optimum for finite alphabets up to size μ, while MMAS is able to do so for any finite alphabet size. Parallel Evolutionary Algorithms on Maze. We prove that while a (1 + λ) EA is unable to track the optimum of the dynamic fitness function Maze for offspring population size up to λ = O(n1-ε...... analysis showing how closely the EA can track the dynamically moving optimum over time. These results are also extended to a finite-alphabet search space....
Evolutionary Design of Both Topologies and Parameters of a Hybrid Dynamical System
Dupuis, Jean-Francois; Fan, Zhun; Goodman, Erik
2012-01-01
This paper investigates the issue of evolutionary design of open-ended plants for hybrid dynamical systems--i.e. both their topologies and parameters. Hybrid bond graphs are used to represent dynamical systems involving both continuous and discrete system dynamics. Genetic programming, with some...... special mechanisms incorporated, is used as a search tool to explore the open-ended design space of hybrid bond graphs. Combination of these two tools--i.e., hybrid bond graphs (HGBs) and genetic programming (GP)--leads to an approach called HBGGP that can automatically generate viable design candidates...... of hybrid dynamical systems that fulfill predefined design specifications. A comprehensive investigation of a case study of DC-DC converter design demonstrates the feasibility and effectiveness of the HBGGP approach. Important characteristics of the approach are also discussed, with some future...
Dynamics of evolutionary rescue in changing environments and the emergence of antibiotic resistance.
Wu, Yue; Saddler, Clare A; Valckenborgh, Frank; Tanaka, Mark M
2014-01-01
Populations can go extinct when their environments deteriorate, but evolutionary rescue occurs when a shrinking population adapts to the new environmental conditions. The emergence of resistance from a drug sensitive bacterial population under treatment can be regarded as an instance of evolutionary rescue. Understanding evolutionary rescue in a particular context such as drug resistance requires knowledge of how the environment changes and how selection coefficients change as a result. In this study, we propose a model for evolutionary rescue under three different scenarios of environmental change: abrupt change, periodic fluctuation and gradual decay. The model makes use of the notion of reaction norms to describe fitness values that depend on both genotype and environmental state. We find that although drug sensitive bacterial populations may be large, allowing them to generate resistant mutants frequently, a harsh abrupt change due to the drug usually drives them extinct. Evolutionary rescue occurs far more frequently under the milder forms of environmental change we investigated. Rescue is favoured when the absolute fitnesses of individuals remain sufficiently high over the range of environment qualities experienced by the population. The minimum environment quality, which is inversely related to drug dose in the antibiotic context, is thus an important factor. Interestingly, in the periodic fluctuation model, the inter-dose period is less influential in promoting rescue through resistance unless the minimum environment quality is in a particular range. We also investigated fitness trade-offs across environments including the case of a resistant allele not subject to any trade-off (a "superbug"). This fitness trade-off affects the probability of rescue in decaying environments, but surprisingly has only a weak effect in the periodic fluctuation scenario. Finally, we use the model to show how niche construction, whereby organisms are the source of environmental
Ancient origin of the tryptophan operon and the dynamics of evolutionary change.
Xie, Gary; Keyhani, Nemat O; Bonner, Carol A; Jensen, Roy A
2003-09-01
The seven conserved enzymatic domains required for tryptophan (Trp) biosynthesis are encoded in seven genetic regions that are organized differently (whole-pathway operons, multiple partial-pathway operons, and dispersed genes) in prokaryotes. A comparative bioinformatics evaluation of the conservation and organization of the genes of Trp biosynthesis in prokaryotic operons should serve as an excellent model for assessing the feasibility of predicting the evolutionary histories of genes and operons associated with other biochemical pathways. These comparisons should provide a better understanding of possible explanations for differences in operon organization in different organisms at a genomics level. These analyses may also permit identification of some of the prevailing forces that dictated specific gene rearrangements during the course of evolution. Operons concerned with Trp biosynthesis in prokaryotes have been in a dynamic state of flux. Analysis of closely related organisms among the Bacteria at various phylogenetic nodes reveals many examples of operon scission, gene dispersal, gene fusion, gene scrambling, and gene loss from which the direction of evolutionary events can be deduced. Two milestone evolutionary events have been mapped to the 16S rRNA tree of Bacteria, one splitting the operon in two, and the other rejoining it by gene fusion. The Archaea, though less resolved due to a lesser genome representation, appear to exhibit more gene scrambling than the Bacteria. The trp operon appears to have been an ancient innovation; it was already present in the common ancestor of Bacteria and Archaea. Although the operon has been subjected, even in recent times, to dynamic changes in gene rearrangement, the ancestral gene order can be deduced with confidence. The evolutionary history of the genes of the pathway is discernible in rough outline as a vertical line of descent, with events of lateral gene transfer or paralogy enriching the analysis as interesting
Beam dynamics calculations for fault-tolerance
The European Transmutation Demonstration requires a high-power proton accelerator operating in CW mode. This accelerator is also expected to have a very limited number of unexpected beam interruptions per year. To reach such an ambitious goal, it is clear that reliability-oriented design practices need to be followed from the early stage of components design and fault-tolerance capabilities have to be introduced to the maximum extent. The goal of this document is precisely to investigate in more details the fault-tolerance capability of the XT-ADS linac. From previous analysis, it appears that if nothing is done, a cavity's failure leads in nearly all the cases to a complete beam loss, due to the non-relativistic varying velocity of the particles. To avoid such a total beam loss, it is clear that some kind of retuning has to be performed to compensate the lack of acceleration due to the faulty cavity. We have to identify and develop fast failure recovery scenarios to ensure that such retuning can be performed in less than 1 second. 2 ways are investigated. The first way is to stop the beam to achieve the retuning (Scenario 1). The other way is to try to perform the retuning without stopping the beam (Scenario 2). The present analysis demonstrates on the beam dynamics point of view that a fast retuning procedure can be envisaged without stopping the beam (Scenario 2). Nevertheless, this Scenario 2 implies stringent specifications, especially on: - the fault detection time, that has to be extremely short (order of magnitude: 100 μs) and - the margins required on the accelerating field and RF power point of view, that are higher than in Scenario 1
On the calculation of dynamic derivatives using computational fluid dynamics
Da Ronch, Andrea
2012-01-01
In this thesis, the exploitation of computational fluid dynamics (CFD) methods for the flight dynamics of manoeuvring aircraft is investigated. It is demonstrated that CFD can now be used in a reasonably routine fashion to generate stability and control databases. Different strategies to create CFD-derived simulation models across the flight envelope are explored, ranging from combined low-fidelity/high-fidelity methods to reduced-order modelling. For the representation of the unsteady aerody...
Monte Carlo dose calculations for dynamic IMRT treatments
Dose calculations for intensity modulated radiation therapy (IMRT) face new challenges due to the complex leaf geometry and time dependent nature of the delivery. A fast method of particle transport through a dynamic multileaf collimator (MLC) geometry that accounts for photon attenuation and first-scattered Compton photon production has been incorporated into an existing Monte Carlo code used for patient dose calculations. Dosimetric agreement between calculation and measurement for two photon energies and MLC types is within experimental error for the sliding window tests. For a patient IMRT field, the Monte Carlo calculations are closer to measured dose than similar superposition or pencil beam calculations. (author)
This paper proposes a new hybrid approach based on nonlinear chaotic dynamics and evolutionary strategy to forecast electricity loads and prices. The main idea is to develop a new training or identification stage in a nonlinear chaotic dynamic based predictor. In the training stage five optimal parameters for a chaotic based predictor are searched through an optimization model based on evolutionary strategy. The objective function of the optimization model is the mismatch minimization between the multi-step-ahead forecasting of predictor and observed data such as it is done in identification problems. The first contribution of this paper is that the proposed approach is capable of capturing the complex dynamic of demand and price time series considered resulting in a more accuracy forecasting. The second contribution is that the proposed approach run on-line manner, i.e. the optimal set of parameters and prediction is executed automatically which can be used to prediction in real-time, it is an advantage in comparison with other models, where the choice of their input parameters are carried out off-line, following qualitative/experience-based recipes. A case study of load and price forecasting is presented using data from New England, Alberta, and Spain. A comparison with other methods such as autoregressive integrated moving average (ARIMA) and artificial neural network (ANN) is shown. The results show that the proposed approach provides a more accurate and effective forecasting than ARIMA and ANN methods. (author)
An experimental investigation of evolutionary dynamics in the Rock-Paper-Scissors game.
Hoffman, Moshe; Suetens, Sigrid; Gneezy, Uri; Nowak, Martin A
2015-01-01
Game theory describes social behaviors in humans and other biological organisms. By far, the most powerful tool available to game theorists is the concept of a Nash Equilibrium (NE), which is motivated by perfect rationality. NE specifies a strategy for everyone, such that no one would benefit by deviating unilaterally from his/her strategy. Another powerful tool available to game theorists are evolutionary dynamics (ED). Motivated by evolutionary and learning processes, ED specify changes in strategies over time in a population, such that more successful strategies typically become more frequent. A simple game that illustrates interesting ED is the generalized Rock-Paper-Scissors (RPS) game. The RPS game extends the children's game to situations where winning or losing can matter more or less relative to tying. Here we investigate experimentally three RPS games, where the NE is always to randomize with equal probability, but the evolutionary stability of this strategy changes. Consistent with the prediction of ED we find that aggregate behavior is far away from NE when it is evolutionarily unstable. Our findings add to the growing literature that demonstrates the predictive validity of ED in large-scale incentivized laboratory experiments with human subjects. PMID:25743257
Qian, Hong
2012-01-01
The study of biological cells in terms of mesoscopic, nonequilibrium, nonlinear, stochastic dynamics of open chemical systems provides a paradigm for other complex, self-organizing systems with ultra-fast stochastic fluctuations, short-time deterministic nonlinear dynamics, and long-time evolutionary behavior with exponentially distributed rare events, discrete jumps among punctuated equilibria, and catastrophe.
Evolutionary dynamics of public goods games with diverse contributions in finite populations
Wang, Jing; Wu, Bin; Chen, Xiaojie; Wang, Long
2010-05-01
The public goods game is a powerful metaphor for exploring the maintenance of social cooperative behavior in a group of interactional selfish players. Here we study the emergence of cooperation in the public goods games with diverse contributions in finite populations. The theory of stochastic process is innovatively adopted to investigate the evolutionary dynamics of the public goods games involving a diversity of contributions. In the limit of rare mutations, the general stationary distribution of this stochastic process can be analytically approximated by means of diffusion theory. Moreover, we demonstrate that increasing the diversity of contributions greatly reduces the probability of finding the population in a homogeneous state full of defectors. This increase also raises the expectation of the total contribution in the entire population and thus promotes social cooperation. Furthermore, by investigating the evolutionary dynamics of optional public goods games with diverse contributions, we find that nonparticipation can assist players who contribute more in resisting invasion and taking over individuals who contribute less. In addition, numerical simulations are performed to confirm our analytical results. Our results may provide insight into the effect of diverse contributions on cooperative behaviors in the real world.
The co-evolutionary dynamics of directed network of spin market agents
Horváth, Denis; Kuscsik, Zoltán; Gmitra, Martin
2006-09-01
The spin market model [S. Bornholdt, Int. J. Mod. Phys. C 12 (2001) 667] is generalized by employing co-evolutionary principles, where strategies of the interacting and competitive traders are represented by local and global couplings between the nodes of dynamic directed stochastic network. The co-evolutionary principles are applied in the frame of Bak-Sneppen self-organized dynamics [P. Bak, K. Sneppen, Phys. Rev. Lett. 71 (1993) 4083] that includes the processes of selection and extinction actuated by the local (node) fitness. The local fitness is related to orientation of spin agent with respect to the instant magnetization. The stationary regime is formed due to the interplay of self-organization and adaptivity effects. The fat tailed distributions of log-price returns are identified numerically. The non-trivial model consequence is the evidence of the long time market memory indicated by the power-law range of the autocorrelation function of volatility with exponent smaller than one. The simulations yield network topology with broad-scale node degree distribution characterized by the range of exponents 1.3<γin<3 coinciding with social networks.
Two-stage evolutionary algorithm for dynamic multicast routing in mesh network
Li ZHU; Zhi-shu LI; Liang-yin CHEN; Yan-hong CHENG
2008-01-01
In order to share multimedia transmissions in mesh networks and optimize the utilization of network resources, this paper presents a Two-stage Evolutionary Algorithm (TEA), i.e., unicast routing evolution and multicast path composition, for dynamic multicast routing. The TEA uses a novel link-duplicate-degree encoding, which can encode a multicast path in the link-duplicate-degree and decode the path as a link vector easily. A dynamic algorithm for adding nodes to or removing nodes from a multicast group and a repairing algorithm are also covered in this paper. As the TEA is based on global evaluation, the quality of the multicast path remains stabilized without degradation when multicast members change over time. Therefore, it is not necessary to rearrange the multicast path during the life cycle of the multicast sessions. Simulation results show that the TEA is efficient and convergent.
Evolutionary systemic risk: Fisher information flow metric in financial network dynamics
Khashanah, Khaldoun; Yang, Hanchao
2016-03-01
Recently the topic of financial network dynamics has gained renewed interest from researchers in the field of empirical systemic risk measurements. We refer to this type of network analysis as information flow networks analysis (IFNA). This paper proposes a new method that applies Fisher information metric to the evolutionary dynamics of financial networks using IFNA. Our paper is the first to apply the Fisher information metric to a set of financial time series. We introduce Evolution Index (EI) as a measure of systemic risk in financial networks. It is shown, for concrete networks with actual data of several stock markets, that the EI can be implemented as a measure of fitness of the stock market and as a leading indicator of systemic risk.
A New Multiobjective Evolutionary Algorithm for Community Detection in Dynamic Complex Networks
Guoqiang Chen
2013-01-01
Full Text Available Community detection in dynamic networks is an important research topic and has received an enormous amount of attention in recent years. Modularity is selected as a measure to quantify the quality of the community partition in previous detection methods. But, the modularity has been exposed to resolution limits. In this paper, we propose a novel multiobjective evolutionary algorithm for dynamic networks community detection based on the framework of nondominated sorting genetic algorithm. Modularity density which can address the limitations of modularity function is adopted to measure the snapshot cost, and normalized mutual information is selected to measure temporal cost, respectively. The characteristics knowledge of the problem is used in designing the genetic operators. Furthermore, a local search operator was designed, which can improve the effectiveness and efficiency of community detection. Experimental studies based on synthetic datasets show that the proposed algorithm can obtain better performance than the compared algorithms.
Multi CPU clusters and calculations by molecular dynamics method
The technical characteristics of multi CPU (Central Processor Unit) clusters in Institute of Ion-Plasma and Laser Technologies AS RUz and Institute of Mathematics and Information Technologies AS RUz are described. There is detail information about cluster s architecture, installed programs and their productivity for decision of molecular dynamics tasks. Molecular dynamics program packages GROMACS, OPENMX and AutoDock-4.2.3 are described. The results of calculations using these program packages are presented. (author)
Quantum parallelism as a tool for ensemble spin dynamics calculations
Alvarez, Gonzalo A.; Danieli, Ernesto P.; Levstein, Patricia R.; Pastawski, Horacio M.
2007-01-01
Efficient simulations of quantum evolutions of spin-1/2 systems are relevant for ensemble quantum computation as well as in typical NMR experiments. We propose an efficient method to calculate the dynamics of an observable provided that the initial excitation is "local". It resorts a single entangled pure initial state built as a superposition, with random phases, of the pure elements that compose the mixture. This ensures self-averaging of any observable, drastically reducing the calculation...
Self-consistent calculation of spin transport and magnetization dynamics
A spin-polarized current transfers its spin-angular momentum to a local magnetization, exciting various types of current-induced magnetization dynamics. So far, most studies in this field have focused on the direct effect of spin transport on magnetization dynamics, but ignored the feedback from the magnetization dynamics to the spin transport and back to the magnetization dynamics. Although the feedback is usually weak, there are situations when it can play an important role in the dynamics. In such situations, simultaneous, self-consistent calculations of the magnetization dynamics and the spin transport can accurately describe the feedback. This review describes in detail the feedback mechanisms, and presents recent progress in self-consistent calculations of the coupled dynamics. We pay special attention to three representative examples, where the feedback generates non-local effective interactions for the magnetization after the spin accumulation has been integrated out. Possibly the most dramatic feedback example is the dynamic instability in magnetic nanopillars with a single magnetic layer. This instability does not occur without non-local feedback. We demonstrate that full self-consistent calculations generate simulation results in much better agreement with experiments than previous calculations that addressed the feedback effect approximately. The next example is for more typical spin valve nanopillars. Although the effect of feedback is less dramatic because even without feedback the current can make stationary states unstable and induce magnetization oscillation, the feedback can still have important consequences. For instance, we show that the feedback can reduce the linewidth of oscillations, in agreement with experimental observations. A key aspect of this reduction is the suppression of the excitation of short wavelength spin waves by the non-local feedback. Finally, we consider nonadiabatic electron transport in narrow domain walls. The non
Fe IX CALCULATIONS FOR THE SOLAR DYNAMICS OBSERVATORY
New calculations of the energy levels, radiative transition rates, and collisional excitation rates of Fe IX have been carried out using the Flexible Atomic Code, paying close attention to experimentally identified levels and extending existing calculations to higher energy levels. For lower levels, R-matrix collisional excitation rates from earlier work have been used. Significant emission is predicted by these calculations in the 5f-3d transitions, which will impact analysis of Solar Dynamics Observatory Atmospheric Imaging Assembly observations using the 94 A filter.
Evolutionary dynamics of birch (Betula aetnensis Rafin coppices on the Mount Etna (Sicily
Bagnato S
2014-04-01
Full Text Available Evolutionary dynamics of birch (Betula aetnensis Rafin coppices on the Mount Etna (Sicily. The aim of this paper is to evaluate the dynamics of Etna birch stands (Betula aetnensis Rafin following the cessation of silvicultural activities in the Etna Regional Park (Sicily. We investigated forest structure, natural regeneration, vegetation and deadwood in different forest types. Our findings highlighted three different dynamics for birch populations: stable birch stands in the high mountain area which might represent an edapho-climax forest; progressive dynamic birch stands in the intermediate mountain area, showing a gradual depletion of birch and a concomitant replacement with monospecific stands (calabrian pine, beech, oaks or mixed ones (with birch; pure birch stands (typical that tend to be regressive - especially under stressful conditions - and to be replaced by xerophilous grasslands. Following the cessation of coppicing and with stand ageing, the stumps transformation into more homogeneous stand structures have been increasing. Within the context of protected areas the restoration of coppice selection system (with appropriate adaptations could help to maintain the traditional forest landscape, acting as a silvicultural technique with low environmental and landscape impact.
2015-03-01
Full Text Available The dynamic economic load dispatch is one of the main problems of power systems generation and operation. The objective is to schedule power generation for units over a certain period of time, while satisfying operating constraints and load demand in each interval. Wind farms, as renewable energy resources are playing an increasing role in electricity generation. In this paper, a computational framework is presented to solve the dynamic economic emission dispatch problem with inclusion of wind farms considering their associated constraints. An optimization algorithm called modified co-evolutionary particle swarm optimization (MCPSO is proposed to solve the problem. In the proposed algorithm, two kinds of swarms evolve interactively where one of them is used to calculate the penalty factors (constraints handling and the other is used for searching good solutions (optimization process. In addition, some modifications such as using an inertia weight that decreases linearly during the simulation are made to improve the performance of the algorithm. Finally, the validity and superiority of the proposed method are demonstrated by simulation results on a modiﬁed IEEE benchmark system including six thermal units and two wind farms.
Yao Yao
Full Text Available One of the important challenges in the field of evolutionary robotics is the development of systems that can adapt to a changing environment. However, the ability to adapt to unknown and fluctuating environments is not straightforward. Here, we explore the adaptive potential of simulated swarm robots that contain a genomic encoding of a bio-inspired gene regulatory network (GRN. An artificial genome is combined with a flexible agent-based system, representing the activated part of the regulatory network that transduces environmental cues into phenotypic behaviour. Using an artificial life simulation framework that mimics a dynamically changing environment, we show that separating the static from the conditionally active part of the network contributes to a better adaptive behaviour. Furthermore, in contrast with most hitherto developed ANN-based systems that need to re-optimize their complete controller network from scratch each time they are subjected to novel conditions, our system uses its genome to store GRNs whose performance was optimized under a particular environmental condition for a sufficiently long time. When subjected to a new environment, the previous condition-specific GRN might become inactivated, but remains present. This ability to store 'good behaviour' and to disconnect it from the novel rewiring that is essential under a new condition allows faster re-adaptation if any of the previously observed environmental conditions is reencountered. As we show here, applying these evolutionary-based principles leads to accelerated and improved adaptive evolution in a non-stable environment.
Matos Margarida
2009-06-01
Full Text Available Abstract Here we present a correction to our article "Evolutionary dynamics of molecular markers during local adaptation: a case study in Drosophila subobscura". We have recently detected an error concerning the application of the Ln RH formula – a test to detect positive selection – to our microsatellite data. Here we provide the corrected data and discuss its implications for our overall findings. The corrections presented here have produced some changes relative to our previous results, namely in a locus (dsub14 that presents indications of being affected by positive selection. In general, our populations present less consistent indications of positive selection for this particular locus in both periods studied – between generations 3 and 14 and between generation 14 and 40 of laboratory adaptation. Despite this, the main findings of our study regarding the possibility of positive selection acting on that particular microsatellite still hold. As previously concluded in our article, further studies should be performed on this specific microsatellite locus (and neighboring areas to elucidate in greater detail the evolutionary forces acting on this specific region of the O chromosome of Drosophila subobscura.
Bardella, Vanessa Bellini; Fernandes, José Antônio Marin; Cabral-de-Mello, Diogo Cavalcanti
2016-10-01
Previous chromosome mapping of multigene families in Pentatomomorpha (Heteroptera) insects, which was restricted to the major rDNA, revealed remarkable conservation of number of clusters and chromosomal positions. Aiming to understand the chromosomal organization and evolutionary patterns of multigene families in karyotypes of Heteroptera, we performed a chromosomal mapping using four distinct multigene families in representatives of Coreidae (ten species) and Pentatomidae (five species). A single pair of the major rDNA cluster (18S rDNA probe) and a single pair of the minor rDNA cluster (5S rDNA probe), both terminally located were primarily observed, being, in most species, located in distinct chromosomes. However, some alternative patterns were also observed. In species in which the U2 snDNA and H4 gene clusters were mapped, they were mainly located in one autosomal pair each, wherein the H4 gene cluster was located in different positions. Our data suggest that the karyotype diversity reported in Coreidae is not reflected in the distribution diversity of multigene families. This contrasts with the data for Pentatomidae, with a conserved gross karyotype but a discrete diversity in the location of the clusters of multigene families, indicating genome dynamics for these markers. The findings are discussed to shed light on the possible causes for the conservation or variation observed and to assist in understanding the chromosomal evolutionary trends in the group. PMID:27380138
Effect of Recombination in the Evolutionary Dynamics of HIV under the Surveillance of Immune System
Peng, Weiqun; Yang, Wenjing; Wang, Guanyu
2009-03-01
Human immunodeficiency virus (HIV) is a retrovirus that causes acquired immunodeficiency syndrome (AIDS), which has become one of the most destructive pandemics in history. The fact that HIV virus evolves very fast plays a central role in AIDS immunopathogenesis and the difficulty we face in finding a cure or a vaccine for AIDS. A distinguishing feature of HIV is its high frequency of recombination. The effect of recombination in the HIV evolution is not clear. We establish a mathematical model of the evolutionary dynamics. This model incorporates both point mutation and recombination for genetic diversity, and employs a fitness function developed by Wang and Deem (PRL 97, 188106, 2006) that accounts for the effect of immune system. Using this model, we explore the role of recombination in the battle between the virus population and the immune system, with a special focus on the condition under which recombination helps the virus population to escape from the immune system.
Rocha, M. C.; Saraiva, J. T.
2012-10-01
The basic objective of Transmission Expansion Planning (TEP) is to schedule a number of transmission projects along an extended planning horizon minimizing the network construction and operational costs while satisfying the requirement of delivering power safely and reliably to load centres along the horizon. This principle is quite simple, but the complexity of the problem and the impact on society transforms TEP on a challenging issue. This paper describes a new approach to solve the dynamic TEP problem, based on an improved discrete integer version of the Evolutionary Particle Swarm Optimization (EPSO) meta-heuristic algorithm. The paper includes sections describing in detail the EPSO enhanced approach, the mathematical formulation of the TEP problem, including the objective function and the constraints, and a section devoted to the application of the developed approach to this problem. Finally, the use of the developed approach is illustrated using a case study based on the IEEE 24 bus 38 branch test system.
Saraiva J. T.
2012-10-01
Full Text Available The basic objective of Transmission Expansion Planning (TEP is to schedule a number of transmission projects along an extended planning horizon minimizing the network construction and operational costs while satisfying the requirement of delivering power safely and reliably to load centres along the horizon. This principle is quite simple, but the complexity of the problem and the impact on society transforms TEP on a challenging issue. This paper describes a new approach to solve the dynamic TEP problem, based on an improved discrete integer version of the Evolutionary Particle Swarm Optimization (EPSO meta-heuristic algorithm. The paper includes sections describing in detail the EPSO enhanced approach, the mathematical formulation of the TEP problem, including the objective function and the constraints, and a section devoted to the application of the developed approach to this problem. Finally, the use of the developed approach is illustrated using a case study based on the IEEE 24 bus 38 branch test system.
An Evolutionary Algorithm Approach to Link Prediction in Dynamic Social Networks
Bliss, Catherine A; Danforth, Christopher M; Dodds, Peter Sheridan
2013-01-01
Many real world, complex phenomena have underlying structures of evolving networks where nodes and links are added and removed over time. A central scientific challenge is the description and explanation of network dynamics, with a key test being the prediction of short and long term changes. For the problem of short-term link prediction, existing methods attempt to determine neighborhood metrics that correlate with the appearance of a link in the next observation period. Recent work has suggested that the incorporation of user-specific metadata and usage patterns can improve link prediction, however methodologies for doing so in a systematic way are largely unexplored in the literature. Here, we provide an approach to predicting future links by applying an evolutionary algorithm to weights which are used in a linear combination of sixteen neighborhood and node similarity indices. We examine Twitter reciprocal reply networks constructed at the time scale of weeks, both as a test of our general method and as a...
Evolutionary genetic optimization of the injector beam dynamics for the ERL test facility at IHEP
Yi, Jiao
2013-01-01
The energy recovery linac test facility (ERL-TF), a compact ERL-FEL (free electron laser) two-purpose machine, was proposed at the Institute of High Energy Physics, Beijing. As one important component of the ERL-TF, the photo-injector started with a photocathode direct-current gun was designed and preliminarily optimized. In this paper an evolutionary genetic method, non-dominated sorting genetic algorithm II, is applied to optimize the injector beam dynamics, especially in the high-charge operation mode. Study shows that using an incident laser with rms transverse size of 1~1.2 mm, the normalized emittance of the electron beam can be kept below 1 mm.mrad at the end of the injector. This work, together with the previous optimization for the low-charge operation mode by using the iterative scan method, provides guidance and confidence for future constructing and commissioning of the ERL-TF injector.
VERIFICATION OF TORSIONAL OSCILLATING MECHANICAL SYSTEM DYNAMIC CALCULATION RESULTS
Peter KAŠŠAY
2014-09-01
Full Text Available On our department we deal with optimization and tuning of torsional oscillating mechanical systems. When solving these problems we often use results of dynamic calculation. The goal of this article is to compare values obtained by computation and experimentally. For this purpose, a mechanical system built in our laboratory was used. At first, classical HARDY type flexible coupling has been applied into the system, then we used a pneumatic flexible shaft coupling developed by us. The main difference of these couplings over conventional flexible couplings is that they can change their dynamic properties during operation, by changing the pressure of the gaseous medium in their flexible elements.
Dissipative Particle Dynamics interaction parameters from ab initio calculations
Sepehr, Fatemeh; Paddison, Stephen J.
2016-02-01
Dissipative Particle Dynamics (DPD) is a commonly employed coarse-grained method to model complex systems. Presented here is a pragmatic approach to connect atomic-scale information to the meso-scale interactions defined between the DPD particles or beads. Specifically, electronic structure calculations were utilized for the calculation of the DPD pair-wise interaction parameters. An implicit treatment of the electrostatic interactions for charged beads is introduced. The method is successfully applied to derive the parameters for a hydrated perfluorosulfonic acid ionomer with absorbed vanadium cations.
Quantum Parallelism as a Tool for Ensemble Spin Dynamics Calculations
Álvarez, Gonzalo A.; Danieli, Ernesto P.; Levstein, Patricia R.; Pastawski, Horacio M.
2008-09-01
Efficient simulations of quantum evolutions of spin-1/2 systems are relevant for ensemble quantum computation as well as in typical NMR experiments. We propose an efficient method to calculate the dynamics of an observable provided that the initial excitation is “local.” It resorts to a single entangled pure initial state built as a superposition, with random phases, of the pure elements that compose the mixture. This ensures self-averaging of any observable, drastically reducing the calculation time. The procedure is tested for two representative systems: a spin star (cluster with random long range interactions) and a spin ladder.
Tian, Yihui; Govindan, Kannan; Zhu, Qinghua
2014-01-01
In this study, a system dynamics (SD) model is developed to guide the subsidy policies to promote the diffusion of green supply chain management (GSCM) in China. The relationships of stakeholders such as government, enterprises and consumers are analyzed through evolutionary game theory. Finally...
Calculation of the dynamic air flow resistivity of fibre materials
Tarnow, Viggo
1997-01-01
The acoustic attenuation of acoustic fiber materials is mainly determined by the dynamic resistivity to an oscillating air flow. The dynamic resistance is calculated for a model with geometry close to the geometry of real fibre material. The model constists of parallel cylinders placed randomly....... Two case are treated: flow perpendicular to the cylinder axes, and flow parallel to the axes. In each case two new approximate procedures were used. In the first procedure, one solves the equation of flow in a Voronoi cell around the fiber, and averages over the distribution of the Voronoi cells.......The second procedure is an extension to oscillating air flow of the Brinkman self-consistent procedure for dc flow. The procedures are valid for volume concentrations of cylinders less than 0.1. The calculations show that for the density of fibers of interest for acoustic fibre materials the simple self...
Andersen, Christian Walther; Sibani, Paolo
2016-05-01
Based on the stochastic dynamics of interacting agents which reproduce, mutate, and die, the tangled nature model (TNM) describes key emergent features of biological and cultural ecosystems' evolution. While trait inheritance is not included in many applications, i.e., the interactions of an agent and those of its mutated offspring are taken to be uncorrelated, in the family of TNMs introduced in this work correlations of varying strength are parametrized by a positive integer K . We first show that the interactions generated by our rule are nearly independent of K . Consequently, the structural and dynamical effects of trait inheritance can be studied independently of effects related to the form of the interactions. We then show that changing K strengthens the core structure of the ecology, leads to population abundance distributions better approximated by log-normal probability densities, and increases the probability that a species extant at time tw also survives at t >tw . Finally, survival probabilities of species are shown to decay as powers of the ratio t /tw , a so-called pure aging behavior usually seen in glassy systems of physical origin. We find a quantitative dynamical effect of trait inheritance, namely, that increasing the value of K numerically decreases the decay exponent of the species survival probability.
Tennessen, Jacob A; Govindarajulu, Rajanikanth; Ashman, Tia-Lynn; Liston, Aaron
2014-12-01
Whole-genome duplications are radical evolutionary events that have driven speciation and adaptation in many taxa. Higher-order polyploids have complex histories often including interspecific hybridization and dynamic genomic changes. This chromosomal reshuffling is poorly understood for most polyploid species, despite their evolutionary and agricultural importance, due to the challenge of distinguishing homologous sequences from each other. Here, we use dense linkage maps generated with targeted sequence capture to improve the diploid strawberry (Fragaria vesca) reference genome and to disentangle the subgenomes of the wild octoploid progenitors of cultivated strawberry, Fragaria virginiana and Fragaria chiloensis. Our novel approach, POLiMAPS (Phylogenetics Of Linkage-Map-Anchored Polyploid Subgenomes), leverages sequence reads to associate informative interhomeolog phylogenetic markers with linkage groups and reference genome positions. In contrast to a widely accepted model, we find that one of the four subgenomes originates with the diploid cytoplasm donor F. vesca, one with the diploid Fragaria iinumae, and two with an unknown ancestor close to F. iinumae. Extensive unidirectional introgression has converted F. iinumae-like subgenomes to be more F. vesca-like, but never the reverse, due either to homoploid hybridization in the F. iinumae-like diploid ancestors or else strong selection spreading F. vesca-like sequence among subgenomes through homeologous exchange. In addition, divergence between homeologous chromosomes has been substantially augmented by interchromosomal rearrangements. Our phylogenetic approach reveals novel aspects of the complicated web of genetic exchanges that occur during polyploid evolution and suggests a path forward for unraveling other agriculturally and ecologically important polyploid genomes. PMID:25477420
Helbing, Dirk; Johansson, Anders
2010-01-01
Evolutionary game theory has been successfully used to investigate the dynamics of systems, in which many entities have competitive interactions. From a physics point of view, it is interesting to study conditions under which a coordination or cooperation of interacting entities will occur, be it spins, particles, bacteria, animals, or humans. Here, we analyze the case, where the entities are heterogeneous, particularly the case of two populations with conflicting interactions and two possible states. For such systems, explicit mathematical formulas will be determined for the stationary solutions and the associated eigenvalues, which determine their stability. In this way, four different types of system dynamics can be classified and the various kinds of phase transitions between them will be discussed. While these results are interesting from a physics point of view, they are also relevant for social, economic, and biological systems, as they allow one to understand conditions for (1) the breakdown of cooperation, (2) the coexistence of different behaviors (“subcultures”), (3) the evolution of commonly shared behaviors (“norms”), and (4) the occurrence of polarization or conflict. We point out that norms have a similar function in social systems that forces have in physics.
Jingjing Ma
2014-01-01
Full Text Available Community structure is one of the most important properties in social networks. In dynamic networks, there are two conflicting criteria that need to be considered. One is the snapshot quality, which evaluates the quality of the community partitions at the current time step. The other is the temporal cost, which evaluates the difference between communities at different time steps. In this paper, we propose a decomposition-based multiobjective community detection algorithm to simultaneously optimize these two objectives to reveal community structure and its evolution in dynamic networks. It employs the framework of multiobjective evolutionary algorithm based on decomposition to simultaneously optimize the modularity and normalized mutual information, which quantitatively measure the quality of the community partitions and temporal cost, respectively. A local search strategy dealing with the problem-specific knowledge is incorporated to improve the effectiveness of the new algorithm. Experiments on computer-generated and real-world networks demonstrate that the proposed algorithm can not only find community structure and capture community evolution more accurately, but also be steadier than the two compared algorithms.
Random and non-random mating populations: Evolutionary dynamics in meiotic drive.
Sarkar, Bijan
2016-01-01
Game theoretic tools are utilized to analyze a one-locus continuous selection model of sex-specific meiotic drive by considering nonequivalence of the viabilities of reciprocal heterozygotes that might be noticed at an imprinted locus. The model draws attention to the role of viability selections of different types to examine the stable nature of polymorphic equilibrium. A bridge between population genetics and evolutionary game theory has been built up by applying the concept of the Fundamental Theorem of Natural Selection. In addition to pointing out the influences of male and female segregation ratios on selection, configuration structure reveals some noted results, e.g., Hardy-Weinberg frequencies hold in replicator dynamics, occurrence of faster evolution at the maximized variance fitness, existence of mixed Evolutionarily Stable Strategy (ESS) in asymmetric games, the tending evolution to follow not only a 1:1 sex ratio but also a 1:1 different alleles ratio at particular gene locus. Through construction of replicator dynamics in the group selection framework, our selection model introduces a redefining bases of game theory to incorporate non-random mating where a mating parameter associated with population structure is dependent on the social structure. Also, the model exposes the fact that the number of polymorphic equilibria will depend on the algebraic expression of population structure. PMID:26524140
Mia T Levine
Full Text Available Heterochromatin is the gene-poor, satellite-rich eukaryotic genome compartment that supports many essential cellular processes. The functional diversity of proteins that bind and often epigenetically define heterochromatic DNA sequence reflects the diverse functions supported by this enigmatic genome compartment. Moreover, heterogeneous signatures of selection at chromosomal proteins often mirror the heterogeneity of evolutionary forces that act on heterochromatic DNA. To identify new such surrogates for dissecting heterochromatin function and evolution, we conducted a comprehensive phylogenomic analysis of the Heterochromatin Protein 1 gene family across 40 million years of Drosophila evolution. Our study expands this gene family from 5 genes to at least 26 genes, including several uncharacterized genes in Drosophila melanogaster. The 21 newly defined HP1s introduce unprecedented structural diversity, lineage-restriction, and germline-biased expression patterns into the HP1 family. We find little evidence of positive selection at these HP1 genes in both population genetic and molecular evolution analyses. Instead, we find that dynamic evolution occurs via prolific gene gains and losses. Despite this dynamic gene turnover, the number of HP1 genes is relatively constant across species. We propose that karyotype evolution drives at least some HP1 gene turnover. For example, the loss of the male germline-restricted HP1E in the obscura group coincides with one episode of dramatic karyotypic evolution, including the gain of a neo-Y in this lineage. This expanded compendium of ovary- and testis-restricted HP1 genes revealed by our study, together with correlated gain/loss dynamics and chromosome fission/fusion events, will guide functional analyses of novel roles supported by germline chromatin.
Electronic Structure and Molecular Dynamics Calculations for KBH4
Papaconstantopoulos, Dimitrios; Shabaev, Andrew; Hoang, Khang; Mehl, Michael; Kioussis, Nicholas
2012-02-01
In the search for hydrogen storage materials, alkali borohydrides MBH4 (M=Li, Na, K) are especially interesting because of their light weight and the high number of hydrogen atoms per metal atom. Electronic structure calculations can give insights into the properties of these complex hydrides and provide understanding of the structural properties and of the bonding of hydrogen. We have performed first-principles density-functional theory (DFT) and tight-binding (TB) calculations for KBH4 in both the high temperature (HT) and low temperature (LT) phases to understand its electronic and structural properties. Our DFT calculations were carried out using the VASP code. The results were then used as a database to develop a tight-binding Hamiltonian using the NRL-TB method. This approach allowed for computationally efficient calculations of phonon frequencies and elastic constants using the static module of the NRL-TB, and also using the molecular dynamics module to calculate mean-square displacements and formation energies of hydrogen vacancies.
A cascade of evolutionary change alters consumer-resource dynamics and ecosystem function
Walsh, Matthew R.; DeLong, John P.; Hanley, Torrance C.; Post, David M
2012-01-01
It is becoming increasingly clear that intraspecific evolutionary divergence influences the properties of populations, communities and ecosystems. The different ecological impacts of phenotypes and genotypes may alter selection on many species and promote a cascade of ecological and evolutionary change throughout the food web. Theory predicts that evolutionary interactions across trophic levels may contribute to hypothesized feedbacks between ecology and evolution. However, the importance of ...
Mueller, Georg P.
2012-01-01
This article analyzes different cultures of corruption with regard to their evolutionary stability, i.e. their ability to annihilate small disturbances in the equilibria between corrupt and noncorrupt agents. The article starts with the development of an evolutionary model of the interactions between corrupt and noncorrupt citizens and functionaries of the state, which is subsequently explored by formal analyses and computer simulation. It turns out that zero-corruption is always evolutionari...
VTT Energy's calculation system for reactor physics and dynamics
VTT Energy has a comprehensive and independent calculation system for reactor physics and dynamics analyses. The system is widely utilized in contract research for the nuclear safety authorities and power companies. The four Finnish reactors, some foreign plants and potential new plant concepts have been studied. The system is being modernized all the time both by own development work and by international cooperation. The reactor physics codes are upgraded and new codes and methods are developed and acquired for calculations and safety evaluations of new, increasingly complicated fuel assembly types and fuel-loading schemes, as well as for criticality and dose rate studies. The reactor dynamics methods are developed and new sophisticated models are created for tasks related to increased safety requirements. The primary aim is to realistically simulate reactor stability and complicated reactivity accidents with three-dimensional core models. For thermal hydraulics calculations, an accurate general flow model based on a new solution method has been developed. (orig.) (30 refs., 3 figs.)
METHOD FOR CALCULATION OF STRESSED STATE SUBSTANTIATED BY DYNAMIC MICROTWIN
V. V. Vlashevich
2015-01-01
Full Text Available Method for calculation of the stressed state in a dynamic twin has been developed on the basis of a non-thin non-coherent micro-twin model with continuous distribution of twinning dislocations at twin boundaries. In this case there is no additional generation with the help of twinning dislocation source. The model takes into account that the twin has coherent and noncoherent boundary sections. The developed model has made it possible to take into consideration a form of non-coherent sections of twinning boundaries in calculations of stressed and deformed state at dynamic twins. It has been established that localized stresses are migrating together with non-coherent sections of the twin. Normal stresses σxx change their sign in relation to direction of the twin development. Shear stresses σxy are alternating in signs in relation to an axis which is perpendicular to the direction of the twin development and which is passing through a mid-point of non-coherent twin section. Distribution of stresses σyy и σyz has similar configuration. Stresses σzx in the second and fourth quarters of XOY plane are negative and the stresses in the first and third quarters are positive. Distribution of stresses σzz practically does not differ from distribution of stresses σyy according to configuration but numerical values of stress tensor component data are different.The results have been obtained without thin twin model that permits to consider only elastic stage of the twinning process. The executed stress calculations at dynamic twin are important for forecasting at the accumulation stage of damage origination which is caused by twinning destruction and permit to improve forecasting accuracy of technical system resources on the basis of twinning materials such as alloys based on iron, copper, zinc, aluminium, titanium.
Advanced Dynamics Analytical and Numerical Calculations with MATLAB
Marghitu, Dan B
2012-01-01
Advanced Dynamics: Analytical and Numerical Calculations with MATLAB provides a thorough, rigorous presentation of kinematics and dynamics while using MATLAB as an integrated tool to solve problems. Topics presented are explained thoroughly and directly, allowing fundamental principles to emerge through applications from areas such as multibody systems, robotics, spacecraft and design of complex mechanical devices. This book differs from others in that it uses symbolic MATLAB for both theory and applications. Special attention is given to solutions that are solved analytically and numerically using MATLAB. The illustrations and figures generated with MATLAB reinforce visual learning while an abundance of examples offer additional support. This book also: Provides solutions analytically and numerically using MATLAB Illustrations and graphs generated with MATLAB reinforce visual learning for students as they study Covers modern technical advancements in areas like multibody systems, robotics, spacecraft and des...
Evolutionary history and spatiotemporal dynamics of dengue virus type 1 in Asia.
Sun, Yan; Meng, Shengli
2013-06-01
Previous studies showed that DENV-1 transmitted from monkeys to humans approximately 125 years ago. However, there is no comprehensive analysis about phylogeography and population dynamics of Asian DENV-1. Here, we adopt a Bayesian phylogeographic approach to investigate the evolutionary history and phylogeography of Asian DENV-1 using envelope (E) protein gene sequences of 450 viruses isolated from 1954 to 2010 throughout 18 Asian countries and regions. Bayesian phylogeographic analyses indicate that the high rates of viral migration possibly follows long-distance travel for humans in Southeast Asia. Our study highlights that Southeast Asian countries have acted as the main viral sources of the dengue epidemics in East Asia. The results reveal that the time to the most recent common ancestor (TMRCA) of Asian DENV-1 is 1906 (95% HPD, years 1897-1915). We show that the spatial dissemination of virus is the major source of DENV-1 outbreaks in the different localities and leads to subsequent establishment and expansion of the virus in these areas. PMID:23395769
Evolutionary dynamics of strategic behavior in a collective-risk dilemma.
Abou Chakra, Maria; Traulsen, Arne
2012-01-01
A collective-risk social dilemma arises when a group must cooperate to reach a common target in order to avoid the risk of collective loss while each individual is tempted to free-ride on the contributions of others. In contrast to the prisoners' dilemma or public goods games, the collective-risk dilemma encompasses the risk that all individuals lose everything. These characteristics have potential relevance for dangerous climate change and other risky social dilemmas. Cooperation is costly to the individual and it only benefits all individuals if the common target is reached. An individual thus invests without guarantee that the investment is worthwhile for anyone. If there are several subsequent stages of investment, it is not clear when individuals should contribute. For example, they could invest early, thereby signaling their willingness to cooperate in the future, constantly invest their fair share, or wait and compensate missing contributions. To investigate the strategic behavior in such situations, we have simulated the evolutionary dynamics of such collective-risk dilemmas in a finite population. Contributions depend individually on the stage of the game and on the sum of contributions made so far. Every individual takes part in many games and successful behaviors spread in the population. It turns out that constant contributors, such as constant fair sharers, quickly lose out against those who initially do not contribute, but compensate this in later stages of the game. In particular for high risks, such late contributors are favored. PMID:22927807
Evolutionary dynamics of strategic behavior in a collective-risk dilemma.
Maria Abou Chakra
Full Text Available A collective-risk social dilemma arises when a group must cooperate to reach a common target in order to avoid the risk of collective loss while each individual is tempted to free-ride on the contributions of others. In contrast to the prisoners' dilemma or public goods games, the collective-risk dilemma encompasses the risk that all individuals lose everything. These characteristics have potential relevance for dangerous climate change and other risky social dilemmas. Cooperation is costly to the individual and it only benefits all individuals if the common target is reached. An individual thus invests without guarantee that the investment is worthwhile for anyone. If there are several subsequent stages of investment, it is not clear when individuals should contribute. For example, they could invest early, thereby signaling their willingness to cooperate in the future, constantly invest their fair share, or wait and compensate missing contributions. To investigate the strategic behavior in such situations, we have simulated the evolutionary dynamics of such collective-risk dilemmas in a finite population. Contributions depend individually on the stage of the game and on the sum of contributions made so far. Every individual takes part in many games and successful behaviors spread in the population. It turns out that constant contributors, such as constant fair sharers, quickly lose out against those who initially do not contribute, but compensate this in later stages of the game. In particular for high risks, such late contributors are favored.
The effect of network structure on innovation initiation process: an evolutionary dynamics approach
Jafari, Afshin; Zolfagharzadeh, Mohammad Mahdi; Mohammadi, Mehdi
2016-01-01
In this paper we have proposed a basic agent-based model based on evolutionary dynamics for investigating innovation initiation process. In our model we suppose each agent will represent a firm which is interacting with other firms through a given network structure. We consider a two-hit process for presenting a potentially successful innovation in this model and therefore at each time step each firm can be in on of three different stages which are respectively, Ordinary, Innovative, and Successful. We design different experiments in order to investigate how different interaction networks may affect the process of presenting a successful innovation to the market. In this experiments, we use five different network structures, i.e. Erd\\H{o}s and R\\'enyi, Ring Lattice, Small World, Scale-Free and Distance-Based networks. According to the results of the simulations, for less frequent innovations like radical innovation, local structures are showing a better performance comparing to Scale-Free and Erd\\H{o}s and R\\...
Evolutionary Dynamics of MERS-CoV: Potential Recombination, Positive Selection and Transmission.
Zhang, Zhao; Shen, Libing; Gu, Xun
2016-01-01
Middle East respiratory syndrome coronavirus (MERS-CoV) belongs to beta group of coronavirus and was first discovered in 2012. MERS-CoV can infect multiple host species and cause severe diseases in human. We conducted a series of phylogenetic and bioinformatic analyses to study the evolution dynamics of MERS-CoV among different host species with genomic data. Our analyses show: 1) 28 potential recombinant sequences were detected and they can be classified into seven potential recombinant types; 2) The spike (S) protein of MERS-CoV was under strong positive selection when MERS-CoV transmitted from their natural host to human; 3) Six out of nine positive selection sites detected in spike (S) protein are located in its receptor-binding domain which is in direct contact with host cells; 4) MERS-CoV frequently transmitted back and forth between human and camel after it had acquired the human-camel infection capability. Together, these results suggest that potential recombination events might have happened frequently during MERS-CoV's evolutionary history and the positive selection sites in MERS-CoV's S protein might enable it to infect human. PMID:27142087
The evolutionary dynamics of the Helena retrotransposon revealed by sequenced Drosophila genomes
Carareto Claudia MA
2009-07-01
Full Text Available Abstract Background Several studies have shown that genomes contain a mixture of transposable elements, some of which are still active and others ancient relics that have degenerated. This is true for the non-LTR retrotransposon Helena, of which only degenerate sequences have been shown to be present in some species (Drosophila melanogaster, whereas putatively active sequences are present in others (D. simulans. Combining experimental and population analyses with the sequence analysis of the 12 Drosophila genomes, we have investigated the evolution of Helena, and propose a possible scenario for the evolution of this element. Results We show that six species of Drosophila have the Helena transposable element at different stages of its evolution. The copy number is highly variable among these species, but most of them are truncated at the 5' ends and also harbor several internal deletions and insertions suggesting that they are inactive in all species, except in D. mojavensis in which quantitative RT-PCR experiments have identified a putative active copy. Conclusion Our data suggest that Helena was present in the common ancestor of the Drosophila genus, which has been vertically transmitted to the derived lineages, but that it has been lost in some of them. The wide variation in copy number and sequence degeneration in the different species suggest that the evolutionary dynamics of Helena depends on the genomic environment of the host species.
Dynamic calculation of structures in seismic zones. 2. ed.
The aims of this book are both didactic and practical. It is therefore addressed to both experienced engineers and students. Some general information about earthquakes and their occurrence is first given. The problem of a simple oscillator is presented. In this way, the reader is provided with an insight into undestanding the dynamic phenomena taking place and is introduced to the concept of response spectra and to an intuitive comprehension of the behavior of structures during earthquakes. The next chapter is devoted to the cases most frequently encountered with multiple oscillator structures. Theoretical studies are based on the usual modal decomposition method. The various practical methods of calculation employed are then examined, emphasis being given to the various different stages involved and to which of them is the best suited for a particular type of structure. Advise is given on how to select the model whose behavior best describes the real structure, both manual and computer methods of calculation being envisaged
Molecular dynamics calculations for sodium using pseudopotential theory
The equation of state of sodium is studied using the molecular dynamics technique whereby the classical motion of a system of ions is solved with the aid of computers. The interaction potential between pairs of sodium ions consists of Coulomb and Born-Mayer repulsion terms and an effective ion-ion interaction derived from pseudopotential theory. This theory includes the effects of electron gas screening, exchange, and correlation. A model pseudopotential with parameters fit to experimental low-temperature data is used. By using this technique, an atomic description of a simple metal proceeds to the calculation of macroscopic thermodynamic properties
Approximate dynamic fault tree calculations for modelling water supply risks
Traditional fault tree analysis is not always sufficient when analysing complex systems. To overcome the limitations dynamic fault tree (DFT) analysis is suggested in the literature as well as different approaches for how to solve DFTs. For added value in fault tree analysis, approximate DFT calculations based on a Markovian approach are presented and evaluated here. The approximate DFT calculations are performed using standard Monte Carlo simulations and do not require simulations of the full Markov models, which simplifies model building and in particular calculations. It is shown how to extend the calculations of the traditional OR- and AND-gates, so that information is available on the failure probability, the failure rate and the mean downtime at all levels in the fault tree. Two additional logic gates are presented that make it possible to model a system's ability to compensate for failures. This work was initiated to enable correct analyses of water supply risks. Drinking water systems are typically complex with an inherent ability to compensate for failures that is not easily modelled using traditional logic gates. The approximate DFT calculations are compared to results from simulations of the corresponding Markov models for three water supply examples. For the traditional OR- and AND-gates, and one gate modelling compensation, the errors in the results are small. For the other gate modelling compensation, the error increases with the number of compensating components. The errors are, however, in most cases acceptable with respect to uncertainties in input data. The approximate DFT calculations improve the capabilities of fault tree analysis of drinking water systems since they provide additional and important information and are simple and practically applicable.
Dynamical calculations of nuclear fission and heavy-ion reactions
With the goal of determining the magnitude and mechanism of nuclear dissipation from comparisons of predictions with experimental data, we describe recent calculations in a unified macroscopic-microscopic approach to large-amplitude collective nuclear motion such as occurs in fission and heavy-ion reactions. We describe the time dependence of the distribution function in phase space of collective coordinates and momenta by a generalized Fokker-Planck equation. The nuclear potential energy of deformation is calculated as the sum of repulsive Coulomb and centrifugal energies and an attractive Yukawa-plus-exponential potential, the inertia tensor is calculated for a superposition of rigid-body rotation and incompressible, nearly irrotational flow by use of the Werner-Wheeler method, and the dissipation ensor that describes the conversion of collective energy into single-particle excitation energy is calculated for two prototype mechanisms that represent opposite extremes of large and small dissipation. We solve the generalized Hamilton equations of motion for the first moments of the distribution function to obtain the mean translational fission-fragment kinetic energy and mass of a third fragment that sometimes forms between the two end fragments, as well as dynamical thresholds, capture cross sections, and ternary events in heavy-ion reactions. 33 references
Neotropical fish-fruit interactions: eco-evolutionary dynamics and conservation.
Correa, Sandra Bibiana; Costa-Pereira, Raul; Fleming, Theodore; Goulding, Michael; Anderson, Jill T
2015-11-01
Frugivorous fish play a prominent role in seed dispersal and reproductive dynamics of plant communities in riparian and floodplain habitats of tropical regions worldwide. In Neotropical wetlands, many plant species have fleshy fruits and synchronize their fruiting with the flood season, when fruit-eating fish forage in forest and savannahs for periods of up to 7 months. We conducted a comprehensive analysis to examine the evolutionary origin of fish-fruit interactions, describe fruit traits associated with seed dispersal and seed predation, and assess the influence of fish size on the effectiveness of seed dispersal by fish (ichthyochory). To date, 62 studies have documented 566 species of fruits and seeds from 82 plant families in the diets of 69 Neotropical fish species. Fish interactions with flowering plants are likely to be as old as 70 million years in the Neotropics, pre-dating most modern bird-fruit and mammal-fruit interactions, and contributing to long-distance seed dispersal and possibly the radiation of early angiosperms. Ichthyochory occurs across the angiosperm phylogeny, and is more frequent among advanced eudicots. Numerous fish species are capable of dispersing small seeds, but only a limited number of species can disperse large seeds. The size of dispersed seeds and the probability of seed dispersal both increase with fish size. Large-bodied species are the most effective seed dispersal agents and remain the primary target of fishing activities in the Neotropics. Thus, conservation efforts should focus on these species to ensure continuity of plant recruitment dynamics and maintenance of plant diversity in riparian and floodplain ecosystems. PMID:25599800
Agoti, Charles N.; Gitahi, Caroline W.; Bett, Ann; Ngama, Mwanajuma; Medley, Graham F.; Cane, Patricia A.; Nokes, D. James
2016-01-01
ABSTRACT The characteristic recurrent epidemics of human respiratory syncytial virus (RSV) within communities may result from the genetic variability of the virus and associated evolutionary adaptation, reducing the efficiency of preexisting immune responses. We analyzed the molecular evolutionary changes in the attachment (G) glycoprotein of RSV-A viruses collected over 13 epidemic seasons (2000 to 2012) in Kilifi (n = 649), Kenya, and contemporaneous sequences (n = 1,131) collected elsewhere within Kenya and 28 other countries. Genetic diversity in the G gene in Kilifi was dynamic both within and between epidemics, characterized by frequent new variant introductions and limited variant persistence between consecutive epidemics. Four RSV-A genotypes were detected in Kilifi: ON1 (11.9%), GA2 (75.5%), GA5 (12.3%), and GA3 (0.3%), with predominant genotype replacement of GA5 by GA2 and then GA2 by ON1. Within these genotypes, there was considerable variation in potential N-glycosylation sites, with GA2 and ON1 viruses showing up to 15 different patterns involving eight possible sites. Further, we identified 15 positively selected and 34 genotype-distinguishing codon sites, with six of these sites exhibiting both characteristics. The mean substitution rate of the G ectodomain for the Kilifi data set was estimated at 3.58 × 10−3 (95% highest posterior density interval = 3.04 to 4.16) nucleotide substitutions/site/year. Kilifi viruses were interspersed in the global phylogenetic tree, clustering mostly with Kenyan and European sequences. Our findings highlight ongoing genetic evolution and high diversity of circulating RSV-A strains, locally and globally, with potential antigenic differences. Taken together, these provide a possible explanation on the nature of recurrent local RSV epidemics. IMPORTANCE The mechanisms underlying recurrent epidemics of RSV are poorly understood. We observe high genetic diversity in circulating strains within and between epidemics in
José Alexandre F. Diniz-Filho
2013-10-01
Full Text Available Macroecology focuses on ecological questions at broad spatial and temporal scales, providing a statistical description of patterns in species abundance, distribution and diversity. More recently, historical components of these patterns have begun to be investigated more deeply. We tentatively refer to the practice of explicitly taking species history into account, both analytically and conceptually, as ‘evolutionary macroecology’. We discuss how the evolutionary dimension can be incorporated into macroecology through two orthogonal and complementary data types: fossils and phylogenies. Research traditions dealing with these data have developed more‐or‐less independently over the last 20–30 years, but merging them will help elucidate the historical components of diversity gradients and the evolutionary dynamics of species’ traits. Here we highlight conceptual and methodological advances in merging these two research traditions and review the viewpoints and toolboxes that can, in combination, help address patterns and unveil processes at temporal and spatial macro‐scales.
In this study either cluster fragmentation, using a time-dependent Hartree-Fock formulation, or cluster deposition, based on classical molecular dynamics, have been studied. An exhaustive analysis has been performed on the many parameters acting on the two processes. Fragmentation calculations show a primary dependence on the input energy whereas the interatomic forces play a primary role in deposition. However the central result of this study is the essential agreement between the classical and quantum mechanical calculation
Biopolymers under large external forces and mean-field RNA virus evolutionary dynamics
Ahsan, Syed Amir
The modeling of the mechanical response of single-molecules of DNA and RNA under large external forces through statistical mechanical methods is central to this thesis with a small portion devoted to modeling the evolutionary dynamics of positive-sense single-stranded RNA viruses. In order to develop and test models of biopolymer mechanics and illuminate the mechanisms underlying biological processes where biopolymers undergo changes in energy on the order of the thermal energy, , entails measuring forces and lengths on the scale of piconewtons (pN) and nanometers (nm), respectively. A capacity achieved in the past two decades at the single-molecule level through the development of micromanipulation techniques such as magnetic and optical tweezers, atomic force microscopy, coupled with advances in micro- and nanofabrication. The statistical mechanical models of biopolymers developed in this dissertation are dependent upon and the outcome of these advancements and resulting experiments. The dissertation begins in chapter 1 with an introduction to the structure and thermodynamics of DNA and RNA, highlighting the importance and effectiveness of simple, two-state models in their description as a prelude to the emergence of two-state models in the research manuscripts. In chapter 2 the standard models of the elasticity of polymers and of a polymer gel are reviewed, characterizing the continuum and mean-field models, including the scaling behavior of DNA in confined spaces. The research manuscript presented in the last section of chapter 2 (section 2.5), subsequent to a review of a Flory gel and in contrast to it, is a model of the elasticity of RNA as a gel, with viral RNA illustrating an instance of such a network, and shown to exhibit anomalous elastic behavior, a negative Poisson ratio, and capable of facilitating viral RNA encapsidation with further context provided in section 5.1. In chapter 3 the experimental methods and behavior of DNA and RNA under mechanical
Wolfe, Alexander P.; Siver, Peter A.
2013-12-01
Chrysophyte algae are common in the plankton of oligotrophic lakes and produce a rich microfossil record of siliceous cysts and scales. Paleolimnological investigations and phytoplankton records suggest that chrysophyte populations are increasing in a wide range of boreal and arctic lakes, ultimately representing one component of the limnological response to contemporary global changes. However, the exact mechanisms responsible for widespread increases of chrysophyte populations remain elusive. We hypothesize that recent increases in chrysophytes are related to rising pCO2 in lakes, in part because these algae lack carbon concentrating mechanisms and therefore rely on diffusive entry of CO2 to Rubisco during photosynthesis. We assessed the abundance of modern sediment chrysophyte microfossils in relation to summer CO2 relative saturation in 46 New England (USA) lakes, revealing significant positive relationships for both cysts and scales. These observations imply that correlations between chrysophytes and limnological conditions including low pH, oligotrophy, and elevated dissolved organic matter are ultimately underscored by the high pCO2 associated with these conditions. In lakes where chrysophyte populations have expanded over recent decades, we infer that increasingly heterotrophic conditions with respect to CO2 have stimulated production by these organisms. This linkage is supported by the remarkable abundance and diversity of chrysophytes from middle Eocene lake sediments, deposited under atmospheric CO2 concentrations significantly higher than present. The Eocene assemblages suggest that any chrysophyte-CO2 connection borne out of results from modern and sub-recent sediments also operated on evolutionary time scales, and thus the absence of carbon concentrating mechanisms appears to be an ancient feature within the group. Chrysophyte microfossils may potentially provide important insights concerning the temporal dynamics of carbon cycling in aquatic
Beam dynamics calculations and particle tracking using massively parallel processors
During the past decade massively parallel processors (MPPs) have slowly gained acceptance within the scientific community. At present these machines typically contain a few hundred to one thousand off-the-shelf microprocessors and a total memory of up to 32 GBytes. The potential performance of these machines is illustrated by the fact that a month long job on a high end workstation might require only a few hours on an MPP. The acceptance of MPPs has been slow for a variety of reasons. For example, some algorithms are not easily parallelizable. Also, in the past these machines were difficult to program. But in recent years the development of Fortran-like languages such as CM Fortran and High Performance Fortran have made MPPs much easier to use. In the following we will describe how MPPs can be used for beam dynamics calculations and long term particle tracking
Flux-vector splitting for unsteady calculations on dynamic meshes
Anderson, W. Kyle; Thomas, James L.; Rumsey, Christopher L.
1989-01-01
The method of flux vector splitting used is that of Van Leer. The fluxes split in this manner have the advantage of being continuously differentiable at eigenvalue sign changes and this allows normal shocks to be captured with at most two interior zones, although in practice only one zone is usually observed. The fluxes as originally derived, however did not include the necessary terms appropriate for calculations on a dynamic mesh. The extension of the splitting to include these terms while retaining the advantages of the original splitting is the main purpose of this investigation. In addition, the use of multiple grids to reduce the computer time is investigated. A subiterative procedure to eliminate factorization and linearization error so that larger time steps can be used is also investigated.
Automating the parallel processing of fluid and structural dynamics calculations
Arpasi, Dale J.; Cole, Gary L.
1987-01-01
The NASA Lewis Research Center is actively involved in the development of expert system technology to assist users in applying parallel processing to computational fluid and structural dynamic analysis. The goal of this effort is to eliminate the necessity for the physical scientist to become a computer scientist in order to effectively use the computer as a research tool. Programming and operating software utilities have previously been developed to solve systems of ordinary nonlinear differential equations on parallel scalar processors. Current efforts are aimed at extending these capabilties to systems of partial differential equations, that describe the complex behavior of fluids and structures within aerospace propulsion systems. This paper presents some important considerations in the redesign, in particular, the need for algorithms and software utilities that can automatically identify data flow patterns in the application program and partition and allocate calculations to the parallel processors. A library-oriented multiprocessing concept for integrating the hardware and software functions is described.
Epidemiological and Evolutionary Dynamics of Influenza B Viruses in Malaysia, 2012-2014.
Xiang Yong Oong
Full Text Available Epidemiological and evolutionary dynamics of influenza B Victoria and Yamagata lineages remained poorly understood in the tropical Southeast Asia region, despite causing seasonal outbreaks worldwide. From 2012-2014, nasopharyngeal swab samples collected from outpatients experiencing acute upper respiratory tract infection symptoms in Kuala Lumpur, Malaysia, were screened for influenza viruses using a multiplex RT-PCR assay. Among 2,010/3,935 (51.1% patients infected with at least one respiratory virus, 287 (14.3% and 183 (9.1% samples were tested positive for influenza A and B viruses, respectively. Influenza-positive cases correlate significantly with meteorological factors-total amount of rainfall, relative humidity, number of rain days, ground temperature and particulate matter (PM10. Phylogenetic reconstruction of haemagglutinin (HA gene from 168 influenza B viruses grouped them into Yamagata Clade 3 (65, 38.7%, Yamagata Clade 2 (48, 28.6% and Victoria Clade 1 (55, 32.7%. With neuraminidase (NA phylogeny, 30 intra-clade (29 within Yamagata Clade 3, 1 within Victoria Clade 1 and 1 inter-clade (Yamagata Clade 2-HA/Yamagata Clade 3-NA reassortants were identified. Study of virus temporal dynamics revealed a lineage shift from Victoria to Yamagata (2012-2013, and a clade shift from Yamagata Clade 2 to Clade 3 (2013-2014. Yamagata Clade 3 predominating in 2014 consisted of intra-clade reassortants that were closely related to a recent WHO vaccine candidate strain (B/Phuket/3073/2013, with the reassortment event occurred approximately 2 years ago based on Bayesian molecular clock estimation. Malaysian Victoria Clade 1 viruses carried H274Y substitution in the active site of neuraminidase, which confers resistance to oseltamivir. Statistical analyses on clinical and demographic data showed Yamagata-infected patients were older and more likely to experience headache while Victoria-infected patients were more likely to experience nasal congestion and
Colin J Worby
Full Text Available Genome sequencing is an increasingly common component of infectious disease outbreak investigations. However, the relationship between pathogen transmission and observed genetic data is complex, and dependent on several uncertain factors. As such, simulation of pathogen dynamics is an important tool for interpreting observed genomic data in an infectious disease outbreak setting, in order to test hypotheses and to explore the range of outcomes consistent with a given set of parameters. We introduce 'seedy', an R package for the simulation of evolutionary and epidemiological dynamics (http://cran.r-project.org/web/packages/seedy/. Our software implements stochastic models for the accumulation of mutations within hosts, as well as individual-level disease transmission. By allowing variables such as the transmission bottleneck size, within-host effective population size and population mixing rates to be specified by the user, our package offers a flexible framework to investigate evolutionary dynamics during disease outbreaks. Furthermore, our software provides theoretical pairwise genetic distance distributions to provide a likelihood of person-to-person transmission based on genomic observations, and using this framework, implements transmission route assessment for genomic data collected during an outbreak. Our open source software provides an accessible platform for users to explore pathogen evolution and outbreak dynamics via simulation, and offers tools to assess observed genomic data in this context.
TRAWA, LWR Dynamic by Coupled Neutron Diffusion and Thermohydraulics Calculation
1 - Description of problem or function: The purpose of the program is to study reactor dynamics in thermal water-cooled reactors. It treats the core as one or a few axially one-dimensional subregions. The two group neutron diffusion equations are solved simultaneously with the heat conduction equations and the two-phase hydraulic equations for one or more channels. Neither thermal nor hydraulic mixing appear between channels. Doppler, coolant density, coolant temperature, and soluble poison density feedbacks due to the thermo- hydraulics of the channels are described by using polynomial expansions for the group constants. The hydraulic circuit outside the reactor core consists of by-pass channels and risers with two- phase flow and of pump lines with incompressible flow. Various transients can be calculated by applying external disturbances. They can affect e.g. on movements of control rods, core inlet hydraulic conditions, system pressure or coefficients of neutronic shape function expansion between subregions. 2 - Method of solution: Nontrivial implicit methods are employed in the discretization of the equations to allow for sparse spatial mesh and flexible choice of time steps. The same spatial and temporal discretization is used for neutronics and thermohydraulics. 3 - Restrictions on the complexity of the problem: The dimensions of the program variable tables can easily be extended. Now the main dimensions are: 52 axial mesh points in core; 3 subregions; 10 axial regions with different fuel compositions; 7 radial mesh points in fuel rod; 6 delayed neutron groups; 6 coupled legs in pressure balance calculation; No flow reversals are allowed
Roaming dynamics in the MgH + H→Mg + H 2 reaction: Quantum dynamics calculations
Takayanagi, Toshiyuki; Tanaka, Tomokazu
2011-03-01
Reaction mechanisms of the MgH + H→Mg + H 2 reaction have been investigated using quantum reactive scattering methods on a global ab initio potential energy surface. There exist two microscopic mechanisms in the dynamics of this reaction. One is a direct hydrogen abstraction reaction and the other proceeds via initial formation of a HMgH complex in the deep potential well. The result of the present quantum dynamics calculations suggests that the HMgH complex formed in the reaction mainly decays into the Mg + H 2 channel via a 'roaming mechanism' without going through the saddle point region.
Self-consistent calculation of spin transport and magnetization dynamics
Lee, Kyung-Jin; Stiles, M. D.; Lee, Hyun-Woo; Moon, Jung-Hwan; Kim, Kyoung-Whan; Lee, Seo-Won
2013-01-01
A spin-polarized current transfers its spin-angular momentum to a local magnetization, exciting current-induced magnetization dynamics. So far, most studies in this field have focused on the direct effect of spin transport on magnetization dynamics, but ignored the feedback from the magnetization dynamics to the spin transport and back to the magnetization dynamics. Although the feedback is usually weak, there are situations when it can play an important role in the dynamics. In such situatio...
Neumann, Frank; Witt, Carsten
combinatorial optimization problem, namely makespan scheduling. We study the model of a strong adversary which is allowed to change one job at regular intervals. Furthermore, we investigate the setting of random changes. Our results show that randomized local search and a simple evolutionary algorithm are very...
Synthesizing mixed H2/H-infinity dynamic controller using evolutionary algorithms
Pedersen, Gerulf; Langballe, A.S.; Wisniewski, Rafal
2001-01-01
This paper covers the design of an Evolutionary Algorithm (EA), which should be able to synthesize a mixed H2/H-infinity. It will be shown how a system can be expressed as Matrix Inequalities (MI) and these will then be used in the design of the EA. The main objective is to examine whether a mixed...
An evolutionary Algorithm for Structural Subjected to Dynamic Loading With Random Excitation
This paper presents an evolutionary algorithm for optimization of structures subjected to random excitation. The new iteration scheme in conjunction with multi-population genetic strategy, entropy-based searching technique with narrowing down space and the quasi-exactness penalty function is developed to ensure rapid and steady convergence. Numerical example shows that proposed method has good accuracy and efficiency
Synthesizing multi-objective H2/H-infinity dynamic controller using evolutionary algorithms
Pedersen, Gerulf; Langballe, A.S.; Wisniewski, Rafal
This paper covers the design of an Evolutionary Algorithm (EA), which should be able to synthesize a mixed H2/H-infinity. It will be shown how a system can be expressed as Matrix Inequalities (MI) and these will then be used in the design of the EA. The main objective is to examine whether a mixed...
Zheng, Chaozhi; Ovaskainen, Otso; Hanski, Ilkka
2009-01-01
Dispersal comprises a complex life-history syndrome that influences the demographic dynamics of especially those species that live in fragmented landscapes, the structure of which may in turn be expected to impose selection on dispersal. We have constructed an individual-based evolutionary sexual model of dispersal for species occurring as metapopulations in habitat patch networks. The model assumes correlated random walk dispersal with edge-mediated behaviour (habitat selection) and spatially correlated stochastic local dynamics. The model is parametrized with extensive data for the Glanville fritillary butterfly. Based on empirical results for a single nucleotide polymorphism (SNP) in the phosphoglucose isomerase (Pgi) gene, we assume that dispersal rate in the landscape matrix, fecundity and survival are affected by a locus with two alleles, A and C, individuals with the C allele being more mobile. The model was successfully tested with two independent empirical datasets on spatial variation in Pgi allele frequency. First, at the level of local populations, the frequency of the C allele is the highest in newly established isolated populations and the lowest in old isolated populations. Second, at the level of sub-networks with dissimilar numbers and connectivities of patches, the frequency of C increases with decreasing network size and hence with decreasing average metapopulation size. The frequency of C is the highest in landscapes where local extinction risk is high and where there are abundant opportunities to establish new populations. Our results indicate that the strength of the coupling of the ecological and evolutionary dynamics depends on the spatial scale and is asymmetric, demographic dynamics having a greater immediate impact on genetic dynamics than vice versa. PMID:19414467
Evolutionary dynamics of human autoimmune disease genes and malfunctioned immunological genes
Podder Soumita
2012-01-01
Full Text Available Abstract Background One of the main issues of molecular evolution is to divulge the principles in dictating the evolutionary rate differences among various gene classes. Immunological genes have received considerable attention in evolutionary biology as candidates for local adaptation and for studying functionally important polymorphisms. The normal structure and function of immunological genes will be distorted when they experience mutations leading to immunological dysfunctions. Results Here, we examined the fundamental differences between the genes which on mutation give rise to autoimmune or other immune system related diseases and the immunological genes that do not cause any disease phenotypes. Although the disease genes examined are analogous to non-disease genes in product, expression, function, and pathway affiliation, a statistically significant decrease in evolutionary rate has been found in autoimmune disease genes relative to all other immune related diseases and non-disease genes. Possible ways of accumulation of mutation in the three steps of the central dogma (DNA-mRNA-Protein have been studied to trace the mutational effects predisposed to disease consequence and acquiring higher selection pressure. Principal Component Analysis and Multivariate Regression Analysis have established the predominant role of single nucleotide polymorphisms in guiding the evolutionary rate of immunological disease and non-disease genes followed by m-RNA abundance, paralogs number, fraction of phosphorylation residue, alternatively spliced exon, protein residue burial and protein disorder. Conclusions Our study provides an empirical insight into the etiology of autoimmune disease genes and other immunological diseases. The immediate utility of our study is to help in disease gene identification and may also help in medicinal improvement of immune related disease.
Jiaqin Shi
Full Text Available Despite their ubiquity and functional importance, microsatellites have been largely ignored in comparative genomics, mostly due to the lack of genomic information. In the current study, microsatellite distribution was characterized and compared in the whole genomes and both the coding and non-coding DNA sequences of the sequenced Brassica, Arabidopsis and other angiosperm species to investigate their evolutionary dynamics in plants. The variation in the microsatellite frequencies of these angiosperm species was much smaller than those for their microsatellite numbers and genome sizes, suggesting that microsatellite frequency may be relatively stable in plants. The microsatellite frequencies of these angiosperm species were significantly negatively correlated with both their genome sizes and transposable elements contents. The pattern of microsatellite distribution may differ according to the different genomic regions (such as coding and non-coding sequences. The observed differences in many important microsatellite characteristics (especially the distribution with respect to motif length, type and repeat number of these angiosperm species were generally accordant with their phylogenetic distance, which suggested that the evolutionary dynamics of microsatellite distribution may be generally consistent with plant divergence/evolution. Importantly, by comparing these microsatellite characteristics (especially the distribution with respect to motif type the angiosperm species (aside from a few species all clustered into two obviously different groups that were largely represented by monocots and dicots, suggesting a complex and generally dichotomous evolutionary pattern of microsatellite distribution in angiosperms. Polyploidy may lead to a slight increase in microsatellite frequency in the coding sequences and a significant decrease in microsatellite frequency in the whole genome/non-coding sequences, but have little effect on the microsatellite
Costa, Gideão Wagner Werneck Félix da; Cioffi, Marcelo de Bello; Bertollo, Luiz Antonio Carlos; Molina, Wagner Franco
2016-03-01
Lutjanidae is a family of primarily marine and carnivorous fishes distributed in the Atlantic, Indian, and Pacific oceans, with enormous economic and ecological importance. In order to better clarify the conservative chromosomal evolution of Lutjanidae, we analyzed the evolutionary dynamics of 5 repetitive DNA classes in 5 Lutjanus and in 1 Ocyurus species from the Western Atlantic. The ribosomal 18S sites were generally located in a single chromosome pair, except for L. jocu and L. alexandrei where they are found in 2 pairs. In turn, the 5S rDNA sites are unique, terminal and nonsyntenic with the 18S rDNA sites. In 3 species analyzed, H3 hisDNA genes were found in 1 chromosomal pair. However, while L. jocu presented 2 H3 sites, O. chrysurus showed a noteworthy dispersion of this gene in almost all chromosomes of the karyotype. Retrotransposons Rex1 and Rex3 do not exhibit any association with the explosive distribution of H3 sequences in O. chrysurus. The low compartmentalization of Rex elements, in addition to the general nondynamic distribution of ribosomal and H3 genes, corroborate the karyotype conservatism in Lutjanidae species, also at the microstructural level. However, some "disturbing evolutionary waves" can break down this conservative scenario, as evidenced by the massive random dispersion of H3 hisDNA in the genome of O. chrysurus. The implication of the genomic expansion of H3 histone genes and their functionality remain unknown, although suggesting that they have higher evolutionary dynamics than previously thought. PMID:26792596
Diversity, abundance, and evolutionary dynamics of Pong-like transposable elements in Triticeae.
Markova, Dragomira N; Mason-Gamer, Roberta J
2015-12-01
Pong-like elements are members of the PIF/Harbinger superfamily of DNA transposons that has been described in many plants, animals, and fungi. Most Pong elements contain two open reading frames (ORFs). One encodes a transposase (ORF2) that catalyzes transposition of Pong and related non-autonomous elements, while the function of the second is unknown. Little is known about the evolutionary history of Pong elements in flowering plants. In this work, we present the first comprehensive analysis of the diversity, abundance, and evolution of the Pong-like transposase gene in the genomes of 21 diploid species from the wheat tribe, Triticeae, and we present the first convincing evidence of horizontal transfer of nuclear-encoded Pong elements in any organism. A phylogenetic analysis of nearly 300 Pong sequences based on a conserved region of the transposase domain revealed a complex evolutionary history of Pong elements that can be best explained by ancestral polymorphism, followed by differential evolutionary success of some transposase lineages, and by occasional horizontal transfer between phylogenetically distant genera. In addition, we used transposon display to estimate the abundance of the transposase gene within Triticeae genomes, and our results revealed varying levels of Pong proliferation, with numbers of transposase copies ranging from 22 to 92. Comparisons of Pong transposase abundance to flow cytometry estimates of genome size revealed that larger Triticeae genome size was not correlated with transposase abundance. PMID:26206730
Dynamic model for calculating heating patterns during microwave sintering
In this paper the authors describe a thermal model that calculates the temperature distribution in ceramic samples and insulation under realistic sintering conditions. The calculation process involves a two-step procedure. The first step is to calculate the microwave power deposition in the sample and surrounding insulation. 3D FDTD calculations described in a companion paper are used for this purpose. The other step involves calculation of the temperature distribution using a 3D finite-difference heat-transfer program developed in our department. Results illustrating the effect of thickness of insulation and the placement of SiC rods in picket-fence arrangement are presented
Hansen, Sara Krogh; Vestergaard, Mikkel; Thøgersen, Lea; Schiøtt, Birgit; Nielsen, Niels Christian; Vosegaard, Thomas
2014-01-01
We present a method to calculate 31P solid-state NMR spectra based on the dynamic input from extended molecular dynamics (MD) simulations. The dynamic information confered by MD simulations is much more comprehensive than the information provided by traditional NMR dynamics models based on, for...... example, order parameters. Therefore, valuable insight into the dynamics of biomolecules may be achieved by the present method. We have applied this method to study the dynamics of lipid bilayers containing the antimicrobial peptide alamethicin, and we show that the calculated 31P spectra obtained with...
A comparison of techniques for calculating protein essential dynamics
van Aalten, D.M.F.; de Groot, B.L.; Findlay, J.B.C.; Berendsen, H.J.C.; Amadei, A
1997-01-01
Recently the basic theory of essential dynamics, a method for extracting large concerted motions from protein molecular dynamics trajectories, was described. Here, we introduce and test new aspects. A method for diagonalizing large covariance matrices is presented. We show that it is possible to per
Evolutionary Dynamics of Floral Homeotic Transcription Factor Protein–Protein Interactions
Bartlett, Madelaine; Thompson, Beth; Brabazon, Holly; Del Gizzi, Robert; Zhang, Thompson; Whipple, Clinton
2016-01-01
Protein–protein interactions (PPIs) have widely acknowledged roles in the regulation of development, but few studies have addressed the timing and mechanism of shifting PPIs over evolutionary history. The B-class MADS-box transcription factors, PISTILLATA (PI) and APETALA3 (AP3) are key regulators of floral development. PI-like (PIL) and AP3-like (AP3L) proteins from a number of plants, including Arabidopsis thaliana (Arabidopsis) and the grass Zea mays (maize), bind DNA as obligate heterodimers. However, a PIL protein from the grass relative Joinvillea can bind DNA as a homodimer. To ascertain whether Joinvillea PIL homodimerization is an anomaly or indicative of broader trends, we characterized PIL dimerization across the Poales and uncovered unexpected evolutionary lability. Both obligate B-class heterodimerization and PIL homodimerization have evolved multiple times in the order, by distinct molecular mechanisms. For example, obligate B-class heterodimerization in maize evolved very recently from PIL homodimerization. A single amino acid change, fixed during domestication, is sufficient to toggle one maize PIL protein between homodimerization and obligate heterodimerization. We detected a signature of positive selection acting on residues preferentially clustered in predicted sites of contact between MADS-box monomers and dimers, and in motifs that mediate MADS PPI specificity in Arabidopsis. Changing one positively selected residue can alter PIL dimerization activity. Furthermore, ectopic expression of a Joinvillea PIL homodimer in Arabidopsis can homeotically transform sepals into petals. Our results provide a window into the evolutionary remodeling of PPIs, and show that novel interactions have the potential to alter plant form in a context-dependent manner. PMID:26908583
Evolutionary Dynamics of Floral Homeotic Transcription Factor Protein-Protein Interactions.
Bartlett, Madelaine; Thompson, Beth; Brabazon, Holly; Del Gizzi, Robert; Zhang, Thompson; Whipple, Clinton
2016-06-01
Protein-protein interactions (PPIs) have widely acknowledged roles in the regulation of development, but few studies have addressed the timing and mechanism of shifting PPIs over evolutionary history. The B-class MADS-box transcription factors, PISTILLATA (PI) and APETALA3 (AP3) are key regulators of floral development. PI-like (PI(L)) and AP3-like (AP3(L)) proteins from a number of plants, including Arabidopsis thaliana (Arabidopsis) and the grass Zea mays (maize), bind DNA as obligate heterodimers. However, a PI(L) protein from the grass relative Joinvillea can bind DNA as a homodimer. To ascertain whether Joinvillea PI(L) homodimerization is an anomaly or indicative of broader trends, we characterized PI(L) dimerization across the Poales and uncovered unexpected evolutionary lability. Both obligate B-class heterodimerization and PI(L) homodimerization have evolved multiple times in the order, by distinct molecular mechanisms. For example, obligate B-class heterodimerization in maize evolved very recently from PI(L) homodimerization. A single amino acid change, fixed during domestication, is sufficient to toggle one maize PI(L) protein between homodimerization and obligate heterodimerization. We detected a signature of positive selection acting on residues preferentially clustered in predicted sites of contact between MADS-box monomers and dimers, and in motifs that mediate MADS PPI specificity in Arabidopsis. Changing one positively selected residue can alter PI(L) dimerization activity. Furthermore, ectopic expression of a Joinvillea PI(L) homodimer in Arabidopsis can homeotically transform sepals into petals. Our results provide a window into the evolutionary remodeling of PPIs, and show that novel interactions have the potential to alter plant form in a context-dependent manner. PMID:26908583
Two sexes, one genome: the evolutionary dynamics of intralocus sexual conflict
Pennell, Tanya M.; Morrow, Edward H.
2013-01-01
As the evolutionary interests of males and females are frequently divergent, a trait value that is optimal for the fitness of one sex is often not optimal for the other. A shared genome also means that the same genes may underlie the same trait in both sexes. This can give rise to a form of sexual antagonism, known as intralocus sexual conflict (IASC). Here, a tug-of-war over allelic expression can occur, preventing the sexes from reaching optimal trait values, thereby causing sex-specific re...
Tassos Patokos
2014-01-01
Full Text Available The paper presents an evolutionary model, based on the assumption that agents may revise their current strategies if they previously failed to attain the maximum level of potential payoffs. We offer three versions of this reflexive mechanism, each one of which describes a distinct type: spontaneous agents, rigid players, and ‘satisficers’. We use simulations to examine the performance of these types. Agents who change their strategies relatively easily tend to perform better in coordination games, but antagonistic games generally lead to more favorable outcomes if the individuals only change their strategies when disappointment from previous rounds surpasses some predefined threshold.
Fabre, Frédéric; Montarry, Josselin; Coville, Jérôme; Senoussi, Rachid; Simon, Vincent; Moury, Benoît
2012-01-01
Uncovering how natural selection and genetic drift shape the evolutionary dynamics of virus populations within their hosts can pave the way to a better understanding of virus emergence. Mathematical models already play a leading role in these studies and are intended to predict future emergences. Here, using high-throughput sequencing, we analyzed the within-host population dynamics of four Potato virus Y (PVY) variants differing at most by two substitutions involved in pathogenicity properties. Model selection procedures were used to compare experimental results to six hypotheses regarding competitiveness and intensity of genetic drift experienced by viruses during host plant colonization. Results indicated that the frequencies of variants were well described using Lotka-Volterra models where the competition coefficients β(ij) exerted by variant j on variant i are equal to their fitness ratio, r(j)/r(i). Statistical inference allowed the estimation of the effect of each mutation on fitness, revealing slight (s = -0.45%) and high (s = -13.2%) fitness costs and a negative epistasis between them. Results also indicated that only 1 to 4 infectious units initiated the population of one apical leaf. The between-host variances of the variant frequencies were described using Dirichlet-multinomial distributions whose scale parameters, closely related to the fixation index F(ST), were shown to vary with time. The genetic differentiation of virus populations among plants increased from 0 to 10 days post-inoculation and then decreased until 35 days. Overall, this study showed that mathematical models can accurately describe both selection and genetic drift processes shaping the evolutionary dynamics of viruses within their hosts. PMID:22532800
Fe IX Calculations for the Solar Dynamics Observatory
Foster, Adam R
2011-01-01
New calculations of the energy levels, radiative transition rates and collisional excitation rates of \\ion{Fe}{ix} have been carried out using the Flexible Atomic Code, paying close attention to experimentally identified levels and extending existing calculations to higher energy levels. For lower levels, R-matrix collisional excitation rates from earlier work have been used. Significant emission is predicted by these calculations in the 5f-3d transitions, which will impact analysis of SDO AIA observations using the 94\\AA\\ filter.
Dynamic Calculation Method of Beam System Under Low Velocity Impact
LI Wen-pei; WANG De-rong; SONG Chun-ming; WANG Ming-yang
2008-01-01
The beating beams and the supporting beams under low velocity impact may be in four different strain stages of deformation depending on the impact intensity and beam structure strength. Based on the different judging conditions of deformation stages, the corresponding calculation models are proposed, the calculation formulae for the determination of the impact force and the beam's lateral displacement are obtained. Calculation shows that the beam's total deflection is small when the flexibility of the supporting component is high and the effect of diminishing deflection disappears almost when the stiffness of the supporting component is high.
Hapuarachchi, Hapuarachchige Chanditha; Koo, Carmen; Kek, Relus; Xu, Helen; Lai, Yee Ling; Liu, Lilac; Kok, Suet Yheng; Shi, Yuan; Chuen, Raphael Lee Tze; Lee, Kim-Sung; Maurer-Stroh, Sebastian; Ng, Lee Ching
2016-01-01
Dengue virus (DENV) is currently the most prevalent mosquito-borne viral pathogen. DENVs naturally exist as highly heterogeneous populations. Even though the descriptions on DENV diversity are plentiful, only a few studies have narrated the dynamics of intra-epidemic virus diversity at a fine scale. Such accounts are important to decipher the reciprocal relationship between viral evolutionary dynamics and disease transmission that shape dengue epidemiology. In the current study, we present a micro-scale genetic analysis of a monophyletic lineage of DENV-1 genotype III (epidemic lineage) detected from November 2012 to May 2014. The lineage was involved in an unprecedented dengue epidemic in Singapore during 2013-2014. Our findings showed that the epidemic lineage was an ensemble of mutants (variants) originated from an initial mixed viral population. The composition of mutant spectrum was dynamic and positively correlated with case load. The close interaction between viral evolution and transmission intensity indicated that tracking genetic diversity through time is potentially a useful tool to infer DENV transmission dynamics and thereby, to assess the epidemic risk in a disease control perspective. Moreover, such information is salient to understand the viral basis of clinical outcome and immune response variations that is imperative to effective vaccine design. PMID:26940650
Jing Chen
2013-01-01
Full Text Available Due to high efficiency and good scalability, hierarchical hybrid P2P architecture has drawn more and more attention in P2P streaming research and application fields recently. The problem about super peer selection, which is the key problem in hybrid heterogeneous P2P architecture, is becoming highly challenging because super peers must be selected from a huge and dynamically changing network. A distributed super peer selection (SPS algorithm for hybrid heterogeneous P2P streaming system based on evolutionary game is proposed in this paper. The super peer selection procedure is modeled based on evolutionary game framework firstly, and its evolutionarily stable strategies are analyzed. Then a distributed Q-learning algorithm (ESS-SPS according to the mixed strategies by analysis is proposed for the peers to converge to the ESSs based on its own payoff history. Compared to the traditional randomly super peer selection scheme, experiments results show that the proposed ESS-SPS algorithm achieves better performance in terms of social welfare and average upload rate of super peers and keeps the upload capacity of the P2P streaming system increasing steadily with the number of peers increasing.
Atienza, David; Baloukas, Christos; Papadopoulos, Lazaros; Poucet, Christophe; Mamagkakis, Stylianos; Hidalgo, Jose I.; Catthoor, Francky; Soudris, Dimitrios; Lanchares, Juan
2007-01-01
Embedded consumer devices are increasing their capabilities and can now implement new multimedia applications reserved only for powerful desktops a few years ago. These applications share complex and intensive dynamic memory use. Thus, dynamic memory optimizations are a requirement when porting these applications. Within these optimizations, the refinement of the Dynamically (de)allocated Data Type (or DDT) implementations is one of the most important and difficult parts for an efficient mapp...
Between Pleasure and Contentment: Evolutionary Dynamics of Some Possible Parameters of Happiness
Gao, Yue; Edelman, Shimon
2016-01-01
We offer and test a simple operationalization of hedonic and eudaimonic well-being (“happiness”) as mediating variables that link outcomes to motivation. In six evolutionary agent-based simulation experiments, we compared the relative performance of agents endowed with different combinations of happiness-related traits (parameter values), under four types of environmental conditions. We found (i) that the effects of attaching more weight to longer-term than to momentary happiness and of extending the memory for past happiness are both stronger in an environment where food is scarce; (ii) that in such an environment “relative consumption,” in which the agent’s well-being is negatively affected by that of its neighbors, is more detrimental to survival when food is scarce; and (iii) that having a positive outlook, under which agents’ longer-term happiness is increased by positive events more than it is decreased by negative ones, is generally advantageous. PMID:27144982
The evolutionary dynamics of a population model with a strong Allee effect.
Cushing, Jim M
2015-08-01
An evolutionary game theoretic model for a population subject to predation and a strong Allee threshold of extinction is analyzed using, among other methods, Poincaré-Bendixson theory. The model is a nonlinear, plane autonomous system whose state variables are population density and the mean of a phenotypic trait, which is subject to Darwinian evolution, that determines the population's inherent (low density) growth rate (fitness). A trade-off is assumed in that an increase in the inherent growth rate results in a proportional increase in the predator's attack rate. The main results are that orbits equilibrate (there are no cycles or cycle chains of saddles), that the extinction set (or Allee basin) shrinks when evolution occurs, and that the meant trait component of survival equilibria occur at maxima of the inherent growth rate (as a function of the trait). PMID:25974340
Between Pleasure and Contentment: Evolutionary Dynamics of Some Possible Parameters of Happiness.
Gao, Yue; Edelman, Shimon
2016-01-01
We offer and test a simple operationalization of hedonic and eudaimonic well-being ("happiness") as mediating variables that link outcomes to motivation. In six evolutionary agent-based simulation experiments, we compared the relative performance of agents endowed with different combinations of happiness-related traits (parameter values), under four types of environmental conditions. We found (i) that the effects of attaching more weight to longer-term than to momentary happiness and of extending the memory for past happiness are both stronger in an environment where food is scarce; (ii) that in such an environment "relative consumption," in which the agent's well-being is negatively affected by that of its neighbors, is more detrimental to survival when food is scarce; and (iii) that having a positive outlook, under which agents' longer-term happiness is increased by positive events more than it is decreased by negative ones, is generally advantageous. PMID:27144982
Dynamic instability of cooperation due to diverse activity patterns in evolutionary social dilemmas
Xia, Cheng-Yi; Perc, Matjaz; Moreno, Yamir
2015-01-01
Individuals might abstain from participating in an instance of an evolutionary game for various reasons, ranging from lack of interest to risk aversion. In order to understand the consequences of such diverse activity patterns on the evolution of cooperation, we study a weak prisoner's dilemma where each player's participation is probabilistic rather than certain. Players that do not participate get a null payoff and are unable to replicate. We show that inactivity introduces cascading failures of cooperation, which are particularly severe on scale-free networks with frequently inactive hubs. The drops in the fraction of cooperators are sudden, while the spatiotemporal reorganization of compact cooperative clusters, and thus the recovery, takes time. Nevertheless, if the activity of players is directly proportional to their degree, or if the interaction network is not strongly heterogeneous, the overall evolution of cooperation is not impaired. This is because inactivity negatively affects the potency of low-...
Dynamic simulation of flash drums using rigorous physical property calculations
F. M. Gonçalves
2007-06-01
Full Text Available The dynamics of flash drums is simulated using a formulation adequate for phase modeling with equations of state (EOS. The energy and mass balances are written as differential equations for the internal energy and the number of moles of each species. The algebraic equations of the model, solved at each time step, are those of a flash with specified internal energy, volume and mole numbers (UVN flash. A new aspect of our dynamic simulations is the use of direct iterations in phase volumes (instead of pressure for solving the algebraic equations. It was also found that an iterative procedure previously suggested in the literature for UVN flashes becomes unreliable close to phase boundaries and a new alternative is proposed. Another unusual aspect of this work is that the model expressions, including the physical properties and their analytical derivatives, were quickly implemented using computer algebra.
A dynamical collective calculation of supernova neutrino signals
Gava, J.; Kneller, J.; Volpe, C.; McLaughlin, G.C.
2009-01-01
We present the first calculations with three flavors of collective and shock wave effects for neutrino propagation in core-collapse supernovae using hydroynamical density profiles and the S matrix formalism. We explore the interplay between the neutrino-neutrino interaction and the effects of multiple resonances upon the time signal of positrons in supernova observatories. A specific signature is found for the inverted hierarchy and a large third neutrino mixing angle and we predict, in this ...
Development of a dynamic calculation tool forsimulation of ditching
Pilorget, Marc
2011-01-01
The present document is the final master thesis report written by Marc PILORGET,student at SUPAERO (home institution) and KTH (Royal Institute of Technology,Exchange University). This six months internship was done at DASSAULT AVIATION(Airframe engineering department) based in Saint-Cloud, France. It spanned from the 5thof July to the 23rd of December. The thesis work aims at developing an SPH (SmoothParticle Hydrodynamics) calculation method for ditching and implementing it in the finiteelem...
Hydraulic calculation and dynamic analysis of columnar reversing gate
You-an HU
2011-09-01
Full Text Available According to the hydraulic calculation principles of the orifice outflow, the discharge capacity of the columnar reversing gate under the partial opening condition was calculated and checked. Using ANSYS, a large finite element analysis software, the discharge process was simulated. The distribution rule of the velocities in the gate chamber and downstream channel was obtained. An FEM model of the columnar reversing gate was built, and the natural vibration properties of the gate were analyzed. Based on the Westergaard added mass method, the added mass caused by the fluid-structure coupling motion was taken into account, and the effects of the coupling interaction were discussed. The results show that the size of the small gates meets the demand for discharge capacity, the current in the gate chamber is quite turbulent, the trunnion and arms are obviously impacted by flow, and the effects of water on vibration characteristics are remarkable. The study provides a reference for the design and calculation of gates of the same type.
In the present paper we report a theoretical study on mechanistic photodissociation of acetaldehyde (CH3CHO). Stationary structures for H2 and CO eliminations in the ground state (S0) have been optimized with density functional theory method, which is followed by the intrinsic reaction coordinate and ab initio molecular dynamics calculations to confirm the elimination mechanism. Equilibrium geometries, transition states, and intersection structures for the C-C and C-H dissociations in excited states were determined by the complete-active-space self-consistent field (CASSCF) method. Based on the CASSCF optimized structures, the potential energy profiles for the dissociations were refined by performing the single-point calculations using the multireference configuration interaction method. Upon the low-energy irradiation of CH3CHO (265 nm1 C-C bond fission following intersystem crossing from the S1 state is the predominant channel and the minor channel, the ground-state elimination to CH4+CO after internal conversion (IC) from S1 to S0, could not be excluded. With the photon energy increasing, another pathway of IC, achieved via an S1/S0 intersection point resulting from the S1 C-C bond fission, becomes accessible and increases the yield of CH4+CO.
A Comprehensive Curation Shows the Dynamic Evolutionary Patterns of Prokaryotic CRISPRs.
Mai, Guoqin; Ge, Ruiquan; Sun, Guoquan; Meng, Qinghan; Zhou, Fengfeng
2016-01-01
Motivation. Clustered regularly interspaced short palindromic repeat (CRISPR) is a genetic element with active regulation roles for foreign invasive genes in the prokaryotic genomes and has been engineered to work with the CRISPR-associated sequence (Cas) gene Cas9 as one of the modern genome editing technologies. Due to inconsistent definitions, the existing CRISPR detection programs seem to have missed some weak CRISPR signals. Results. This study manually curates all the currently annotated CRISPR elements in the prokaryotic genomes and proposes 95 updates to the annotations. A new definition is proposed to cover all the CRISPRs. The comprehensive comparison of CRISPR numbers on the taxonomic levels of both domains and genus shows high variations for closely related species even in the same genus. The detailed investigation of how CRISPRs are evolutionarily manipulated in the 8 completely sequenced species in the genus Thermoanaerobacter demonstrates that transposons act as a frequent tool for splitting long CRISPRs into shorter ones along a long evolutionary history. PMID:27195295
A Comprehensive Curation Shows the Dynamic Evolutionary Patterns of Prokaryotic CRISPRs
Guoqin Mai
2016-01-01
Full Text Available Motivation. Clustered regularly interspaced short palindromic repeat (CRISPR is a genetic element with active regulation roles for foreign invasive genes in the prokaryotic genomes and has been engineered to work with the CRISPR-associated sequence (Cas gene Cas9 as one of the modern genome editing technologies. Due to inconsistent definitions, the existing CRISPR detection programs seem to have missed some weak CRISPR signals. Results. This study manually curates all the currently annotated CRISPR elements in the prokaryotic genomes and proposes 95 updates to the annotations. A new definition is proposed to cover all the CRISPRs. The comprehensive comparison of CRISPR numbers on the taxonomic levels of both domains and genus shows high variations for closely related species even in the same genus. The detailed investigation of how CRISPRs are evolutionarily manipulated in the 8 completely sequenced species in the genus Thermoanaerobacter demonstrates that transposons act as a frequent tool for splitting long CRISPRs into shorter ones along a long evolutionary history.
Evolutionary history of chordate PAX genes: dynamics of change in a complex gene family.
Vanessa Rodrigues Paixão-Côrtes
Full Text Available Paired box (PAX genes are transcription factors that play important roles in embryonic development. Although the PAX gene family occurs in animals only, it is widely distributed. Among the vertebrates, its 9 genes appear to be the product of complete duplication of an original set of 4 genes, followed by an additional partial duplication. Although some studies of PAX genes have been conducted, no comprehensive survey of these genes across the entire taxonomic unit has yet been attempted. In this study, we conducted a detailed comparison of PAX sequences from 188 chordates, which revealed restricted variation. The absence of PAX4 and PAX8 among some species of reptiles and birds was notable; however, all 9 genes were present in all 74 mammalian genomes investigated. A search for signatures of selection indicated that all genes are subject to purifying selection, with a possible constraint relaxation in PAX4, PAX7, and PAX8. This result indicates asymmetric evolution of PAX family genes, which can be associated with the emergence of adaptive novelties in the chordate evolutionary trajectory.
Calculated dynamical scattering in a fibrous protein structure
The effects of dynamical electron scattering on the amplitudes and phases diffracted by the microcrystals of an alpha-helical coiled-coil protein have been investigated by computer simulation using a model structure derived from relatively low resolution data recorded at liquid nitrogen temperature. The projected potential of this model show only small departure from the weak-phase-object approximation predictions, with amplitudes increasing almost linearly with thickness and phases nearly constant up to thickness of 200 Angstroems. 4 refs., 3 figs
Non-linear calculation of PCRV using dynamic relaxation
A brief review is presented of a numerical method called the dynamic relaxation method for stress analysis of the concrete in prestressed concrete pressure vessels. By this method the three-dimensional elliptic differential equations of the continuum are changed into the four-dimensional hyperbolic differential equations known as wave equations. The boundary value problem of the static system is changed into an initial and boundary value problem for which a solution exists if the physical system is defined at time t=0. The effect of non-linear stress-strain behaviour of the material as well as creep and cracking are considered
Parallel beam dynamics calculations on high performance computers
Faced with a backlog of nuclear waste and weapons plutonium, as well as an ever-increasing public concern about safety and environmental issues associated with conventional nuclear reactors, many countries are studying new, accelerator-driven technologies that hold the promise of providing safe and effective solutions to these problems. Proposed projects include accelerator transmutation of waste (ATW), accelerator-based conversion of plutonium (ABC), accelerator-driven energy production (ADEP), and accelerator production of tritium (APT). Also, next-generation spallation neutron sources based on similar technology will play a major role in materials science and biological science research. The design of accelerators for these projects will require a major advance in numerical modeling capability. For example, beam dynamics simulations with approximately 100 million particles will be needed to ensure that extremely stringent beam loss requirements (less than a nanoampere per meter) can be met. Compared with typical present-day modeling using 10,000-100,000 particles, this represents an increase of 3-4 orders of magnitude. High performance computing (HPC) platforms make it possible to perform such large scale simulations, which require 10's of GBytes of memory. They also make it possible to perform smaller simulations in a matter of hours that would require months to run on a single processor workstation. This paper will describe how HPC platforms can be used to perform the numerically intensive beam dynamics simulations required for development of these new accelerator-driven technologies
Sumathy, K; Ella, Krishna M
2012-03-01
The genetic diversity of Chikungunya virus (CHIKV) causing recurring outbreaks in India since 2006 was studied. The 2006 epidemic was caused by a virus strain of the East, Central and South African (ECSA) genotype with 226A in the E1 glycoprotein. The variant strain with E1-A226V mutation caused outbreaks since 2007 in the state of Kerala where Aedes albopictus is the abundant mosquito vector. Molecular epidemiology data since 2007 is scarce from other regions of the country. RT-PCR, sequencing and phylogenetic analyses of CHIKV isolates from the 2009 to 2010 epidemics in the States of Tamil Nadu and Andhra Pradesh placed them in a separate clade within the ECSA lineage. The isolates of the study had 226A in the E1 glycoprotein. The isolates had a novel E1-K211E mutation that was under significant positive selection. E1-211E is highly conserved in the Asian genotype of the virus circulated by Aedes aegypti. Unique mutations in E2 glycoprotein were identified. The two sub-lineages of ECSA genotype circulating in India parallel the abundance of Ae. albopictus and Ae. aegypti. Novel mutations in the envelope glycoproteins suggest adaptive evolution of the virus to local vector abundance. Cross neutralization of the virus isolates from recurring Indian epidemics indicated that no distinct serotypes had evolved. The study has provided insights into the origin, distribution and evolutionary adaptation of the virus to local vector abundance in the region that has reportedly, the highest incidence of CHIKV infection in the world. PMID:22246833
Toward a Mechanics of Adaptive Behavior: Evolutionary Dynamics and Matching Theory Statics
McDowell, J. J.; Popa, Andrei
2010-01-01
One theory of behavior dynamics instantiates the idea that behavior evolves in response to selection pressure from the environment in the form of reinforcement. This computational theory implements Darwinian principles of selection, reproduction, and mutation, which operate on a population of potential behaviors by means of a genetic algorithm.…
Ordway, Stephen; King, Dawn; Bahar, Sonya
Reaction-diffusion processes, such as branching-coalescing random walks, can be used to describe the underlying dynamics of nonequilibrium phase transitions. In an agent-based, neutral model of evolutionary dynamics, we have previously shown that our system undergoes a continuous, nonequilibrium phase transition, from extinction to survival, as various system parameters were tuned. This model was shown to belong to the directed percolation (DP) universality class, by measuring the critical exponents corresponding to correlation length ξ⊥, correlation time ξ| |, and particle density β. The fourth critical exponent that defines the DP universality class is β', which measures the survival probability of growth from a single seed organism. Since DP universality is theorized to have time-reversal symmetry, it is assumed that β = β '. In order to confirm the existence of time-reversal symmetry in our model, we evaluate the system growth from a single asexually reproducing organism. Importantly, the critical exponent β' could be useful for comparison to experimental studies of phase transitions in biological systems, since observing growth of microbial populations is significantly easier than observing death. This research was supported by funding from the James S. McDonnell Foundation.
From Binding-Induced Dynamic Effects in SH3 Structures to Evolutionary Conserved Sectors.
Ana Zafra Ruano
2016-05-01
Full Text Available Src Homology 3 domains are ubiquitous small interaction modules known to act as docking sites and regulatory elements in a wide range of proteins. Prior experimental NMR work on the SH3 domain of Src showed that ligand binding induces long-range dynamic changes consistent with an induced fit mechanism. The identification of the residues that participate in this mechanism produces a chart that allows for the exploration of the regulatory role of such domains in the activity of the encompassing protein. Here we show that a computational approach focusing on the changes in side chain dynamics through ligand binding identifies equivalent long-range effects in the Src SH3 domain. Mutation of a subset of the predicted residues elicits long-range effects on the binding energetics, emphasizing the relevance of these positions in the definition of intramolecular cooperative networks of signal transduction in this domain. We find further support for this mechanism through the analysis of seven other publically available SH3 domain structures of which the sequences represent diverse SH3 classes. By comparing the eight predictions, we find that, in addition to a dynamic pathway that is relatively conserved throughout all SH3 domains, there are dynamic aspects specific to each domain and homologous subgroups. Our work shows for the first time from a structural perspective, which transduction mechanisms are common between a subset of closely related and distal SH3 domains, while at the same time highlighting the differences in signal transduction that make each family member unique. These results resolve the missing link between structural predictions of dynamic changes and the domain sectors recently identified for SH3 domains through sequence analysis.
From Binding-Induced Dynamic Effects in SH3 Structures to Evolutionary Conserved Sectors.
Zafra Ruano, Ana; Cilia, Elisa; Couceiro, José R; Ruiz Sanz, Javier; Schymkowitz, Joost; Rousseau, Frederic; Luque, Irene; Lenaerts, Tom
2016-05-01
Src Homology 3 domains are ubiquitous small interaction modules known to act as docking sites and regulatory elements in a wide range of proteins. Prior experimental NMR work on the SH3 domain of Src showed that ligand binding induces long-range dynamic changes consistent with an induced fit mechanism. The identification of the residues that participate in this mechanism produces a chart that allows for the exploration of the regulatory role of such domains in the activity of the encompassing protein. Here we show that a computational approach focusing on the changes in side chain dynamics through ligand binding identifies equivalent long-range effects in the Src SH3 domain. Mutation of a subset of the predicted residues elicits long-range effects on the binding energetics, emphasizing the relevance of these positions in the definition of intramolecular cooperative networks of signal transduction in this domain. We find further support for this mechanism through the analysis of seven other publically available SH3 domain structures of which the sequences represent diverse SH3 classes. By comparing the eight predictions, we find that, in addition to a dynamic pathway that is relatively conserved throughout all SH3 domains, there are dynamic aspects specific to each domain and homologous subgroups. Our work shows for the first time from a structural perspective, which transduction mechanisms are common between a subset of closely related and distal SH3 domains, while at the same time highlighting the differences in signal transduction that make each family member unique. These results resolve the missing link between structural predictions of dynamic changes and the domain sectors recently identified for SH3 domains through sequence analysis. PMID:27213566
From Binding-Induced Dynamic Effects in SH3 Structures to Evolutionary Conserved Sectors
Ruiz Sanz, Javier; Schymkowitz, Joost; Rousseau, Frederic
2016-01-01
Src Homology 3 domains are ubiquitous small interaction modules known to act as docking sites and regulatory elements in a wide range of proteins. Prior experimental NMR work on the SH3 domain of Src showed that ligand binding induces long-range dynamic changes consistent with an induced fit mechanism. The identification of the residues that participate in this mechanism produces a chart that allows for the exploration of the regulatory role of such domains in the activity of the encompassing protein. Here we show that a computational approach focusing on the changes in side chain dynamics through ligand binding identifies equivalent long-range effects in the Src SH3 domain. Mutation of a subset of the predicted residues elicits long-range effects on the binding energetics, emphasizing the relevance of these positions in the definition of intramolecular cooperative networks of signal transduction in this domain. We find further support for this mechanism through the analysis of seven other publically available SH3 domain structures of which the sequences represent diverse SH3 classes. By comparing the eight predictions, we find that, in addition to a dynamic pathway that is relatively conserved throughout all SH3 domains, there are dynamic aspects specific to each domain and homologous subgroups. Our work shows for the first time from a structural perspective, which transduction mechanisms are common between a subset of closely related and distal SH3 domains, while at the same time highlighting the differences in signal transduction that make each family member unique. These results resolve the missing link between structural predictions of dynamic changes and the domain sectors recently identified for SH3 domains through sequence analysis. PMID:27213566
Phase diagram of kaolinite from Molecular Dynamics calculations
Structural and thermal behaviors of uncharged 1:1 phyllosilicates kaolinite were investigated from molecular dynamics simulations based on the CLAYFF force field. The focus is on the variation of structural properties including bulk modulus with pressure from 0 to 20 GPa under various range of temperature. The largest bulk modulus between the pressures of 200 and 800 MPa varies from 80 GPa at 298 K to 50 GPa at 1473 K. The obtained value of Cp varies between 7.8 and 13.6 Kcal mol-1 K-1 in the pressure range of 0.1 MPa-20 GPa. Besides, a huge difference was noticed regarding the computed properties at the superheating point. Finally, we show the relationship between superheating point temperature and pressure leading to a phase diagram of kaolinite.
Parallel beam dynamics calculations on high performance computers
Faced with a backlog of nuclear waste and weapons plutonium, as well as an ever-increasing public concern about safety and environmental issues associated with conventional nuclear reactors, many countries are studying new, accelerator-driven technologies that hold the promise of providing safe and effective solutions to these problems. Proposed projects include accelerator transmutation of waste (ATW), accelerator-based conversion of plutonium (ABC), accelerator-driven energy production (ADEP), and accelerator production of tritium (APT). Also, next-generation spallation neutron sources based on similar technology will play a major role in materials science and biological science research. The design of accelerators for these projects will require a major advance in numerical modeling capability. For example, beam dynamics simulations with approximately 100 million particles will be needed to ensure that extremely stringent beam loss requirements (less than a nanoampere per meter) can be met. Compared with typical present-day modeling using 10,000 endash 100,000 particles, this represents an increase of 3 endash 4 orders of magnitude. High performance computing (HPC) platforms make it possible to perform such large scale simulations, which require 10 close-quote s of GBytes of memory. They also make it possible to perform smaller simulations in a matter of hours that would require months to run on a single processor workstation. This paper will describe how HPC platforms can be used to perform the numerically intensive beam dynamics simulations required for development of these new accelerator-driven technologies. copyright 1997 American Institute of Physics
Body drop into a fluid tank and dynamic loads calculation
Komarov Aleksandr Andreevich
2014-05-01
Full Text Available The theory of a body striking a fluid began intensively developing due to the tasks of hydroplanes landing. For the recent years the study of a stroke and submersion of bodies into fluid became even more current. We face them in the process of strength calculation of ship hulls and other structures in modern technology. These tasks solution represents great mathematical difficulty even in case of the mentioned simplifications. These difficulties emerge due to the unsteady character of fluid motion in case of body submersion, and also jet and spray phenomena, which lead to discontinuous motions. On the basis of G.V. Logvinovich’s concept the problem of loads determination with consideration for air gap is solved for both a body and reservoir enclosing structures when a body falls into a fluid. Numerical method is based on the decay of an arbitrary discontinuity.
Calculating Free Energies Using Scaled-Force Molecular Dynamics Algorithm
Darve, Eric; Wilson, Micahel A.; Pohorille, Andrew
2000-01-01
One common objective of molecular simulations in chemistry and biology is to calculate the free energy difference between different states of the system of interest. Examples of problems that have such an objective are calculations of receptor-ligand or protein-drug interactions, associations of molecules in response to hydrophobic, and electrostatic interactions or partition of molecules between immiscible liquids. Another common objective is to describe evolution of the system towards a low energy (possibly the global minimum energy), 'native' state. Perhaps the best example of such a problem is folding of proteins or short RNA molecules. Both types of problems share the same difficulty. Often, different states of the system are separated by high energy barriers, which implies that transitions between these states are rare events. This, in turn, can greatly impede exploration of phase space. In some instances this can lead to 'quasi non-ergodicity', whereby a part of phase space is inaccessible on timescales of the simulation. A host of strategies has been developed to improve efficiency of sampling the phase space. For example, some Monte Carlo techniques involve large steps which move the system between low-energy regions in phase space without the need for sampling the configurations corresponding to energy barriers (J-walking). Most strategies, however, rely on modifying probabilities of sampling low and high-energy regions in phase space such that transitions between states of interest are encouraged. Perhaps the simplest implementation of this strategy is to increase the temperature of the system. This approach was successfully used to identify denaturation pathways in several proteins, but it is clearly not applicable to protein folding. It is also not a successful method for determining free energy differences. Finally, the approach is likely to fail for systems with co-existing phases, such as water-membrane systems, because it may lead to spontaneous
Sorokin, Valery A.; Gelfand, Mikhail S.; Artamonova, Irena I.
2010-01-01
Clustered regularly interspaced short palindromic repeats (CRISPRs) form a recently characterized type of prokaryotic antiphage defense system. The phage-host interactions involving CRISPRs have been studied in experiments with selected bacterial or archaeal species and, computationally, in completely sequenced genomes. However, these studies do not allow one to take prokaryotic population diversity and phage-host interaction dynamics into account. This gap can be filled by using metagenomic ...
Lande, Russell; Engen, Steinar; Sæther, Bernt-Erik
2009-01-01
The evolution of population dynamics in a stochastic environment is analysed under a general form of density-dependence with genetic variation in r and K, the intrinsic rate of increase and carrying capacity in the average environment, and in σe2, the environmental variance of population growth rate. The continuous-time model assumes a large population size and a stationary distribution of environments with no autocorrelation. For a given population density, N, and genotype frequency, p, the ...
Sexual selection and the evolutionary dynamics of the major histocompatibility complex
Jan Ejsmond, Maciej; Radwan, Jacek; Wilson, Anthony B
2014-01-01
The genes of the major histocompatibility complex (MHC) are a key component of the adaptive immune system and among the most variable loci in the vertebrate genome. Pathogen-mediated natural selection and MHC-based disassortative mating are both thought to structure MHC polymorphism, but their effects have proven difficult to discriminate in natural systems. Using the first model of MHC dynamics incorporating both survival and reproduction, we demonstrate that natural and sexual selection pro...
Co-Evolutionary Mechanisms of Emotional Bursts in Online Social Dynamics and Networks
Bosiljka Tadić
2013-11-01
Full Text Available Collective emotional behavior of users is frequently observed on various Web portals; however, its complexity and the role of emotions in the acting mechanisms are still not thoroughly understood. In this work, using the empirical data and agent-based modeling, a parallel analysis is performed of two archetypal systems—Blogs and Internet-Relayed-Chats—both of which maintain self-organized dynamics but not the same communication rules and time scales. The emphasis is on quantifying the collective emotions by means of fractal analysis of the underlying processes as well as topology of social networks, which arise and co-evolve in these stochastic processes. The results reveal that two distinct mechanisms, which are based on different use of emotions (an emotion is characterized by two components, arousal and valence, are intrinsically associated with two classes of emergent social graphs. Their hallmarks are the evolution of communities in accordance with the excess of the negative emotions on popular Blogs, on one side, and smooth spreading of the Bot’s emotional impact over the entire hierarchical network of chats, on the other. Another emphasis of this work is on the understanding of nonextensivity of the emotion dynamics; it was found that, in its own way, each mechanism leads to a reduced phase space of the emotion components when the collective dynamics takes place. That a non-additive entropy describes emotion dynamics, is further confirmed by computing the q-generalized Kolmogorov-Sinai entropy rate in the empirical data of chats as well as in the simulations of interacting emotional agents and Bots.
Ab initio molecular dynamics calculations of ion hydration free energies
We apply ab initio molecular dynamics (AIMD) methods in conjunction with the thermodynamic integration or 'λ-path' technique to compute the intrinsic hydration free energies of Li+, Cl-, and Ag+ ions. Using the Perdew-Burke-Ernzerhof functional, adapting methods developed for classical force field applications, and with consistent assumptions about surface potential (φ) contributions, we obtain absolute AIMD hydration free energies (ΔGhyd) within a few kcal/mol, or better than 4%, of Tissandier et al.'s [J. Phys. Chem. A 102, 7787 (1998)] experimental values augmented with the SPC/E water model φ predictions. The sums of Li+/Cl- and Ag+/Cl- AIMD ΔGhyd, which are not affected by surface potentials, are within 2.6% and 1.2 % of experimental values, respectively. We also report the free energy changes associated with the transition metal ion redox reaction Ag++Ni+→Ag+Ni2+ in water. The predictions for this reaction suggest that existing estimates of ΔGhyd for unstable radiolysis intermediates such as Ni+ may need to be extensively revised.
Pion dose distribution calculations and measurements for dynamic radiotherapy
Routine three dimensional conformation therapy with negative pions is done with the PIOTRON at SIN since two years. More than 60 patients have been treated by spot scan with the 60 converging beams for deep seated tumors in the pelvic region. Extensive measurements have been performed on various phantoms, homogeneous and anthropomorphic, to investigate the influence of tissue inhomogeneities and verify treatment planning calculations. Total dose has been measured by T.E. ionization chambers and TLD, two dimensional stop distributions exposed in planes between phantom slices. In vivo measurements with ionization chambers, as well as catheters filled with /sup 7/LiF TLD's and rolled Al foils, introduced in bladder or rectum, have been used to confirm dose distributions in patients. To check predictions of differences of RBE due to variations in treatment volumes or beam configuration, treatment plans, reflecting typical situation in therapy, have been created for radiobiological investigations. Various user groups have measured biological effects by cell survival experiments with mammalian cells or with mouse intestinal crypt cell assay
Evolutionary dynamics in the Bak-Sneppen model on small-world networks
Kulkarni, R. V.; Almaas, E.; Stroud, D.
1999-01-01
We study the dynamics of the Bak-Sneppen model on small-world networks. For each site in the network, we define a ``connectance,'' which measures the distance to all other sites. We find radically different patterns of activity for different sites, depending on their connectance and also on the topology of the network. For a given network, the site with the minimal connectance shows long periods of stasis interrupted by much smaller periods of activity. In contrast, the activity pattern for t...