Time evolution of the Partridge-Barton Model

The time evolution of the Partridge-Barton model in the presence of the pleiotropic constraint and deleterious somatic mutations is exactly solved for arbitrary fecundity in the context of a matricial formalism. Analytical expressions for the time dependence of the mean survival probabilities are derived. Using the fact that the asymptotic behavior for large time $t$ is controlled by the largest matrix eigenvalue, we obtain the steady state values for the mean survival probabilities and the Malthusian growth exponent. The mean age of the population exhibits a $t^{-1}$ power law decayment. Some Monte Carlo simulations were also performed and they corroborated our theoretical results.Comment: 10 pages, Latex, 1 postscript figure, published in Phys. Rev. E 61, 5664 (2000

Critical behavior for mixed site-bond directed percolation

We study mixed site-bond directed percolation on 2D and 3D lattices by using time-dependent simulations. Our results are compared with rigorous bounds recently obtained by Liggett and by Katori and Tsukahara. The critical fractions $p_{site}^c$ and $p_{bond}^c$ of sites and bonds are extremely well approximated by a relationship reported earlier for isotropic percolation, $(\log p_{site}^c/\log p_{site}^{c^*}+\log p_{bond}^c/\log p_{bond}^{c^*} = 1)$, where $p_{site}^{c^*}$ and $p_{bond}^{c^*}$ are the critical fractions in pure site and bond directed percolation.Comment: 10 pages, figures available on request from [email protected]

Series expansions of the percolation probability for directed square and honeycomb lattices

We have derived long series expansions of the percolation probability for site and bond percolation on directed square and honeycomb lattices. For the square bond problem we have extended the series from 41 terms to 54, for the square site problem from 16 terms to 37, and for the honeycomb bond problem from 13 terms to 36. Analysis of the series clearly shows that the critical exponent $\beta$ is the same for all the problems confirming expectations of universality. For the critical probability and exponent we find in the square bond case, $q_c = 0.3552994\pm 0.0000010$, $\beta = 0.27643\pm 0.00010$, in the square site case $q_c = 0.294515 \pm 0.000005$, $\beta = 0.2763 \pm 0.0003$, and in the honeycomb bond case $q_c = 0.177143 \pm 0.000002$, $\beta = 0.2763 \pm 0.0002$. In addition we have obtained accurate estimates for the critical amplitudes. In all cases we find that the leading correction to scaling term is analytic, i.e., the confluent exponent $\Delta = 1$.Comment: LaTex with epsf, 26 pages, 2 figures and 2 tables in Postscript format included (uufiled). LaTeX version of tables also included for the benefit of those without access to PS printers (note that the tables should be printed in landscape mode). Accepted by J. Phys.

Ac transport studies in polymers by a resistor network and transfer matrix approaches: application to polyaniline

A statistical model of resistor network is proposed to describe a polymer structure and to simulate the real and imaginary components of its ac resistivity. It takes into account the polydispersiveness of the material as well as intrachain and interchain charge transport processes. By the application of a transfer matrix technique, it reproduces ac resistivity measurements carried out with polyaniline films in different doping degrees and at different temperatures. Our results indicate that interchain processes govern the resistivity behavior in the low frequency region while, for higher frequencies, intrachain mechanisms are dominant.Comment: LaTeX file, 15 pages, 5 ps figures, to appear in Phys. Rev.

Coevolution of Glauber-like Ising dynamics on typical networks

We consider coevolution of site status and link structures from two different initial networks: a one dimensional Ising chain and a scale free network. The dynamics is governed by a preassigned stability parameter $S$, and a rewiring factor $\phi$, that determines whether the Ising spin at the chosen site flips or whether the node gets rewired to another node in the system. This dynamics has also been studied with Ising spins distributed randomly among nodes which lie on a network with preferential attachment. We have observed the steady state average stability and magnetisation for both kinds of systems to have an idea about the effect of initial network topology. Although the average stability shows almost similar behaviour, the magnetisation depends on the initial condition we start from. Apart from the local dynamics, the global effect on the dynamics has also been studied. These parameters show interesting variations for different values of $S$ and $\phi$, which helps in determining the steady-state condition for a given substrate.Comment: 8 pages, 10 figure

Exact Solution of an Evolutionary Model without Ageing

We introduce an age-structured asexual population model containing all the relevant features of evolutionary ageing theories. Beneficial as well as deleterious mutations, heredity and arbitrary fecundity are present and managed by natural selection. An exact solution without ageing is found. We show that fertility is associated with generalized forms of the Fibonacci sequence, while mutations and natural selection are merged into an integral equation which is solved by Fourier series. Average survival probabilities and Malthusian growth exponents are calculated indicating that the system may exhibit mutational meltdown. The relevance of the model in the context of fissile reproduction groups as many protozoa and coelenterates is discussed.Comment: LaTeX file, 15 pages, 2 ps figures, to appear in Phys. Rev.

Dynamics of tournaments: the soccer case

A random walk-like model is considered to discuss statistical aspects of tournaments. The model is applied to soccer leagues with emphasis on the scores. This competitive system was computationally simulated and the results are compared with empirical data from the English, the German and the Spanish leagues and showed a good agreement with them. The present approach enabled us to characterize a diffusion where the scores are not normally distributed, having a short and asymmetric tail extending towards more positive values. We argue that this non-Gaussian behavior is related with the difference between the teams and with the asymmetry of the scores system. In addition, we compared two tournament systems: the all-play-all and the elimination tournaments.Comment: To appear in EPJ