Applications and Sexual Version of a Simple Model for Biological Ageing

We use a simple model for biological ageing to study the mortality of the population, obtaining a good agreement with the Gompertz law. We also simulate the same model on a square lattice, considering different strategies of parental care. The results are in agreement with those obtained earlier with the more complicated Penna model for biological ageing. Finally, we present the sexual version of this simple model.Comment: For Int.J.Mod.Phys.C Dec. 2001; 11 pages including 6 fig

Number of spanning clusters at the high-dimensional percolation thresholds

A scaling theory is used to derive the dependence of the average number of spanning clusters at threshold on the lattice size L. This number should become independent of L for dimensions d<6, and vary as log L at d=6. The predictions for d>6 depend on the boundary conditions, and the results there may vary between L^{d-6} and L^0. While simulations in six dimensions are consistent with this prediction (after including corrections of order loglog L), in five dimensions the average number of spanning clusters still increases as log L even up to L = 201. However, the histogram P(k) of the spanning cluster multiplicity does scale as a function of kX(L), with X(L)=1+const/L, indicating that for sufficiently large L the average will approach a finite value: a fit of the 5D multiplicity data with a constant plus a simple linear correction to scaling reproduces the data very well. Numerical simulations for d>6 and for d=4 are also presented.Comment: 8 pages, 11 figures. Final version to appear on Physical Review

Non-Local Product Rules for Percolation

Despite original claims of a first-order transition in the product rule model proposed by Achlioptas et al. [Science 323, 1453 (2009)], recent studies indicate that this percolation model, in fact, displays a continuous transition. The distinctive scaling properties of the model at criticality, however, strongly suggest that it should belong to a different universality class than ordinary percolation. Here we introduce a generalization of the product rule that reveals the effect of non-locality on the critical behavior of the percolation process. Precisely, pairs of unoccupied bonds are chosen according to a probability that decays as a power-law of their Manhattan distance, and only that bond connecting clusters whose product of their sizes is the smallest, becomes occupied. Interestingly, our results for two-dimensional lattices at criticality shows that the power-law exponent of the product rule has a significant influence on the finite-size scaling exponents for the spanning cluster, the conducting backbone, and the cutting bonds of the system. In all three cases, we observe a continuous variation from ordinary to (non-local) explosive percolation exponents.Comment: 5 pages, 4 figure

Resonances in low-energy positron-alkali scattering

Close-coupling calculations were performed with up to five target states at energies in the excitation threshold region for positron scattering from Li, Na and K. Resonances were discovered in the L = 0, 1 and 2 channels in the vicinity of the atomic excitation thresholds. The widths of these resonances vary between 0.2 and 130 MeV. Evidence was found for the existence of positron-alkali bound states in all cases

Positron-alkali atom scattering

Positron-alkali atom scattering was recently investigated both theoretically and experimentally in the energy range from a few eV up to 100 eV. On the theoretical side calculations of the integrated elastic and excitation cross sections as well as total cross sections for Li, Na and K were based upon either the close-coupling method or the modified Glauber approximation. These theoretical results are in good agreement with experimental measurements of the total cross section for both Na and K. Resonance structures were also found in the L = 0, 1 and 2 partial waves for positron scattering from the alkalis. The structure of these resonances appears to be quite complex and, as expected, they occur in conjunction with the atomic excitation thresholds. Currently both theoretical and experimental work is in progress on positron-Rb scattering in the same energy range

On Spatial Consensus Formation: Is the Sznajd Model Different from a Voter Model?

In this paper, we investigate the so-called ``Sznajd Model'' (SM) in one dimension, which is a simple cellular automata approach to consensus formation among two opposite opinions (described by spin up or down). To elucidate the SM dynamics, we first provide results of computer simulations for the spatio-temporal evolution of the opinion distribution $L(t)$, the evolution of magnetization $m(t)$, the distribution of decision times $P(\tau)$ and relaxation times $P(\mu)$. In the main part of the paper, it is shown that the SM can be completely reformulated in terms of a linear VM, where the transition rates towards a given opinion are directly proportional to frequency of the respective opinion of the second-nearest neighbors (no matter what the nearest neighbors are). So, the SM dynamics can be reduced to one rule, ``Just follow your second-nearest neighbor''. The equivalence is demonstrated by extensive computer simulations that show the same behavior between SM and VM in terms of $L(t)$, $m(t)$, $P(\tau)$, $P(\mu)$, and the final attractor statistics. The reformulation of the SM in terms of a VM involves a new parameter $\sigma$, to bias between anti- and ferromagnetic decisions in the case of frustration. We show that $\sigma$ plays a crucial role in explaining the phase transition observed in SM. We further explore the role of synchronous versus asynchronous update rules on the intermediate dynamics and the final attractors. Compared to the original SM, we find three additional attractors, two of them related to an asymmetric coexistence between the opposite opinions.Comment: 22 pages, 20 figures. For related publications see http://www.ais.fraunhofer.de/~fran

Residue network in protein native structure belongs to the universality class of three dimensional critical percolation cluster

A single protein molecule is regarded as a contact network of amino-acid residues. Some studies have indicated that this network is a small world network (SWN), while other results have implied that this is a fractal network (FN). However, SWN and FN are essentially different in the dependence of the shortest path length on the number of nodes. In this paper, we investigate this dependence in the residue contact networks of proteins in native structures, and show that the networks are not SWN but FN. FN is generally characterized by several dimensions. Among them, we focus on three dimensions; the network topological dimension $D_c$, the fractal dimension $D_f$, and the spectral dimension $D_s$. We find that proteins universally yield $D_c \approx 1.9$, $D_f \approx 2.5$ and $Ds \approx 1.3$. These values are in surprisingly good coincidence with those in three dimensional critical percolation cluster. Hence the residue contact networks in the protein native structures belong to the universality class of three dimensional percolation cluster. The criticality is relevant to the ambivalent nature of the protein native structures, i.e., the coexistence of stability and instability, both of which are necessary for a protein to function as a molecular machine or an allosteric enzyme.Comment: 4 pages, 3 figure

Square lattice site percolation at increasing ranges of neighbor interactions

We report site percolation thresholds for square lattice with neighbor interactions at various increasing ranges. Using Monte Carlo techniques we found that nearest neighbors (N$^2$), next nearest neighbors (N$^3$), next next nearest neighbors (N$^4$) and fifth nearest neighbors (N$^6$) yield the same $p_c=0.592...$. At odds, fourth nearest neighbors (N$^5$) give $p_c=0.298...$. These results are given an explanation in terms of symmetry arguments. We then consider combinations of various ranges of interactions with (N$^2$+N$^3$), (N$^2$+N$^4$), (N$^2$+N$^3$+N$^4$) and (N$^2$+N$^5$). The calculated associated thresholds are respectively $p_c=0.407..., 0.337..., 0.288..., 0.234...$. The existing Galam--Mauger universal formula for percolation thresholds does not reproduce the data showing dimension and coordination number are not sufficient to build a universal law which extends to complex lattices.Comment: 4 pages, revtex

Percolation in high dimensions is not understood

The number of spanning clusters in four to nine dimensions does not fully follow the expected size dependence for random percolation.Comment: 9-dimensional data and more points for large lattices added; statistics improved, text expanded, table of exponents inserte