A lattice refinement scheme for finding periodic orbits

A lattice refinement scheme based on the principle of linearized stability is introduced to locate periodic orbits in a two-dimensional map. The method locates all periodic orbits of a specified order within a given starting window and it can be equally well applied when the map is only known implicitly, e.g., as a two-dimensional surface of section arising from a three-dimensional flow. Periodic orbits in the Henon Map, the Predator-Prey Map, the Rossler Flow, and the Lorenz Flow are constructed as illustrations of the method

A certain class of Laplace transforms with applications to reaction and reaction-diffusion equations

A class of Laplace transforms is examined to show that particular cases of this class are associated with production-destruction and reaction-diffusion problems in physics, study of differences of independently distributed random variables and the concept of Laplacianness in statistics, alpha-Laplace and Mittag-Leffler stochastic processes, the concepts of infinite divisibility and geometric infinite divisibility problems in probability theory and certain fractional integrals and fractional derivatives. A number of applications are pointed out with special reference to solutions of fractional reaction and reaction-diffusion equations and their generalizations.Comment: LaTeX, 12 pages, corrected typo

Mean field analysis of Williams-Bjerknes type growth

We investigate a class of stochastic growth models involving competition between two phases in which one of the phases has a competitive advantage. The equilibrium populations of the competing phases are calculated using a mean field analysis. Regression probabilities for the extinction of the advantaged phase are calculated in a leading order approximation. The results of the calculations are in good agreement with simulations carried out on a square lattice with periodic boundaries. The class of models are variants of the Williams- Bjerknes model for the growth of tumours in the basal layer of an epithelium. In the limit in which only one of the phases is unstable the class of models reduces to the well known variants of the Eden model.Comment: 21 pages, Latex2e, Elsevier style, 5 figure

Reaction-diffusion systems and nonlinear waves

The authors investigate the solution of a nonlinear reaction-diffusion equation connected with nonlinear waves. The equation discussed is more general than the one discussed recently by Manne, Hurd, and Kenkre (2000). The results are presented in a compact and elegant form in terms of Mittag-Leffler functions and generalized Mittag-Leffler functions, which are suitable for numerical computation. The importance of the derived results lies in the fact that numerous results on fractional reaction, fractional diffusion, anomalous diffusion problems, and fractional telegraph equations scattered in the literature can be derived, as special cases, of the results investigated in this article.Comment: LaTeX, 16 pages, corrected typo

Solution of generalized fractional reaction-diffusion equations

This paper deals with the investigation of a closed form solution of a generalized fractional reaction-diffusion equation. The solution of the proposed problem is developed in a compact form in terms of the H-function by the application of direct and inverse Laplace and Fourier transforms. Fractional order moments and the asymptotic expansion of the solution are also obtained.Comment: LaTeX, 18 pages, corrected typo

Statistical physics and stromatolite growth: new perspectives on an ancient dilemma

This paper outlines our recent attempts to model the growth and form of microbialites from the perspective of the statistical physics of evolving surfaces. Microbialites arise from the environmental interactions of microbial communities (microbial mats). The mats evolve over time to form internally laminated organosedimentary structures (stromatolites). Modern day stromatolites exist in only a few locations, whereas ancient stromatolitic microbialites were the only form of life for much of the Earth's history. They existed in a wide variety of growth forms, ranging from almost perfect cones to branched columnar structures. The coniform structures are central to the heated debate on the oldest evidence of life. We proposed a biotic model which considers the relationship between upward growth of a phototropic or phototactic biofilm and mineral accretion normal to the surface. These processes are sufficient to account for the growth and form of many ancient stromatolities. These include domical stromatolites and coniform structures with thickened apical zones typical of Conophyton. More angular coniform structures, similar to the stromatolites claimed as the oldest macroscopic evidence of life, form when the photic effects dominate over mineral accretion.Comment: 8 pages, 3 figures. To be published in Proceedings of StatPhys-Taiwan 2004: Biologically Motivated Statistical Physics and Related Problems, 22-26 June 200

Fractional reaction-diffusion equations

In a series of papers, Saxena, Mathai, and Haubold (2002, 2004a, 2004b) derived solutions of a number of fractional kinetic equations in terms of generalized Mittag-Leffler functions which provide the extension of the work of Haubold and Mathai (1995, 2000). The subject of the present paper is to investigate the solution of a fractional reaction-diffusion equation. The results derived are of general nature and include the results reported earlier by many authors, notably by Jespersen, Metzler, and Fogedby (1999) for anomalous diffusion and del-Castillo-Negrete, Carreras, and Lynch (2003) for reaction-diffusion systems with L\'evy flights. The solution has been developed in terms of the H-function in a compact form with the help of Laplace and Fourier transforms. Most of the results obtained are in a form suitable for numerical computation.Comment: LaTeX, 17 pages, corrected typo

A case for biotic morphogenesis of coniform stromatolites

Mathematical models have recently been used to cast doubt on the biotic origin of stromatolites. Here by contrast we propose a biotic model for stromatolite morphogenesis which considers the relationship between upward growth of a phototropic or phototactic biofilm ($v$) and mineral accretion normal to the surface ($\lambda$). These processes are sufficient to account for the growth and form of many ancient stromatolities. Domical stromatolites form when $v$ is less than or comparable to $\lambda$. Coniform structures with thickened apical zones, typical of Conophyton, form when $v >> \lambda$. More angular coniform structures, similar to the stromatolites claimed as the oldest macroscopic evidence of life, form when $v >>> \lambda$.Comment: 10 pages, 3 figures, to appear in Physica

Reaction Front in an A+B -> C Reaction-Subdiffusion Process

We study the reaction front for the process A+B -> C in which the reagents move subdiffusively. Our theoretical description is based on a fractional reaction-subdiffusion equation in which both the motion and the reaction terms are affected by the subdiffusive character of the process. We design numerical simulations to check our theoretical results, describing the simulations in some detail because the rules necessarily differ in important respects from those used in diffusive processes. Comparisons between theory and simulations are on the whole favorable, with the most difficult quantities to capture being those that involve very small numbers of particles. In particular, we analyze the total number of product particles, the width of the depletion zone, the production profile of product and its width, as well as the reactant concentrations at the center of the reaction zone, all as a function of time. We also analyze the shape of the product profile as a function of time, in particular its unusual behavior at the center of the reaction zone

The resistive state in a superconducting wire: Bifurcation from the normal state

We study formally and rigorously the bifurcation to steady and time-periodic states in a model for a thin superconducting wire in the presence of an imposed current. Exploiting the PT-symmetry of the equations at both the linearized and nonlinear levels, and taking advantage of the collision of real eigenvalues leading to complex spectrum, we obtain explicit asymptotic formulas for the stationary solutions, for the amplitude and period of the bifurcating periodic solutions and for the location of their zeros or "phase slip centers" as they are known in the physics literature. In so doing, we construct a center manifold for the flow and give a complete description of the associated finite-dimensional dynamics