Two site self consistent method for front propagation in reaction-diffusion system

We study front propagation in the reaction diffusion process $A\leftrightarrow2A$ on one dimensional lattice with hard core interaction between the particles. We propose a two site self consistent method (TSSCM) to make analytic estimates for the front velocity and are in excellent agreement with the simulation results for all parameter regimes. We expect that the simplicity of the method will allow one to use this technique for estimating the front velocity in other reaction diffusion processes as well.Comment: 6 figure

Exact Results for a Three-Body Reaction-Diffusion System

A system of particles hopping on a line, singly or as merged pairs, and annihilating in groups of three on encounters, is solved exactly for certain symmetrical initial conditions. The functional form of the density is nearly identical to that found in two-body annihilation, and both systems show non-mean-field, ~1/t**(1/2) instead of ~1/t, decrease of particle density for large times.Comment: 10 page

Universality and tree structure of high energy QCD

Using non-trivial mathematical properties of a class of nonlinear evolution equations, we obtain the universal terms in the asymptotic expansion in rapidity of the saturation scale and of the unintegrated gluon density from the Balitsky-Kovchegov equation. These terms are independent of the initial conditions and of the details of the equation. The last subasymptotic terms are new results and complete the list of all possible universal contributions. Universality is interpreted in a general qualitative picture of high energy scattering, in which a scattering process corresponds to a tree structure probed by a given source.Comment: 4 pages, 3 figure

Model of Cluster Growth and Phase Separation: Exact Results in One Dimension

We present exact results for a lattice model of cluster growth in 1D. The growth mechanism involves interface hopping and pairwise annihilation supplemented by spontaneous creation of the stable-phase, +1, regions by overturning the unstable-phase, -1, spins with probability p. For cluster coarsening at phase coexistence, p=0, the conventional structure-factor scaling applies. In this limit our model falls in the class of diffusion-limited reactions A+A->inert. The +1 cluster size grows diffusively, ~t**(1/2), and the two-point correlation function obeys scaling. However, for p>0, i.e., for the dynamics of formation of stable phase from unstable phase, we find that structure-factor scaling breaks down; the length scale associated with the size of the growing +1 clusters reflects only the short-distance properties of the two-point correlations.Comment: 12 page

A phenomenological theory giving the full statistics of the position of fluctuating pulled fronts

We propose a phenomenological description for the effect of a weak noise on the position of a front described by the Fisher-Kolmogorov-Petrovsky-Piscounov equation or any other travelling wave equation in the same class. Our scenario is based on four hypotheses on the relevant mechanism for the diffusion of the front. Our parameter-free analytical predictions for the velocity of the front, its diffusion constant and higher cumulants of its position agree with numerical simulations.Comment: 10 pages, 3 figure

Symmetry and species segregation in diffusion-limited pair annihilation

We consider a system of q diffusing particle species A_1,A_2,...,A_q that are all equivalent under a symmetry operation. Pairs of particles may annihilate according to A_i + A_j -> 0 with reaction rates k_{ij} that respect the symmetry, and without self-annihilation (k_{ii} = 0). In spatial dimensions d > 2 mean-field theory predicts that the total particle density decays as n(t) ~ 1/t, provided the system remains spatially uniform. We determine the conditions on the matrix k under which there exists a critical segregation dimension d_{seg} below which this uniformity condition is violated; the symmetry between the species is then locally broken. We argue that in those cases the density decay slows down to n(t) ~ t^{-d/d_{seg}} for 2 < d < d_{seg}. We show that when d_{seg} exists, its value can be expressed in terms of the ratio of the smallest to the largest eigenvalue of k. The existence of a conservation law (as in the special two-species annihilation A + B -> 0), although sufficient for segregation, is shown not to be a necessary condition for this phenomenon to occur. We work out specific examples and present Monte Carlo simulations compatible with our analytical results.Comment: latex, 19 pages, 3 eps figures include

Random Energy Model with complex replica number, complex temperatures and classification of the string's phases

The results by E. Gardner and B.Derrida have been enlarged for the complex temperatures and complex numbers of replicas. The phase structure is found. There is a connection with string models and their phase structure is analyzed from the REM's point of view.Comment: 11 pages,revte

Applying bioinformatics for antibody epitope prediction using affinity-selected mimotopes – relevance for vaccine design

To properly characterize protective polyclonal antibody responses, it is necessary to examine epitope specificity. Most antibody epitopes are conformational in nature and, thus, cannot be identified using synthetic linear peptides. Cyclic peptides can function as mimetics of conformational epitopes (termed mimotopes), thereby providing targets, which can be selected by immunoaffinity purification. However, the management of large collections of random cyclic peptides is cumbersome. Filamentous bacteriophage provides a useful scaffold for the expression of random peptides (termed phage display) facilitating both the production and manipulation of complex peptide libraries. Immunoaffinity selection of phage displaying random cyclic peptides is an effective strategy for isolating mimotopes with specificity for a given antiserum. Further epitope prediction based on mimotope sequence is not trivial since mimotopes generally display only small homologies with the target protein. Large numbers of unique mimotopes are required to provide sufficient sequence coverage to elucidate the target epitope. We have developed a method based on pattern recognition theory to deal with the complexity of large collections of conformational mimotopes. The analysis consists of two phases: 1) The learning phase where a large collection of epitope-specific mimotopes is analyzed to identify epitope specific “signs” and 2) The identification phase where immunoaffinity-selected mimotopes are interrogated for the presence of the epitope specific “signs” and assigned to specific epitopes. We are currently using computational methods to define epitope “signs” without the need for prior knowledge of specific mimotopes. This technology provides an important tool for characterizing the breadth of antibody specificities within polyclonal antisera

Exact asymptotics of the freezing transition of a logarithmically correlated random energy model

We consider a logarithmically correlated random energy model, namely a model for directed polymers on a Cayley tree, which was introduced by Derrida and Spohn. We prove asymptotic properties of a generating function of the partition function of the model by studying a discrete time analogy of the KPP-equation - thus translating Bramson's work on the KPP-equation into a discrete time case. We also discuss connections to extreme value statistics of a branching random walk and a rescaled multiplicative cascade measure beyond the critical point

Approach to Asymptotic Behaviour in the Dynamics of the Trapping Reaction

We consider the trapping reaction A + B -> B in space dimension d=1, where the A and B particles have diffusion constants D_A, D_B respectively. We calculate the probability, Q(t), that a given A particle has not yet reacted at time t. Exploiting a recent formulation in which the B particles are eliminated from the problem we find, for t -> \infty, $Q(t) \sim \exp[-(4/\sqrt{\pi})(\rho^2 D_Bt)^{1/2} - (C \rho^2 D_A t)^{1/3} + ...]$, where $\rho$ is the density of B particles and $C \propto D_A/D_B$ for $D_A/D_B << 1$.Comment: 8 pages, 2 figures; minor change