-
Universal asymptotic clone size distribution for general population growth
Authors:
Michael D. Nicholson,
Tibor Antal
Abstract:
Deterministically growing (wild-type) populations which seed stochastically developing mutant clones have found an expanding number of applications from microbial populations to cancer. The special case of exponential wild-type population growth, usually termed the Luria-Delbrück or Lea-Coulson model, is often assumed but seldom realistic. In this article we generalise this model to different type…
▽ More
Deterministically growing (wild-type) populations which seed stochastically developing mutant clones have found an expanding number of applications from microbial populations to cancer. The special case of exponential wild-type population growth, usually termed the Luria-Delbrück or Lea-Coulson model, is often assumed but seldom realistic. In this article we generalise this model to different types of wild-type population growth, with mutants evolving as a birth-death branching process. Our focus is on the size distribution of clones - that is the number of progeny of a founder mutant - which can be mapped to the total number of mutants. Exact expressions are derived for exponential, power-law and logistic population growth. Additionally for a large class of population growth we prove that the long time limit of the clone size distribution has a general two-parameter form, whose tail decays as a power-law. Considering metastases in cancer as the mutant clones, upon analysing a data-set of their size distribution, we indeed find that a power-law tail is more likely than an exponential one.
△ Less
Submitted 26 October, 2016; v1 submitted 17 April, 2016;
originally announced April 2016.
-
Mutant number distribution in an exponentially growing population
Authors:
Peter Keller,
Tibor Antal
Abstract:
We present an explicit solution to a classic model of cell-population growth introduced by Luria and Delbrueck 70 years ago to study the emergence of mutations in bacterial populations. In this model a wild-type population is assumed to grow exponentially in a deterministic fashion. Proportional to the wild-type population size, mutants arrive randomly and initiate new sub-populations of mutants t…
▽ More
We present an explicit solution to a classic model of cell-population growth introduced by Luria and Delbrueck 70 years ago to study the emergence of mutations in bacterial populations. In this model a wild-type population is assumed to grow exponentially in a deterministic fashion. Proportional to the wild-type population size, mutants arrive randomly and initiate new sub-populations of mutants that grows stochastically according to a supercritical birth and death process. We give an exact expression for the generating function of the total number of mutants at a given wild type population size. We present a simple expression for the probability of finding no mutants, and a recursion formula for the probability of finding a given number of mutants. In the "large population-small mutation"-limit we recover recent results of Kessler and Levin for a fully stochastic version of the process.
△ Less
Submitted 13 October, 2014;
originally announced October 2014.
-
Spatial evolution of tumors with successive driver mutations
Authors:
Tibor Antal,
P. L. Krapivsky,
M. A. Nowak
Abstract:
We study the spatial evolutionary dynamics of solid tumors as they obtain additional driver mutations. We start with a cancer clone that expands uniformly in three dimensions giving rise to a spherical shape. We assume that cell division occurs on the surface of the growing tumor. Each cell division has a chance to give rise to a mutation that activates an additional driver gene. The resulting clo…
▽ More
We study the spatial evolutionary dynamics of solid tumors as they obtain additional driver mutations. We start with a cancer clone that expands uniformly in three dimensions giving rise to a spherical shape. We assume that cell division occurs on the surface of the growing tumor. Each cell division has a chance to give rise to a mutation that activates an additional driver gene. The resulting clone has an enhanced growth rate, which generates a local ensemble of faster growing cells, thereby distorting the spherical shape of the tumor. We derive analytic formulas for the geometric boundary that separates the original cancer clone from the new mutant as well as the expanding frontier of the new mutant. The total number of original cancer cells converges to a constant as time goes to infinity, because this clone becomes enveloped by mutants. We derive formulas for the abundance and diversity of additional driver mutations as function of time. Our model is semi-deterministic: the spatial growth of the various cancer clones follows deterministic equations, but the arrival of a new mutant is a stochastic event.
△ Less
Submitted 7 August, 2013;
originally announced August 2013.
-
Molecular Spiders on the Plane
Authors:
Tibor Antal,
P. L. Krapivsky
Abstract:
Synthetic bio-molecular spiders with "legs" made of single-stranded segments of DNA can move on a surface covered by single-stranded segments of DNA called substrates when the substrate DNA is complementary to the leg DNA. If the motion of a spider does not affect the substrates, the spider behaves asymptotically as a random walk. We study the diffusion coefficient and the number of visited sites…
▽ More
Synthetic bio-molecular spiders with "legs" made of single-stranded segments of DNA can move on a surface covered by single-stranded segments of DNA called substrates when the substrate DNA is complementary to the leg DNA. If the motion of a spider does not affect the substrates, the spider behaves asymptotically as a random walk. We study the diffusion coefficient and the number of visited sites for spiders moving on the square lattice with a substrate in each lattice site. The spider's legs hop to nearest-neighbor sites with the constraint that the distance between any two legs cannot exceed a maximal span. We establish analytic results for bipedal spiders, and investigate multileg spiders numerically. In experimental realizations legs usually convert substrates into products (visited sites). The binding of legs to products is weaker, so the hopping rate from the substrates is smaller. This makes the problem non-Markovian and we investigate it numerically. We demonstrate the emergence of a counter-intuitive behavior - the more spiders are slowed down on unvisited sites, the more motile they become.
△ Less
Submitted 20 July, 2012; v1 submitted 28 April, 2012;
originally announced April 2012.
-
Outbreak size distributions in epidemics with multiple stages
Authors:
Tibor Antal,
P. L. Krapivsky
Abstract:
Multiple-type branching processes that model the spread of infectious diseases are investigated. In these stochastic processes, the disease goes through multiple stages before it eventually disappears. We mostly focus on the critical multistage Susceptible-Infected-Recovered (SIR) infection process. In the infinite population limit, we compute the outbreak size distributions and show that asymptot…
▽ More
Multiple-type branching processes that model the spread of infectious diseases are investigated. In these stochastic processes, the disease goes through multiple stages before it eventually disappears. We mostly focus on the critical multistage Susceptible-Infected-Recovered (SIR) infection process. In the infinite population limit, we compute the outbreak size distributions and show that asymptotic results apply to more general multiple-type critical branching processes. Finally using heuristic arguments and simulations we establish scaling laws for a multistage SIR model in a finite population.
△ Less
Submitted 20 July, 2012; v1 submitted 18 April, 2012;
originally announced April 2012.
-
Exact solution of a two-type branching process: Models of tumor progression
Authors:
Tibor Antal,
P. L. Krapivsky
Abstract:
An explicit solution for a general two-type birth-death branching process with one way mutation is presented. This continuous time process mimics the evolution of resistance to treatment, or the onset of an extra driver mutation during tumor progression. We obtain the exact generating function of the process at arbitrary times, and derive various large time scaling limits. In the simultaneous smal…
▽ More
An explicit solution for a general two-type birth-death branching process with one way mutation is presented. This continuous time process mimics the evolution of resistance to treatment, or the onset of an extra driver mutation during tumor progression. We obtain the exact generating function of the process at arbitrary times, and derive various large time scaling limits. In the simultaneous small mutation rate and large time scaling limit, the distribution of the mutant cells develops some atypical properties, including a power law tail and diverging average.
△ Less
Submitted 5 May, 2011;
originally announced May 2011.
-
First Passage Properties of the Polya Urn Process
Authors:
Tibor Antal,
E. Ben-Naim,
P. L. Krapivsky
Abstract:
We study first passage statistics of the Polya urn model. In this random process, the urn contains two types of balls. In each step, one ball is drawn randomly from the urn, and subsequently placed back into the urn together with an additional ball of the same type. We derive the probability G_n that the two types of balls are equal in number, for the first time, when there is a total of 2n balls.…
▽ More
We study first passage statistics of the Polya urn model. In this random process, the urn contains two types of balls. In each step, one ball is drawn randomly from the urn, and subsequently placed back into the urn together with an additional ball of the same type. We derive the probability G_n that the two types of balls are equal in number, for the first time, when there is a total of 2n balls. This first passage probability decays algebraically, G_n ~ n^{-2}, when n is large. We also derive the probability that a tie ever happens. This probability is between zero and one, so that a tie may occur in some realizations but not in others. The likelihood of a tie is appreciable only if the initial difference in the number balls is of the order of the square-root of the total number of balls.
△ Less
Submitted 5 May, 2010;
originally announced May 2010.
-
Galaxy distribution and extreme value statistics
Authors:
Tibor Antal,
Francesco Sylos Labini,
Nikolay L. Vasilyev,
Yurij V. Baryshev
Abstract:
We consider the conditional galaxy density around each galaxy, and study its fluctuations in the newest samples of the Sloan Digital Sky Survey Data Release 7. Over a large range of scales, both the average conditional density and its variance show a nontrivial scaling behavior, which resembles to criticality. The density depends, for 10 < r < 80 Mpc/h, only weakly (logarithmically) on the syste…
▽ More
We consider the conditional galaxy density around each galaxy, and study its fluctuations in the newest samples of the Sloan Digital Sky Survey Data Release 7. Over a large range of scales, both the average conditional density and its variance show a nontrivial scaling behavior, which resembles to criticality. The density depends, for 10 < r < 80 Mpc/h, only weakly (logarithmically) on the system size. Correspondingly, we find that the density fluctuations follow the Gumbel distribution of extreme value statistics. This distribution is clearly distinguishable from a Gaussian distribution, which would arise for a homogeneous spatial galaxy configuration. We also point out similarities between the galaxy distribution and critical systems of statistical physics.
△ Less
Submitted 8 September, 2009;
originally announced September 2009.
-
Exact solution of a two-type branching process: Clone size distribution in cell division kinetics
Authors:
Tibor Antal,
P. L. Krapivsky
Abstract:
We study a two-type branching process which provides excellent description of experimental data on cell dynamics in skin tissue (Clayton et al., 2007). The model involves only a single type of progenitor cell, and does not require support from a self-renewed population of stem cells. The progenitor cells divide and may differentiate into post-mitotic cells. We derive an exact solution of this mo…
▽ More
We study a two-type branching process which provides excellent description of experimental data on cell dynamics in skin tissue (Clayton et al., 2007). The model involves only a single type of progenitor cell, and does not require support from a self-renewed population of stem cells. The progenitor cells divide and may differentiate into post-mitotic cells. We derive an exact solution of this model in terms of generating functions for the total number of cells, and for the number of cells of different types. We also deduce large time asymptotic behaviors drawing on our exact results, and on an independent diffusion approximation.
△ Less
Submitted 4 August, 2009;
originally announced August 2009.
-
Shepherd Model for Knot-Limited Polymer Ejection from a Capsid
Authors:
T. Antal,
P. L. Krapivsky,
S. Redner
Abstract:
We construct a tractable model to describe the rate at which a knotted polymer is ejected from a spherical capsid via a small pore. Knots are too large to fit through the pore and must reptate to the end of the polymer for ejection to occur. The reptation of knots is described by symmetric exclusion on the line, with the internal capsid pressure represented by an additional biased particle that…
▽ More
We construct a tractable model to describe the rate at which a knotted polymer is ejected from a spherical capsid via a small pore. Knots are too large to fit through the pore and must reptate to the end of the polymer for ejection to occur. The reptation of knots is described by symmetric exclusion on the line, with the internal capsid pressure represented by an additional biased particle that drives knots to the end of the chain. We compute the exact ejection speed for a finite number of knots L and find that it scales as 1/L. We also construct a continuum theory for many knots that matches the exact discrete theory for large L.
△ Less
Submitted 11 August, 2009; v1 submitted 6 May, 2009;
originally announced May 2009.
-
Logarithmic current fluctuations in non-equilibrium quantum spin chains
Authors:
Tibor Antal,
P. L. Krapivsky,
A. Rákos
Abstract:
We study zero-temperature quantum spin chains which are characterized by a non-vanishing current. For the XX model starting from the initial state |... + + + - - - ...> we derive an exact expression for the variance of the total spin current. We show that asymptotically the variance exhibits an anomalously slow logarithmic growth; we also extract the sub-leading constant term. We then argue that…
▽ More
We study zero-temperature quantum spin chains which are characterized by a non-vanishing current. For the XX model starting from the initial state |... + + + - - - ...> we derive an exact expression for the variance of the total spin current. We show that asymptotically the variance exhibits an anomalously slow logarithmic growth; we also extract the sub-leading constant term. We then argue that the logarithmic growth remains valid for the XXZ model in the critical region.
△ Less
Submitted 9 January, 2009; v1 submitted 26 August, 2008;
originally announced August 2008.
-
Front Propagation with Rejuvenation in Flipping Processes
Authors:
T. Antal,
D. ben-Avraham,
E. Ben-Naim,
P. L. Krapivsky
Abstract:
We study a directed flipping process that underlies the performance of the random edge simplex algorithm. In this stochastic process, which takes place on a one-dimensional lattice whose sites may be either occupied or vacant, occupied sites become vacant at a constant rate and simultaneously cause all sites to the right to change their state. This random process exhibits rich phenomenology. Fir…
▽ More
We study a directed flipping process that underlies the performance of the random edge simplex algorithm. In this stochastic process, which takes place on a one-dimensional lattice whose sites may be either occupied or vacant, occupied sites become vacant at a constant rate and simultaneously cause all sites to the right to change their state. This random process exhibits rich phenomenology. First, there is a front, defined by the position of the left-most occupied site, that propagates at a nontrivial velocity. Second, the front involves a depletion zone with an excess of vacant sites. The total excess D_k increases logarithmically, D_k ~ ln k, with the distance k from the front. Third, the front exhibits rejuvenation -- young fronts are vigorous but old fronts are sluggish. We investigate these phenomena using a quasi-static approximation, direct solutions of small systems, and numerical simulations.
△ Less
Submitted 1 August, 2008;
originally announced August 2008.
-
Exciting Hard Spheres
Authors:
T. Antal,
P. L. Krapivsky,
S. Redner
Abstract:
We investigate the collision cascade that is generated by a single moving incident particle on a static hard-sphere gas. We argue that the number of moving particles at time t grows as t^{xi} and the number collisions up to time t grows as t^{eta}, with xi=2d/(d+2) and eta=2(d+1)/(d+2) and d the spatial dimension. These growth laws are the same as those from a hydrodynamic theory for the shock w…
▽ More
We investigate the collision cascade that is generated by a single moving incident particle on a static hard-sphere gas. We argue that the number of moving particles at time t grows as t^{xi} and the number collisions up to time t grows as t^{eta}, with xi=2d/(d+2) and eta=2(d+1)/(d+2) and d the spatial dimension. These growth laws are the same as those from a hydrodynamic theory for the shock wave emanating from an explosion. Our predictions are verified by molecular dynamics simulations in d=1 and 2. For a particle incident on a static gas in a half-space, the resulting backsplatter ultimately contains almost all the initial energy.
△ Less
Submitted 25 May, 2008;
originally announced May 2008.
-
Voter Models on Heterogeneous Networks
Authors:
V. Sood,
Tibor Antal,
S. Redner
Abstract:
We study simple interacting particle systems on heterogeneous networks, including the voter model and the invasion process. These are both two-state models in which in an update event an individual changes state to agree with a neighbor. For the voter model, an individual "imports" its state from a randomly-chosen neighbor. Here the average time T_N to reach consensus for a network of N nodes wi…
▽ More
We study simple interacting particle systems on heterogeneous networks, including the voter model and the invasion process. These are both two-state models in which in an update event an individual changes state to agree with a neighbor. For the voter model, an individual "imports" its state from a randomly-chosen neighbor. Here the average time T_N to reach consensus for a network of N nodes with an uncorrelated degree distribution scales as N mu_1^2/mu_2, where mu_k is the kth moment of the degree distribution. Quick consensus thus arises on networks with broad degree distributions. We also identify the conservation law that characterizes the route by which consensus is reached. Parallel results are derived for the invasion process, in which the state of an agent is "exported" to a random neighbor. We further generalize to biased dynamics in which one state is favored. The probability for a single fitter mutant located at a node of degree k to overspread the population--the fixation probability--is proportional to k for the voter model and to 1/k for the invasion process.
△ Less
Submitted 22 April, 2008; v1 submitted 27 December, 2007;
originally announced December 2007.
-
Guiding-fields for phase-separation: Controlling Liesegang patterns
Authors:
Tibor Antal,
Ioana Bena,
Michel Droz,
Kirsten Martens,
Zoltan Racz
Abstract:
Liesegang patterns emerge from precipitation processes and may be used to build bulk structures at submicron lengthscales. Thus they have significant potential for technological applications provided adequate methods of control can be devised. Here we describe a simple, physically realizable pattern-control based on the notion of driven precipitation, meaning that the phase-separation is governe…
▽ More
Liesegang patterns emerge from precipitation processes and may be used to build bulk structures at submicron lengthscales. Thus they have significant potential for technological applications provided adequate methods of control can be devised. Here we describe a simple, physically realizable pattern-control based on the notion of driven precipitation, meaning that the phase-separation is governed by a guiding field such as, for example, a temperature or a pH field. The phase-separation is modeled through a non-autonomous Cahn-Hilliard equation whose spinodal is determined by the evolving guiding field. Control over the dynamics of the spinodal gives control over the velocity of the instability front which separates the stable and unstable regions of the system. Since the wavelength of the pattern is largely determined by this velocity, the distance between successive precipitation bands becomes controllable. We demonstrate the above ideas by numerical studies of a 1D system with diffusive guiding field. We find that the results can be accurately described by employing a linear stability analysis (pulled-front theory) for determining the velocity -- local-wavelength relationship. From the perspective of the Liesegang theory, our results indicate that the so-called revert patterns may be naturally generated by diffusive guiding fields.
△ Less
Submitted 12 August, 2007; v1 submitted 5 June, 2007;
originally announced June 2007.
-
Molecular Spiders with Memory
Authors:
Tibor Antal,
P. L. Krapivsky
Abstract:
Synthetic bio-molecular spiders with "legs" made of single-stranded segments of DNA can move on a surface which is also covered by single-stranded segments of DNA complementary to the leg DNA. In experimental realizations, when a leg detaches from a segment of the surface for the first time it alters that segment, and legs subsequently bound to these altered segments more weakly. Inspired by the…
▽ More
Synthetic bio-molecular spiders with "legs" made of single-stranded segments of DNA can move on a surface which is also covered by single-stranded segments of DNA complementary to the leg DNA. In experimental realizations, when a leg detaches from a segment of the surface for the first time it alters that segment, and legs subsequently bound to these altered segments more weakly. Inspired by these experiments we investigate spiders moving along a one-dimensional substrate, whose legs leave newly visited sites at a slower rate than revisited sites. For a random walk (one-leg spider) the slowdown does not effect the long time behavior. For a bipedal spider, however, the slowdown generates an effective bias towards unvisited sites, and the spider behaves similarly to the excited walk. Surprisingly, the slowing down of the spider at new sites increases the diffusion coefficient and accelerates the growth of the number of visited sites.
△ Less
Submitted 17 May, 2007;
originally announced May 2007.
-
Molecular Spiders in One Dimension
Authors:
Tibor Antal,
P. L. Krapivsky,
Kirone Mallick
Abstract:
Molecular spiders are synthetic bio-molecular systems which have "legs" made of short single-stranded segments of DNA. Spiders move on a surface covered with single-stranded DNA segments complementary to legs. Different mappings are established between various models of spiders and simple exclusion processes. For spiders with simple gait and varying number of legs we compute the diffusion coeffi…
▽ More
Molecular spiders are synthetic bio-molecular systems which have "legs" made of short single-stranded segments of DNA. Spiders move on a surface covered with single-stranded DNA segments complementary to legs. Different mappings are established between various models of spiders and simple exclusion processes. For spiders with simple gait and varying number of legs we compute the diffusion coefficient; when the hopping is biased we also compute their velocity.
△ Less
Submitted 18 May, 2007;
originally announced May 2007.
-
Dynamics of Microtubule Instabilities
Authors:
T. Antal,
P. L. Krapivsky,
S. Redner
Abstract:
We investigate the dynamics of an idealized model of microtubule growth that evolves by: (i) attachment of guanosine triphosphate (GTP) at rate lambda, (ii) conversion of GTP to guanosine diphosphate (GDP) at rate 1, and (iii) detachment of GDP at rate mu. As a function of these rates, a microtubule can grow steadily or its length can fluctuate wildly. For mu=0, we find the exact tubule and GTP ca…
▽ More
We investigate the dynamics of an idealized model of microtubule growth that evolves by: (i) attachment of guanosine triphosphate (GTP) at rate lambda, (ii) conversion of GTP to guanosine diphosphate (GDP) at rate 1, and (iii) detachment of GDP at rate mu. As a function of these rates, a microtubule can grow steadily or its length can fluctuate wildly. For mu=0, we find the exact tubule and GTP cap length distributions, and power-law length distributions of GTP and GDP islands. For mu=infinity, we argue that the time between catastrophes, where the microtubule shrinks to zero length, scales as exp(lambda). We also find the phase boundary between a growing and shrinking microtubule.
△ Less
Submitted 23 November, 2010; v1 submitted 14 March, 2007;
originally announced March 2007.
-
Dynamics of Microtubule Growth and Catastrophe
Authors:
T. Antal,
P. L. Krapivsky,
S. Redner,
M. Mailman,
B. Chakraborty
Abstract:
We investigate a simple model of microtubule dynamics in which a microtubule evolves by: (i) attachment of guanosine triphosphate (GTP) to its end at rate lambda, (ii) GTP converting irreversibly to guanosine diphosphate (GDP) at rate 1, and (iii) detachment of GDP from the end of a microtubule at rate mu. As a function of these elemental rates, the microtubule can grow steadily or its length ca…
▽ More
We investigate a simple model of microtubule dynamics in which a microtubule evolves by: (i) attachment of guanosine triphosphate (GTP) to its end at rate lambda, (ii) GTP converting irreversibly to guanosine diphosphate (GDP) at rate 1, and (iii) detachment of GDP from the end of a microtubule at rate mu. As a function of these elemental rates, the microtubule can grow steadily or its length can fluctuate wildly. A master equation approach is developed to characterize these intriguing features. For mu=0, we find exact expressions for tubule and GTP cap length distributions, as well as a power-law length distributions of GTP and GDP islands. For mu=oo, we find the average time between catastrophes, where the microtubule shrinks to zero length, and extend this approach to also determine the size distribution of avalanches (sequence of consecutive GDP detachment events). We obtain the phase diagram for general rates and verify our predictions by numerical simulations.
△ Less
Submitted 22 April, 2008; v1 submitted 1 March, 2007;
originally announced March 2007.
-
Exact steady-state velocity of ratchets driven by random sequential adsorption
Authors:
Maria R. D'Orsogna,
Tom Chou,
Tibor Antal
Abstract:
We solve the problem of discrete translocation of a polymer through a pore, driven by the irreversible, random sequential adsorption of particles on one side of the pore. Although the kinetics of the wall motion and the deposition are coupled, we find the exact steady-state distribution for the gap between the wall and the nearest deposited particle. This result enables us to construct the mean…
▽ More
We solve the problem of discrete translocation of a polymer through a pore, driven by the irreversible, random sequential adsorption of particles on one side of the pore. Although the kinetics of the wall motion and the deposition are coupled, we find the exact steady-state distribution for the gap between the wall and the nearest deposited particle. This result enables us to construct the mean translocation velocity demonstrating that translocation is faster when the adsorbing particles are smaller. Monte-Carlo simulations also show that smaller particles gives less dispersion in the ratcheted motion. We also define and compare the relative efficiencies of ratcheting by deposition of particles with different sizes and we describe an associated "zone-refinement" process.
△ Less
Submitted 18 December, 2008; v1 submitted 28 March, 2007;
originally announced March 2007.
-
Aging and Immortality in a Cell Proliferation Model
Authors:
T. Antal,
K. B. Blagoev,
S. A. Trugman,
S. Redner
Abstract:
We investigate a model of cell division in which the length of telomeres within the cell regulate their proliferative potential. At each cell division the ends of linear chromosomes change and a cell becomes senescent when one or more of its telomeres become shorter than a critical length. In addition to this systematic shortening, exchange of telomere DNA between the two daughter cells can occu…
▽ More
We investigate a model of cell division in which the length of telomeres within the cell regulate their proliferative potential. At each cell division the ends of linear chromosomes change and a cell becomes senescent when one or more of its telomeres become shorter than a critical length. In addition to this systematic shortening, exchange of telomere DNA between the two daughter cells can occur at each cell division. We map this telomere dynamics onto a biased branching diffusion process with an absorbing boundary condition whenever any telomere reaches the critical length. As the relative effects of telomere shortening and cell division are varied, there is a phase transition between finite lifetime and infinite proliferation of the cell population. Using simple first-passage ideas, we quantify the nature of this transition.
△ Less
Submitted 22 April, 2008; v1 submitted 25 September, 2006;
originally announced September 2006.
-
Flows on Graphs with Random Capacities
Authors:
T. Antal,
P. L. Krapivsky
Abstract:
We investigate flows on graphs whose links have random capacities. For binary trees we derive the probability distribution for the maximal flow from the root to a leaf, and show that for infinite trees it vanishes beyond a certain threshold that depends on the distribution of capacities. We then examine the maximal total flux from the root to the leaves. Our methods generalize to simple graphs w…
▽ More
We investigate flows on graphs whose links have random capacities. For binary trees we derive the probability distribution for the maximal flow from the root to a leaf, and show that for infinite trees it vanishes beyond a certain threshold that depends on the distribution of capacities. We then examine the maximal total flux from the root to the leaves. Our methods generalize to simple graphs with loops, e.g., to hierarchical lattices and to complete graphs.
△ Less
Submitted 6 September, 2006;
originally announced September 2006.
-
Evolutionary dynamics on degree-heterogeneous graphs
Authors:
T. Antal,
S. Redner,
V. Sood
Abstract:
The evolution of two species with different fitness is investigated on degree-heterogeneous graphs. The population evolves either by one individual dying and being replaced by the offspring of a random neighbor (voter model (VM) dynamics) or by an individual giving birth to an offspring that takes over a random neighbor node (invasion process (IP) dynamics). The fixation probability for one spec…
▽ More
The evolution of two species with different fitness is investigated on degree-heterogeneous graphs. The population evolves either by one individual dying and being replaced by the offspring of a random neighbor (voter model (VM) dynamics) or by an individual giving birth to an offspring that takes over a random neighbor node (invasion process (IP) dynamics). The fixation probability for one species to take over a population of N individuals depends crucially on the dynamics and on the local environment. Starting with a single fitter mutant at a node of degree k, the fixation probability is proportional to k for VM dynamics and to 1/k for IP dynamics.
△ Less
Submitted 11 May, 2006; v1 submitted 21 January, 2006;
originally announced January 2006.
-
Escape of a Uniform Random Walk from an Interval
Authors:
T. Antal,
S. Redner
Abstract:
We study the first-passage properties of a random walk in the unit interval in which the length of a single step is uniformly distributed over the finite range [-a,a]. For a of the order of one, the exit probabilities to each edge of the interval and the exit time from the interval exhibit anomalous properties stemming from the change in the minimum number of steps to escape the interval as a fu…
▽ More
We study the first-passage properties of a random walk in the unit interval in which the length of a single step is uniformly distributed over the finite range [-a,a]. For a of the order of one, the exit probabilities to each edge of the interval and the exit time from the interval exhibit anomalous properties stemming from the change in the minimum number of steps to escape the interval as a function of the starting point. As a decreases, first-passage properties approach those of continuum diffusion, but non-diffusive effects remain because of residual discreteness effects
△ Less
Submitted 22 December, 2005;
originally announced December 2005.
-
Fixation of strategies for an evolutionary game in finite populations
Authors:
Tibor Antal,
Istvan Scheuring
Abstract:
A stochastic evolutionary dynamics of two strategies given by 2 x 2 matrix games is studied in finite populations. We focus on stochastic properties of fixation: how a strategy represented by a single individual wins over the entire population. The process is discussed in the framework of a random walk with arbitrary hopping rates. The time of fixation is found to be identical for both strategie…
▽ More
A stochastic evolutionary dynamics of two strategies given by 2 x 2 matrix games is studied in finite populations. We focus on stochastic properties of fixation: how a strategy represented by a single individual wins over the entire population. The process is discussed in the framework of a random walk with arbitrary hopping rates. The time of fixation is found to be identical for both strategies in any particular game. The asymptotic behavior of the fixation time and fixation probabilities in the large population size limit is also discussed. We show that fixation is fast when there is at least one pure evolutionary stable strategy (ESS) in the infinite population size limit, while fixation is slow when the ESS is the coexistence of the two strategies.
△ Less
Submitted 8 September, 2005;
originally announced September 2005.
-
Dynamics of Social Balance on Networks
Authors:
T. Antal,
P. L. Krapivsky,
S. Redner
Abstract:
We study the evolution of social networks that contain both friendly and unfriendly pairwise links between individual nodes. The network is endowed with dynamics in which the sense of a link in an imbalanced triad--a triangular loop with 1 or 3 unfriendly links--is reversed to make the triad balanced. With this dynamics, an infinite network undergoes a dynamic phase transition from a steady stat…
▽ More
We study the evolution of social networks that contain both friendly and unfriendly pairwise links between individual nodes. The network is endowed with dynamics in which the sense of a link in an imbalanced triad--a triangular loop with 1 or 3 unfriendly links--is reversed to make the triad balanced. With this dynamics, an infinite network undergoes a dynamic phase transition from a steady state to "paradise"--all links are friendly--as the propensity p for friendly links in an update event passes through 1/2. A finite network always falls into a socially-balanced absorbing state where no imbalanced triads remain. If the additional constraint that the number of imbalanced triads in the network does not increase in an update is imposed, then the network quickly reaches a balanced final state.
△ Less
Submitted 19 June, 2005;
originally announced June 2005.
-
A "Burnt Bridge'' Brownian Ratchet
Authors:
T. Antal,
P. L. Krapivsky
Abstract:
Motivated by a biased diffusion of molecular motors with the bias dependent on the state of the substrate, we investigate a random walk on a one-dimensional lattice that contains weak links (called "bridges'') which are affected by the walker. Namely, a bridge is destroyed with probability p when the walker crosses it; the walker is not allowed to cross it again and this leads to a directed moti…
▽ More
Motivated by a biased diffusion of molecular motors with the bias dependent on the state of the substrate, we investigate a random walk on a one-dimensional lattice that contains weak links (called "bridges'') which are affected by the walker. Namely, a bridge is destroyed with probability p when the walker crosses it; the walker is not allowed to cross it again and this leads to a directed motion. The velocity of the walker is determined analytically for equidistant bridges. The special case of p=1 is more tractable -- both the velocity and the diffusion constant are calculated for uncorrelated locations of bridges, including periodic and random distributions.
△ Less
Submitted 25 April, 2005;
originally announced April 2005.
-
Excited Random Walk in One Dimension
Authors:
T. Antal,
S. Redner
Abstract:
We study the excited random walk, in which a walk that is at a site that contains cookies eats one cookie and then hops to the right with probability p and to the left with probability q=1-p. If the walk hops onto an empty site, there is no bias. For the 1-excited walk on the half-line (one cookie initially at each site), the probability of first returning to the starting point at time t scales…
▽ More
We study the excited random walk, in which a walk that is at a site that contains cookies eats one cookie and then hops to the right with probability p and to the left with probability q=1-p. If the walk hops onto an empty site, there is no bias. For the 1-excited walk on the half-line (one cookie initially at each site), the probability of first returning to the starting point at time t scales as t^{-(2-p)}. Although the average return time to the origin is infinite for all p, the walk eats, on average, only a finite number of cookies until this first return when p<1/2. For the infinite line, the probability distribution for the 1-excited walk has an unusual anomaly at the origin. The positions of the leftmost and rightmost uneaten cookies can be accurately estimated by probabilistic arguments and their corresponding distributions have power-law singularities near the origin. The 2-excited walk on the infinite line exhibits peculiar features in the regime p>3/4, where the walk is transient, including a mean displacement that grows as t^{nu}, with nu>1/2 dependent on p, and a breakdown of scaling for the probability distribution of the walk.
△ Less
Submitted 20 December, 2004;
originally announced December 2004.
-
Weight-driven growing networks
Authors:
T. Antal,
P. L. Krapivsky
Abstract:
We study growing networks in which each link carries a certain weight (randomly assigned at birth and fixed thereafter). The weight of a node is defined as the sum of the weights of the links attached to the node, and the network grows via the simplest weight-driven rule: A newly-added node is connected to an already existing node with the probability which is proportional to the weight of that…
▽ More
We study growing networks in which each link carries a certain weight (randomly assigned at birth and fixed thereafter). The weight of a node is defined as the sum of the weights of the links attached to the node, and the network grows via the simplest weight-driven rule: A newly-added node is connected to an already existing node with the probability which is proportional to the weight of that node. We show that the node weight distribution n(w) has a universal, that is independent on the link weight distribution, tail: n(w) ~ w^-3 as w->oo. Results are particularly neat for the exponential link weight distribution when n(w) is algebraic over the entire weight range.
△ Less
Submitted 12 August, 2004;
originally announced August 2004.
-
Probability distribution of magnetization in the one-dimensional Ising model: Effects of boundary conditions
Authors:
T. Antal,
M. Droz,
Z. Racz
Abstract:
Finite-size scaling functions are investigated both for the mean-square magnetization fluctuations and for the probability distribution of the magnetization in the one-dimensional Ising model. The scaling functions are evaluated in the limit of the temperature going to zero (T -> 0), the size of the system going to infinity (N -> oo) while N[1-tanh(J/k_BT)] is kept finite (J being the nearest ne…
▽ More
Finite-size scaling functions are investigated both for the mean-square magnetization fluctuations and for the probability distribution of the magnetization in the one-dimensional Ising model. The scaling functions are evaluated in the limit of the temperature going to zero (T -> 0), the size of the system going to infinity (N -> oo) while N[1-tanh(J/k_BT)] is kept finite (J being the nearest neighbor coupling). Exact calculations using various boundary conditions (periodic, antiperiodic, free, block) demonstrate explicitly how the scaling functions depend on the boundary conditions. We also show that the block (small part of a large system) magnetization distribution results are identical to those obtained for free boundary conditions.
△ Less
Submitted 21 August, 2003;
originally announced August 2003.
-
Exponential velocity tails in a driven inelastic Maxwell model
Authors:
Tibor Antal,
Michel Droz,
Adam Lipowski
Abstract:
The problem of the steady-state velocity distribution in a driven inelastic Maxwell model of shaken granular material is revisited. Numerical solution of the master equation and analytical arguments show that the model has bilateral exponential velocity tails ($P(v)\sim e^{-|v|/\sqrt D}$), where $D$ is the amplitude of the noise. Previous study of this model predicted Gaussian tails (…
▽ More
The problem of the steady-state velocity distribution in a driven inelastic Maxwell model of shaken granular material is revisited. Numerical solution of the master equation and analytical arguments show that the model has bilateral exponential velocity tails ($P(v)\sim e^{-|v|/\sqrt D}$), where $D$ is the amplitude of the noise. Previous study of this model predicted Gaussian tails ($P(v)\sim e^{-av^2}$).
△ Less
Submitted 3 July, 2002;
originally announced July 2002.
-
On the Role of External Constraints in a Spatially Extended Evolutionary Prisoner's Dilemma Game
Authors:
Gyorgy Szabo,
Tibor Antal,
Peter Szabo,
Michel Droz
Abstract:
We study the emergency of mutual cooperation in evolutionary prisoner's dilemma games when the players are located on a square lattice. The players can choose one of the three strategies: cooperation (C), defection (D) or "tit for tat" (T), and their total payoffs come from games with the nearest neighbors. During the random sequential updates the players adopt one of their neighboring strategie…
▽ More
We study the emergency of mutual cooperation in evolutionary prisoner's dilemma games when the players are located on a square lattice. The players can choose one of the three strategies: cooperation (C), defection (D) or "tit for tat" (T), and their total payoffs come from games with the nearest neighbors. During the random sequential updates the players adopt one of their neighboring strategies if the chosen neighbor has higher payoff. We compare the effect of two types of external constraints added to the Darwinian evolutionary processes. In both cases the strategy of a randomly chosen player is replaced with probability P by another strategy. In the first case, the strategy is replaced by a randomly chosen one among the two others, while in the second case the new strategy is always C. Using generalized mean-field approximations and Monte Carlo simulations the strategy concentrations are evaluated in the stationary state for different strength of external constraints characterized by the probability P.
△ Less
Submitted 28 May, 2002;
originally announced May 2002.
-
Repeated bond traversal probabilities for the simple random walk
Authors:
T. Antal,
H. J. Hilhorst,
R. K. P. Zia
Abstract:
We consider the average number B_m(t) of bonds traversed exactly m times by a t step simple random walk. We determine B_m(t) explicitly in the scaling limit t -> oo with m/sqrt(t) fixed in dimension d=1 and m/log(t) fixed in dimension d=2. The scaling function is an erfc in d=1 and an exponential in d=2.
We consider the average number B_m(t) of bonds traversed exactly m times by a t step simple random walk. We determine B_m(t) explicitly in the scaling limit t -> oo with m/sqrt(t) fixed in dimension d=1 and m/log(t) fixed in dimension d=2. The scaling function is an erfc in d=1 and an exponential in d=2.
△ Less
Submitted 28 May, 2002;
originally announced May 2002.
-
Roughness distributions for 1/f^alpha signals
Authors:
T. Antal,
M. Droz,
G. Gyorgyi,
Z. Racz
Abstract:
The probability density function (PDF) of the roughness, i.e., of the temporal variance, of 1/f^alpha noise signals is studied. Our starting point is the generalization of the model of Gaussian, time-periodic, 1/f noise, discussed in our recent Letter [T. Antal et al., PRL, vol. 87, 240601 (2001)], to arbitrary power law. We investigate three main scaling regions, distinguished by the scaling of…
▽ More
The probability density function (PDF) of the roughness, i.e., of the temporal variance, of 1/f^alpha noise signals is studied. Our starting point is the generalization of the model of Gaussian, time-periodic, 1/f noise, discussed in our recent Letter [T. Antal et al., PRL, vol. 87, 240601 (2001)], to arbitrary power law. We investigate three main scaling regions, distinguished by the scaling of the cumulants in terms of the microscopic scale and the total length of the period. Various analytical representations of the PDF allow for a precise numerical evaluation of the scaling function of the PDF for any alpha. A simulation of the periodic process makes it possible to study also non-periodic signals on short intervals embedded in the full period. We find that for alpha=<1/2 the scaled PDF-s in both the periodic and the non-periodic cases are Gaussian, but for alpha>1/2 they differ from the Gaussian and from each other. Both deviations increase with growing alpha. That conclusion, based on numerics, is reinforced by analytic results for alpha=2 and alpha->infinity. We suggest that our theoretical and numerical results open a new perspective on the data analysis of 1/f^alpha processes.
△ Less
Submitted 15 December, 2001; v1 submitted 14 December, 2001;
originally announced December 2001.
-
1/f Noise and Extreme Value Statistics
Authors:
T. Antal,
M. Droz,
G. Gyorgyi,
Z. Racz
Abstract:
We study the finite-size scaling of the roughness of signals in systems displaying Gaussian 1/f power spectra. It is found that one of the extreme value distributions (Gumbel distribution) emerges as the scaling function when the boundary conditions are periodic. We provide a realistic example of periodic 1/f noise, and demonstrate by simulations that the Gumbel distribution is a good approximat…
▽ More
We study the finite-size scaling of the roughness of signals in systems displaying Gaussian 1/f power spectra. It is found that one of the extreme value distributions (Gumbel distribution) emerges as the scaling function when the boundary conditions are periodic. We provide a realistic example of periodic 1/f noise, and demonstrate by simulations that the Gumbel distribution is a good approximation for the case of nonperiodic boundary conditions as well. Experiments on voltage fluctuations in GaAs films are analyzed and excellent agreement is found with the theory.
△ Less
Submitted 30 May, 2001;
originally announced May 2001.
-
On the critical behavior of a lattice prey-predator model
Authors:
Tibor Antal,
Michel Droz,
Adam Lipowski,
Geza Odor
Abstract:
The critical properties of a simple prey-predator model are revisited. For some values of the control parameters, the model exhibits a line of directed percolation like transitions to a single absorbing state. For other values of the control parameters one finds a second line of continuous transitions toward infinite number of absorbing states, and the corresponding steady-state exponents are me…
▽ More
The critical properties of a simple prey-predator model are revisited. For some values of the control parameters, the model exhibits a line of directed percolation like transitions to a single absorbing state. For other values of the control parameters one finds a second line of continuous transitions toward infinite number of absorbing states, and the corresponding steady-state exponents are mean-field like. The critical behavior of the special point T (bicritical point), where the two transition lines meet, belongs to a different universality class. The use of dynamical Monte-Carlo method shows that a particular strategy for preparing the initial state should be devised to correctly describe the physics of the system near the second transition line. Relationships with a forest fire model with immunization are also discussed.
△ Less
Submitted 11 April, 2001;
originally announced April 2001.
-
Formation of Liesegang patterns: Simulations using a kinetic Ising model
Authors:
T. Antal,
M. Droz,
J. Magnin,
A. Pekalski,
Z. Racz
Abstract:
A kinetic Ising model description of Liesegang phenomena is studied using Monte Carlo simulations. The model takes into account thermal fluctuations, contains noise in the chemical reactions, and its control parameters are experimentally accessible. We find that noisy, irregular precipitation takes place in dimension d=2 while, depending on the values of the control parameters, either irregular…
▽ More
A kinetic Ising model description of Liesegang phenomena is studied using Monte Carlo simulations. The model takes into account thermal fluctuations, contains noise in the chemical reactions, and its control parameters are experimentally accessible. We find that noisy, irregular precipitation takes place in dimension d=2 while, depending on the values of the control parameters, either irregular patterns or precipitation bands satisfying the regular spacing law emerge in d=3.
△ Less
Submitted 29 September, 2000;
originally announced September 2000.
-
Phase Transitions and Oscillations in a Lattice Prey-Predator Model
Authors:
Tibor Antal,
Michel Droz
Abstract:
A coarse grained description of a two-dimensional prey-predator system is given in terms of a 3-state lattice model containing two control parameters: the spreading rates of preys and predators. The properties of the model are investigated by dynamical mean-field approximations and extensive numerical simulations. It is shown that the stationary state phase diagram is divided into two phases: a…
▽ More
A coarse grained description of a two-dimensional prey-predator system is given in terms of a 3-state lattice model containing two control parameters: the spreading rates of preys and predators. The properties of the model are investigated by dynamical mean-field approximations and extensive numerical simulations. It is shown that the stationary state phase diagram is divided into two phases: a pure prey phase and a coexistence phase of preys and predators in which temporal and spatial oscillations can be present. The different type of phase transitions occuring at the boundary of the prey absorbing phase, as well as the crossover phenomena occuring between the oscillatory and non-oscillatory domains of the coexistence phase are studied. The importance of finite size effects are discussed and scaling relations between different quantities are established. Finally, physical arguments, based on the spatial structure of the model, are given to explain the underlying mechanism leading to oscillations.
△ Less
Submitted 28 September, 2000;
originally announced September 2000.
-
Asymmetric exclusion process with next-nearest-neighbor interaction: some comments on traffic flow and a nonequilibrium reentrance transition
Authors:
T. Antal,
G. M. Schütz
Abstract:
We study the steady-state behavior of a driven non-equilibrium lattice gas of hard-core particles with next-nearest-neighbor interaction. We calculate the exact stationary distribution of the periodic system and for a particular line in the phase diagram of the system with open boundaries where particles can enter and leave the system. For repulsive interactions the dynamics can be interpreted a…
▽ More
We study the steady-state behavior of a driven non-equilibrium lattice gas of hard-core particles with next-nearest-neighbor interaction. We calculate the exact stationary distribution of the periodic system and for a particular line in the phase diagram of the system with open boundaries where particles can enter and leave the system. For repulsive interactions the dynamics can be interpreted as a two-speed model for traffic flow. The exact stationary distribution of the periodic continuous-time system turns out to coincide with that of the asymmetric exclusion process (ASEP) with discrete-time parallel update. However, unlike in the (single-speed) ASEP, the exact flow diagram for the two-speed model resembles in some important features the flow diagram of real traffic. The stationary phase diagram of the open system obtained from Monte Carlo simulations can be understood in terms of a shock moving through the system and an overfeeding effect at the boundaries, thus confirming theoretical predictions of a recently developed general theory of boundary-induced phase transitions. In the case of attractive interaction we observe an unexpected reentrance transition due to boundary effects.
△ Less
Submitted 20 January, 2000;
originally announced January 2000.
-
Spatial evolutionary prisoner's dilemma game with three strategies and external constraints
Authors:
Gyorgy Szabo,
Tibor Antal,
Peter Szabo,
Michel Droz
Abstract:
The emergence of mutual cooperation is studied in a spatially extended evolutionary prisoner's dilemma game in which the players are located on the sites of cubic lattices for dimensions d=1, 2, and 3. Each player can choose one of the three strategies: cooperation (C), defection (D) or Tit for Tat (T). During the evolutionary process the randomly chosen players adopt one of their neighboring st…
▽ More
The emergence of mutual cooperation is studied in a spatially extended evolutionary prisoner's dilemma game in which the players are located on the sites of cubic lattices for dimensions d=1, 2, and 3. Each player can choose one of the three strategies: cooperation (C), defection (D) or Tit for Tat (T). During the evolutionary process the randomly chosen players adopt one of their neighboring strategies if the chosen neighbor has higher payoff. Morover, an external constraint imposes that the players always cooperate with probability p. The stationary state phase diagram are computed by both using generalized mean-field approximations and Monte Carlo simulations. Nonequilibrium second order phase transitions assosiated with the extinction of one of the possible strategies are found and the corresponding critical exponents belong to the directed percolation universality class. It is shown that forcing externally the collaboration does not always produce the desired result.
△ Less
Submitted 5 April, 2000; v1 submitted 25 November, 1999;
originally announced November 1999.
-
Formation of Liesegang patterns: A spinodal decomposition scenario
Authors:
T. Antal,
M. Droz,
J. Magnin,
Z. Rácz
Abstract:
Spinodal decomposition in the presence of a moving particle source is proposed as a mechanism for the formation of Liesegang bands. This mechanism yields a sequence of band positions x_n that obeys the spacing law x_n~Q(1+p)^n. The dependence of the parameters p and Q on the initial concentration of the reagents is determined and we find that the functional form of p is in agreement with the exp…
▽ More
Spinodal decomposition in the presence of a moving particle source is proposed as a mechanism for the formation of Liesegang bands. This mechanism yields a sequence of band positions x_n that obeys the spacing law x_n~Q(1+p)^n. The dependence of the parameters p and Q on the initial concentration of the reagents is determined and we find that the functional form of p is in agreement with the experimentally observed Matalon-Packter law.
△ Less
Submitted 29 March, 1999;
originally announced March 1999.
-
Transport in the XX chain at zero temperature: Emergence of flat magnetization profiles
Authors:
T. Antal,
Z. Rácz,
A. Rákos,
G. M. Schütz
Abstract:
We study the connection between magnetization transport and magnetization profiles in zero-temperature XX chains. The time evolution of the transverse magnetization, m(x,t), is calculated using an inhomogeneous initial state that is the ground state at fixed magnetization but with m reversed from -m_0 for x<0 to m_0 for x>0. In the long-time limit, the magnetization evolves into a scaling form m…
▽ More
We study the connection between magnetization transport and magnetization profiles in zero-temperature XX chains. The time evolution of the transverse magnetization, m(x,t), is calculated using an inhomogeneous initial state that is the ground state at fixed magnetization but with m reversed from -m_0 for x<0 to m_0 for x>0. In the long-time limit, the magnetization evolves into a scaling form m(x,t)=P(x/t) and the profile develops a flat part (m=P=0) in the |x/t|< c(m_0) region. The flat region shrinks to zero if m_0->1/2 while it expands with the maximum velocity, c_0=1, for m_0->0. The states emerging in the scaling limit are compared to those of a homogeneous system where the same magnetization current is driven by a bulk field, and we find that the expectation values of various quantities (energy, occupation number in the fermionic representation) agree in the two systems.
△ Less
Submitted 14 December, 1998;
originally announced December 1998.
-
Derivation of the Matalon-Packter law for Liesegang patterns
Authors:
T. Antal,
M. Droz,
J. Magnin,
Z. Racz,
M. Zrinyi
Abstract:
Theoretical models of the Liesegang phenomena are studied and simple expressions for the spacing coefficients characterizing the patterns are derived. The emphasis is on displaying the explicit dependences on the concentrations of the inner- and the outer-electrolytes. Competing theories (ion-product supersaturation, nucleation and droplet growth, induced sol- coagulation) are treated with the a…
▽ More
Theoretical models of the Liesegang phenomena are studied and simple expressions for the spacing coefficients characterizing the patterns are derived. The emphasis is on displaying the explicit dependences on the concentrations of the inner- and the outer-electrolytes. Competing theories (ion-product supersaturation, nucleation and droplet growth, induced sol- coagulation) are treated with the aim of finding the distinguishing features of the theories. The predictions are compared with experiments and the results suggest that the induced sol-coagulation theory is the best candidate for describing the experimental observations embodied in the Matalon-Packter law.
△ Less
Submitted 17 July, 1998;
originally announced July 1998.
-
Isotropic Transverse XY Chain with Energy- and Magnetization Currents
Authors:
T. Antal,
Z. Racz,
A. Rakos,
G. M. Schutz
Abstract:
The ground-state correlations are investigated for an isotropic transverse XY chain which is constrained to carry either a current of magnetization J_M or a current of energy J_E. We find that the effect of nonzero J_M on the large-distance decay of correlations is twofold: i) oscillations are introduced and ii) the amplitude of the power law decay increases with increasing current. The effect o…
▽ More
The ground-state correlations are investigated for an isotropic transverse XY chain which is constrained to carry either a current of magnetization J_M or a current of energy J_E. We find that the effect of nonzero J_M on the large-distance decay of correlations is twofold: i) oscillations are introduced and ii) the amplitude of the power law decay increases with increasing current. The effect of energy current is more complex. Generically, correlations in current carrying states are found to decay faster than in the J_E=0 states, contrary to expectations that correlations are increased by the presence of currents. However, increasing the current, one reaches a special line where the correlations become comparable to those of the J_E=0 states. On this line, the symmetry of the ground state is enhanced and the transverse magnetization vanishes. Further increase of the current destroys the extra symmetry but the transverse magnetization remains at the high-symmetry, zero value.
△ Less
Submitted 12 November, 1997;
originally announced November 1997.
-
Nonequilibrium Steady State in a Quantum System: One-dimensional Transverse Ising Model with Energy Current
Authors:
T. Antal,
Z. Racz,
L. Sasvari
Abstract:
We study the nonequilibrium steady states of an Ising chain in a transverse field, h, by investigating the effect of a field, lambda, which drives the current of energy. The zero-temperature, h-lambda phase diagram is determined exactly and it is found that the energy current appears continuously above a threshold, lambda > lambda_c(h). The long-range magnetic order (present for h<h_c, lambda <…
▽ More
We study the nonequilibrium steady states of an Ising chain in a transverse field, h, by investigating the effect of a field, lambda, which drives the current of energy. The zero-temperature, h-lambda phase diagram is determined exactly and it is found that the energy current appears continuously above a threshold, lambda > lambda_c(h). The long-range magnetic order (present for h<h_c, lambda < lambda_c) is destroyed in the current-carrying state where the correlations are characterized by power-law, oscillatory decay. The mechanism which generates the power-law correlations in the current-carrying state is discussed.
△ Less
Submitted 4 December, 1996;
originally announced December 1996.
-
Critical domain size in a driven diffusive system
Authors:
Attila Szolnoki,
Tibor Antal,
Gyorgy Szabo
Abstract:
The homogeneous ordered state transforms into a polydomain state via a nucleation mechanism in two-dimensional lattice gas if the particle jumps are biased by an external field $E$. A simple phenomenological model is used to describe the time evolution of a circular interface separating the ordered regions. It is shown that the area of a domain increases if its radius exceeds a critical value pr…
▽ More
The homogeneous ordered state transforms into a polydomain state via a nucleation mechanism in two-dimensional lattice gas if the particle jumps are biased by an external field $E$. A simple phenomenological model is used to describe the time evolution of a circular interface separating the ordered regions. It is shown that the area of a domain increases if its radius exceeds a critical value proportional to $1/E$ which agrees qualitatively with Monte Carlo simulations.
△ Less
Submitted 7 March, 1996;
originally announced March 1996.
-
Dynamic Scaling of Width Distribution in Edwards--Wilkinson Type Models of Interface Dynamics
Authors:
T. Antal,
Z. Racz
Abstract:
Edwards--Wilkinson type models are studied in 1+1 dimensions and the time-dependent distribution, P_L(w^2,t), of the square of the width of an interface, w^2, is calculated for systems of size L. We find that, using a flat interface as an initial condition, P_L(w^2,t) can be calculated exactly and it obeys scaling in the form <w^2>_\infty P_L(w^2,t) = Phi(w^2 / <w^2>_\infty, t/L^2) where <w^2>_\…
▽ More
Edwards--Wilkinson type models are studied in 1+1 dimensions and the time-dependent distribution, P_L(w^2,t), of the square of the width of an interface, w^2, is calculated for systems of size L. We find that, using a flat interface as an initial condition, P_L(w^2,t) can be calculated exactly and it obeys scaling in the form <w^2>_\infty P_L(w^2,t) = Phi(w^2 / <w^2>_\infty, t/L^2) where <w^2>_\infty is the stationary value of w^2. For more complicated initial states, scaling is observed only in the large- time limit and the scaling function depends on the initial amplitude of the longest wavelength mode. The short-time limit is also interesting since P_L(w^2,t) is found to closely approximate the log-normal distribution. These results are confirmed by Monte Carlo simulations on a `roof-top' model of surface evolution.
△ Less
Submitted 31 October, 1995;
originally announced October 1995.
-
Self-organizing domain structure in a driven lattice gas
Authors:
Gyorgy Szabo,
Attila Szolnoki,
Tibor Antal,
Istvan Borsos
Abstract:
Instead of the homogeneous ordered particle distributions characteristic to equilibrium systems a self-organizing polydomain structure is found to be stable at low temperatures in a square lattice-gas model with repulsive nearest neighbor interaction when the particle jumps are biased by a uniform electric field. We have performed Monte Carlo simulations to study the domain structure varying the…
▽ More
Instead of the homogeneous ordered particle distributions characteristic to equilibrium systems a self-organizing polydomain structure is found to be stable at low temperatures in a square lattice-gas model with repulsive nearest neighbor interaction when the particle jumps are biased by a uniform electric field. We have performed Monte Carlo simulations to study the domain structure varying the temperatures, fields and system sizes. Revision of the traditional picture on this system raises interesting problems as discussed.
△ Less
Submitted 20 September, 1995;
originally announced September 1995.