J. S. Sá Martins

WHY SEX? MONTE CARLO SIMULATIONS OF SURVIVAL AFTER CATASTROPHES

Using the Penna bit-string model for biological ageing we compare two kinds of reproductive regimes: Sexual reproduction (SR) and meiotic parthenogenesis (MP). The last one is a common type of asexual reproduction with recombination, found in diploid organisms. We show that although both regimes present roughly the same survival rates, the diversity generated by SR is much larger, and can prevent the extinction of a population submitted to a natural disaster. The fixation of bad genes inside an MP population, after many generations, explains our results. We also study the consequences of cloning (simple copy) on population diversity.

Ergodicity in natural earthquake fault networks

Numerical simulations have shown that certain driven nonlinear systems can be characterized by mean-field statistical properties often associated with ergodic dynamics [C. D. Ferguson, W. Klein, and J. B. Rundle, Phys. Rev. E 60, 1359 (1999); D. Egolf, Science 287, 101 (2000)]. These driven mean-field threshold systems feature long-range interactions and can be treated as equilibriumlike systems with statistically stationary dynamics over long time intervals. Recently the equilibrium property of ergodicity was identified in an earthquake fault system, a natural driven threshold system, by means of the Thirumalai-Mountain (TM) fluctuation metric developed in the study of diffusive systems [K. F. Tiampo, J. B. Rundle, W. Klein, J. S. Sá Martins, and C. D. Ferguson, Phys. Rev. Lett. 91, 238501 (2003)]. We analyze the seismicity of three naturally occurring earthquake fault networks from a variety of tectonic settings in an attempt to investigate the range of applicability of effective ergodicity, using the TM metric and other related statistics. Results suggest that, once variations in the catalog data resulting from technical and network issues are accounted for, all of these natural earthquake systems display stationary periods of metastable equilibrium and effective ergodicity that are disrupted by large events. We conclude that a constant rate of events is an important prerequisite for these periods of punctuated ergodicity and that, while the level of temporal variability in the spatial statistics is the controlling factor in the ergodic behavior of seismic networks, no single statistic is sufficient to ensure quantification of ergodicity. Ergodicity in this application not only requires that the system be stationary for these networks at the applicable spatial and temporal scales, but also implies that they are in a state of metastable equilibrium, one in which the ensemble averages can be substituted for temporal averages in studying their spatiotemporal evolution.

Simulated coevolution in a mutating ecology.

The bit-string Penna model is used to simulate the competition between an asexual parthenogenetic and a sexual population sharing the same environment. A newborn of either population can mutate and become a part of the other with some probability. In a stable environment the sexual population soon dies out. When an infestation by rapidly mutation genetically coupled parasites is introduced, however, sexual reproduction prevails, as predicted by the so-called Red Queen hypothesis for the evolution of sex.

Monte Carlo simulation of magnetic systems in the Tsallis statistics

We apply the broad histogram method to an Ising system in the context of the recently reformulated Generalized Thermostatistics, and we claim it to be a very efficient simulation tool for this non-extensive statistics. Results are obtained for the nearest-neighbour version of the Ising model for a range of values of the q parameter of Generalized Thermostatistics. We found an evidence that the 2D-Ising model does not undergo, in the thermodynamical limit, phase transitions at finite temperatures except for the extensive case q=1.

Models of earthquake faults with long-range stress transfer

Seismologists, and more recently computational and condensed-matter physicists, have made extensive use of computer modeling to investigate the physics of earthquakes. The authors describe a simple cellular automaton model that explores the possible relationship between Gutenberg-Richter scaling law and critical phenomena. The presence of scale invariance coupled with the ubiquity of this law in diverse geological conditions indicates that the large-scale statistical properties of earthquakes may be insensitive to environmental details. In this case, relatively simple models will properly capture the essential physics responsible for these behaviors.

Computer simulation of sympatric speciation with Penna ageing model

One species is simulated to split into two separate species via random mutations, even if both populations live together in the same environment. This speciation is achieved in the Penna bitstring model of biological ageing, with modified Verhulst factors, and in part by additional bitstrings regulating phenotype and mate selection.

Justification of sexual reproduction by modified Penna model of ageing

We generalize the standard Penna bit-string model of biological ageing by assuming that each deleterious mutation diminishes the survival probability in every time interval by a small percentage. This effect is added to the usual lethal but age-dependent effect of the same mutation. We then find strong advantages or disadvantages of sexual reproduction (with males and females) compared to asexual cloning, depending on parameters.

Evolution of biodiversity and sympatric speciation through competition in a unimodal distribution of resources

A microscopic agent-based dynamical model for diploid age-structured populations is used to study the evolution of biodiversity and sympatric speciation. The underlying ecology is represented by a unimodal distribution of resources of some width. Competition among individuals is also described by a similar distribution, and its strength is maximum for individuals with the same phenotype and decreases with distance in phenotype space as a Gaussian, with some width. These two widths define the models phase space, in which we identify the regions where an autonomous emergence of stable biodiversity or speciation is more likely.

Simulated ecology-driven sympatric speciation

We introduce a multilocus genetically acquired phenotype, submitted to mutations and with selective value, in an age-structured model for biological aging. This phenotype describes a single-trait effect of the environment on an individual, and we study the resulting distribution of this trait among the population. In particular, our simulations show that the appearance of a double phenotypic attractor in the ecology induces the emergence of a stable polymorphism, as observed in the Galapagos finches. In the presence of this polymorphism, the simulations generate short-term speciation, when mating preferences are also allowed to suffer mutations and acquire selective value.

Thermodynamic behavior of a phase transition in a model for sympatric speciation

We investigate the macroscopic effects of the ingredients that drive the origin of species through sympatric speciation. In our model, sympatric speciation is obtained as we tune up the strength of competition between individuals with different phenotypes. As a function of this control parameter, we can characterize, through the behavior of a macroscopic order parameter, a phase transition from a nonspeciation to a speciation state of the system. The behavior of the first derivative of the order parameter with respect to the control parameter is consistent with a phase transition and exhibits a sharp peak at the transition point. For different resources distribution, the transition point is shifted, an effect similar to pressure in a PVT system. The inverse of the parameter related to a sexual selection strength behaves like an external field in the system and, as thus, is also a control parameter. The macroscopic effects of the biological parameters used in our model are a reminiscent of the behavior of thermodynamic quantities in a phase transition of an equilibrium physical system.