Source author record

Joachim Krug

Joachim Krug appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

39works
17topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

39 published item(s)

preprint2021arXiv

From genotypes to organisms: State-of-the-art and perspectives of a cornerstone in evolutionary dynamics

Understanding how genotypes map onto phenotypes, fitness, and eventually organisms is arguably the next major missing piece in a fully predictive theory of evolution. We refer to this generally as the problem of the genotype-phenotype map. Though we are still far from achieving a complete picture of these relationships, our current understanding of simpler questions, such as the structure induced in the space of genotypes by sequences mapped to molecular structures, has revealed important facts that deeply affect the dynamical description of evolutionary processes. Empirical evidence supporting the fundamental relevance of features such as phenotypic bias is mounting as well, while the synthesis of conceptual and experimental progress leads to questioning current assumptions on the nature of evolutionary dynamics-cancer progression models or synthetic biology approaches being notable examples. This work delves into a critical and constructive attitude in our current knowledge of how genotypes map onto molecular phenotypes and organismal functions, and discusses theoretical and empirical avenues to broaden and improve this comprehension. As a final goal, this community should aim at deriving an updated picture of evolutionary processes soundly relying on the structural properties of genotype spaces, as revealed by modern techniques of molecular and functional analysis.

preprint2021arXiv

Non-equilibrium Roughness Evolution of Small Molecule Mixed Films Reflecting Equilibrium Phase Behavior

Understanding non-equilibrium phenomena, such as growth, and connecting them to equilibrium phase behavior is a major challenge, in particular for complex multicomponent materials. We use X-ray reflectivity to determine the surface roughness of binary mixtures of several prototypical organic compounds. By analyzing the roughness as a function of composition, we find a systematic behavior depending on the bulk phase behavior in terms of intermixing, co-crystallization or phase separation. Supported by kinetic Monte Carlo simulations, we provide evidence that the growth behavior can be rationalized by a lowered step edge barrier in the mixed films which is induced by reduced in-plane crystallinity.

preprint2020arXiv

Distribution of the number of fitness maxima in Fisher's Geometric Model

Fisher's geometric model describes biological fitness landscapes by combining a linear map from the discrete space of genotypes to an $n$-dimensional Euclidean phenotype space with a nonlinear, single-peaked phenotype-fitness map. Genotypes are represented by binary sequences of length $L$, and the phenotypic effects of mutations at different sites are represented by $L$ random vectors drawn from an isotropic Gaussian distribution. Recent work has shown that the interplay between the genotypic and phenotypic levels gives rise to a range of different landscape topographies that can be characterised by the number of local fitness maxima. Extending our previous study of the mean number of local maxima, here we focus on the distribution of the number of maxima when the limit $L \to \infty$ is taken at finite $n$. We identify the typical scale of the number of maxima for general $n$, and determine the full scaled probability density and two point correlation function of maxima for the one-dimensional case. We also elaborate on the close relation of the model to the anti-ferromagnetic Hopfield model with $n$ random continuous pattern vectors, and show that many of our results carry over to this setting. More generally, we expect that our analysis can help to elucidate the fluctuation structure of metastable states in various spin glass problems.

preprint2016arXiv

$δ$-exceedance records and random adaptive walks

We study a modified record process where the $k$'th record in a series of independent and identically distributed random variables is defined recursively through the condition $Y_k > Y_{k-1} - δ_{k-1}$ with a deterministic sequence $δ_k > 0$ called the handicap. For constant $δ_k \equiv δ$ and exponentially distributed random variables it has been shown in previous work that the process displays a phase transition as a function of $δ$ between a normal phase where the mean record value increases indefinitely and a stationary phase where the mean record value remains bounded and a finite fraction of all entries are records (Park \textit{et al} 2015 {\it Phys. Rev.} E \textbf{91} 042707). Here we explore the behavior for general probability distributions and decreasing and increasing sequences $δ_k$, focusing in particular on the case when $δ_k$ matches the typical spacing between subsequent records in the underlying simple record process without handicap. We find that a continuous phase transition occurs only in the exponential case, but a novel kind of first order transition emerges when $δ_k$ is increasing. The problem is partly motivated by the dynamics of evolutionary adaptation in biological fitness landscapes, where $δ_k$ corresponds to the change of the deterministic fitness component after $k$ mutational steps. The results for the record process are used to compute the mean number of steps that a population performs in such a landscape before being trapped at a local fitness maximum.

preprint2016arXiv

Greedy adaptive walks on a correlated fitness landscape

We study adaptation of a haploid asexual population on a fitness landscape defined over binary genotype sequences of length $L$. We consider greedy adaptive walks in which the population moves to the fittest among all single mutant neighbors of the current genotype until a local fitness maximum is reached. The landscape is of the rough mount Fuji type, which means that the fitness value assigned to a sequence is the sum of a random and a deterministic component. The random components are independent and identically distributed random variables, and the deterministic component varies linearly with the distance to a reference sequence. The deterministic fitness gradient $c$ is a parameter that interpolates between the limits of an uncorrelated random landscape ($c = 0$) and an effectively additive landscape ($c \to \infty$). When the random fitness component is chosen from the Gumbel distribution, explicit expressions for the distribution of the number of steps taken by the greedy walk are obtained, and it is shown that the walk length varies non-monotonically with the strength of the fitness gradient when the starting point is sufficiently close to the reference sequence. Asymptotic results for general distributions of the random fitness component are obtained using extreme value theory, and it is found that the walk length attains a non-trivial limit for $L \to \infty$, different from its values for $c=0$ and $c = \infty$, if $c$ is scaled with $L$ in an appropriate combination.

preprint2015arXiv

Analysis of adaptive walks on NK fitness landscapes with different interaction schemes

Fitness landscapes are genotype to fitness mappings commonly used in evolutionary biology and computer science which are closely related to spin glass models. In this paper, we study the NK model for fitness landscapes where the interaction scheme between genes can be explicitly defined. The focus is on how this scheme influences the overall shape of the landscape. Our main tool for the analysis are adaptive walks, an idealized dynamics by which the population moves uphill in fitness and terminates at a local fitness maximum. We use three different types of walks and investigate how their length (the number of steps required to reach a local peak) and height (the fitness at the endpoint of the walk) depend on the dimensionality and structure of the landscape. We find that the distribution of local maxima over the landscape is particularly sensitive to the choice of interaction pattern. Most quantities that we measure are simply correlated to the rank of the scheme, which is equal to the number of nonzero coefficients in the expansion of the fitness landscape in terms of Walsh functions.

preprint2015arXiv

Phase transition in random adaptive walks on correlated fitness landscapes

We study biological evolution on a random fitness landscape where correlations are introduced through a linear fitness gradient of strength $c$. When selection is strong and mutations rare the dynamics is a directed uphill walk that terminates at a local fitness maximum. We analytically calculate the dependence of the walk length on the genome size $L$. When the distribution of the random fitness component has an exponential tail we find a phase transition of the walk length $D$ between a phase at small $c$ where walks are short $(D \sim \ln L)$ and a phase at large $c$ where walks are long $(D \sim L)$. For all other distributions only a single phase exists for any $c > 0$. The considered process is equivalent to a zero temperature Metropolis dynamics for the random energy model in an external magnetic field, thus also providing insight into the aging dynamics of spin glasses.

preprint2015arXiv

Symmetry-based coarse-graining of evolved dynamical networks

Networks with a prescribed power-law scaling in the spectrum of the graph Laplacian can be generated by evolutionary optimization. The Laplacian spectrum encodes the dynamical behavior of many important processes. Here, the networks are evolved to exhibit subdiffusive dynamics. Under the additional constraint of degree-regularity, the evolved networks display an abundance of symmetric motifs arranged into loops and long linear segments. Exploiting results from algebraic graph theory on symmetric networks, we find the underlying backbone structures and how they contribute to the spectrum. The resulting coarse-grained networks provide an intuitive view of how the anomalous diffusive properties can be realized in the evolved structures.

preprint2014arXiv

Adaptation in tunably rugged fitness landscapes: The Rough Mount Fuji Model

Much of the current theory of adaptation is based on Gillespie's mutational landscape model (MLM), which assumes that the fitness values of genotypes linked by single mutational steps are independent random variables. On the other hand, a growing body of empirical evidence shows that real fitness landscapes, while possessing a considerable amount of ruggedness, are smoother than predicted by the MLM. In the present article we propose and analyse a simple fitness landscape model with tunable ruggedness based on the Rough Mount Fuji (RMF) model originally introduced by Aita et al. [Biopolymers 54:64-79 (2000)] in the context of protein evolution. We provide a comprehensive collection of results pertaining to the topographical structure of RMF landscapes, including explicit formulae for the expected number of local fitness maxima, the location of the global peak, and the fitness correlation function. The statistics of single and multiple adaptive steps on the RMF landscape are explored mainly through simulations, and the results are compared to the known behavior in the MLM model. Finally, we show that the RMF model can explain the large number of second-step mutations observed on a highly-fit first step backgound in a recent evolution experiment with a microvirid bacteriophage [Miller et al., Genetics 187:185-202 (2011)].

preprint2014arXiv

Clonal interference and Muller's ratchet in spatial habitats

Competition between independently arising beneficial mutations is enhanced in spatial populations due to the linear rather than exponential growth of clones. Recent theoretical studies have pointed out that the resulting fitness dynamics is analogous to a surface growth process, where new layers nucleate and spread stochastically, leading to the build up of scale-invariant roughness. This scenario differs qualitatively from the standard view of adaptation in that the speed of adaptation becomes independent of population size while the fitness variance does not. Here we exploit recent progress in the understanding of surface growth processes to obtain precise predictions for the universal, non-Gaussian shape of the fitness distribution for one-dimensional habitats, which are verified by simulations. When the mutations are deleterious rather than beneficial the problem becomes a spatial version of Muller's ratchet. In contrast to the case of well-mixed populations, the rate of fitness decline remains finite even in the limit of an infinite habitat, provided the ratio $U_d/s^2$ between the deleterious mutation rate and the square of the (negative) selection coefficient is sufficiently large. Using again an analogy to surface growth models we show that the transition between the stationary and the moving state of the ratchet is governed by directed percolation.

preprint2014arXiv

Evolutionary accessibility of modular fitness landscapes

A fitness landscape is a mapping from the space of genetic sequences, which is modeled here as a binary hypercube of dimension $L$, to the real numbers. We consider random models of fitness landscapes, where fitness values are assigned according to some probabilistic rule, and study the statistical properties of pathways to the global fitness maximum along which fitness increases monotonically. Such paths are important for evolution because they are the only ones that are accessible to an adapting population when mutations occur at a low rate. The focus of this work is on the block model introduced by A.S. Perelson and C.A. Macken [Proc. Natl. Acad. Sci. USA 92:9657 (1995)] where the genome is decomposed into disjoint sets of loci (`modules') that contribute independently to fitness, and fitness values within blocks are assigned at random. We show that the number of accessible paths can be written as a product of the path numbers within the blocks, which provides a detailed analytic description of the path statistics. The block model can be viewed as a special case of Kauffman's NK-model, and we compare the analytic results to simulations of the NK-model with different genetic architectures. We find that the mean number of accessible paths in the different versions of the model are quite similar, but the distribution of the path number is qualitatively different in the block model due to its multiplicative structure. A similar statement applies to the number of local fitness maxima in the NK-models, which has been studied extensively in previous works. The overall evolutionary accessibility of the landscape, as quantified by the probability to find at least one accessible path to the global maximum, is dramatically lowered by the modular structure.

preprint2014arXiv

Multidimensional epistasis and the transitory advantage of sex

Identifying and quantifying the benefits of sex and recombination is a long standing problem in evolutionary theory. In particular, contradictory claims have been made about the existence of a benefit of recombination on high dimensional fitness landscapes in the presence of sign epistasis. Here we present a comparative numerical study of sexual and asexual evolutionary dynamics of haploids on tunably rugged model landscapes under strong selection, paying special attention to the temporal development of the evolutionary advantage of recombination and the link between population diversity and the rate of adaptation. We show that the adaptive advantage of recombination on static rugged landscapes is strictly transitory. At early times, an advantage of recombination arises through the possibility to combine individually occurring beneficial mutations, but this effect is reversed at longer times by the much more efficient trapping of recombining populations at local fitness peaks. These findings are explained by means of well established results for a setup with only two loci. In accordance with the Red Queen hypothesis the transitory advantage can be prolonged indefinitely in fluctuating environments, and it is maximal when the environment fluctuates on the same time scale on which trapping at local optima typically occurs.

preprint2013arXiv

Exact Results for Amplitude Spectra of Fitness Landscapes

Starting from fitness correlation functions, we calculate exact expressions for the amplitude spectra of fitness landscapes as defined by P.F. Stadler [J. Math. Chem. 20, 1 (1996)] for common landscape models, including Kauffman's NK-model, rough Mount Fuji landscapes and general linear superpositions of such landscapes. We further show that correlations decaying exponentially with the Hamming distance yield exponentially decaying spectra similar to those reported recently for a model of molecular signal transduction. Finally, we compare our results for the model systems to the spectra of various experimentally measured fitness landscapes. We claim that our analytical results should be helpful when trying to interpret empirical data and guide the search for improved fitness landscape models.

preprint2013arXiv

Rate of adaptation in sexuals and asexuals: A solvable model of the Fisher-Muller effect

The adaptation of large asexual populations is hampered by the competition between independently arising beneficial mutations in different individuals, which is known as clonal interference. Fisher and Muller proposed that recombination provides an evolutionary advantage in large populations by alleviating this competition. Based on recent progress in quantifying the speed of adaptation in asexual populations undergoing clonal interference, we present a detailed analysis of the Fisher-Muller mechanism for a model genome consisting of two loci with an infinite number of beneficial alleles each and multiplicative fitness effects. We solve the infinite population dynamics exactly and show that, for a particular, natural mutation scheme, the speed of adaptation in sexuals is twice as large as in asexuals. Guided by the infinite population result and by previous work on asexual adaptation, we postulate an expression for the speed of adaptation in finite sexual populations that agrees with numerical simulations over a wide range of population sizes and recombination rates. The ratio of the sexual to asexual adaptation speed is a function of population size that increases in the clonal interference regime and approaches 2 for extremely large populations. The simulations also show that the imbalance between the numbers of accumulated mutations at the two loci is strongly suppressed even by a small amount of recombination. The generalization of the model to an arbitrary number $L$ of loci is briefly discussed. If each offspring samples the alleles at each locus from the gene pool of the whole population rather than from two parents, the ratio of the sexual to asexual adaptation speed is approximately equal to $L$ in large populations. A possible realization of this scenario is the reassortment of genetic material in RNA viruses with $L$ genomic segments.

preprint2013arXiv

Record occurrence and record values in daily and monthly temperatures

We analyze the occurrence and the values of record-breaking temperatures in daily and monthly temperature observations. Our aim is to better understand and quantify the statistics of temperature records in the context of global warming. Similar to earlier work we employ a simple mathematical model of independent and identically distributed random variables with a linearly growing expectation value. This model proved to be useful in predicting the increase (decrease) in upper (lower) temperature records in a warming climate. Using both station and re-analysis data from Europe and the United States we further investigate the statistics of temperature records and the validity of this model. The most important new contribution in this article is an analysis of the statistics of record values for our simple model and European reanalysis data. We estimate how much the mean values and the distributions of record temperatures are affected by the large scale warming trend. In this context we consider both the values of records that occur at a certain time and the values of records that have a certain record number in the series of record events. We compare the observational data both to simple analytical computations and numerical simulations. We find that it is more difficult to describe the values of record breaking temperatures within the framework of our linear drift model. Observations from the summer months fit well into the model with Gaussian random variables under the observed linear warming, in the sense that record breaking temperatures are more extreme in the summer. In winter however a significant asymmetry of the daily temperature distribution hides the effect of the slow warming trends. Therefore very extreme cold records are still possible in winter. This effect is even more pronounced if one considers only data from subpolar regions.

preprint2012arXiv

Evolutionary accessibility in tunably rugged fitness landscapes

The adaptive evolution of a population under the influence of mutation and selection is strongly influenced by the structure of the underlying fitness landscape, which encodes the interactions between mutations at different genetic loci. Theoretical studies of such landscapes have been carried out for several decades, but only recently experimental fitness measurements encompassing all possible combinations of small sets of mutations have become available. The empirical studies have spawned new questions about the accessibility of optimal genotypes under natural selection. Depending on population dynamic parameters such as mutation rate and population size, evolutionary accessibility can be quantified through the statistics of accessible mutational pathways (along which fitness increases monotoncially), or through the study of the basin of attraction of the optimal genotype under greedy (steepest ascent) dynamics. Here we investigate these two measures of accessibility in the framework of Kauffman's $LK$-model, a paradigmatic family of random fitness landscapes with tunable ruggedness. The key parameter governing the strength of genetic interactions is the number $K$ of interaction partners of each of the $L$ sites in the genotype sequence. In general, accessibility increases with increasing genotype dimensionality $L$ and decreases with increasing number of interactions $K$. Remarkably, however, we find that some measures of accessibility behave non-monotonically as a function of $K$, indicating a special role of the most sparsely connected, non-trivial cases K=1 and 2. The relation between models for fitness landscapes and spin glasses is also addressed.

preprint2012arXiv

Non-conserved dynamics of steps on vicinal surfaces during electromigration-induced step bunching

We report new results on the non-conserved dynamics of parallel steps on vicinal surfaces in the case of sublimation with electromigration and step-step interactions. The derived equations are valid in the quasistatic approximation and in the limit $f^{-1}\gg l_D\gg l_{\pm} \gg l_i$, where $f$ is the inverse electromigration length, $l_D$ the diffusion length, $l_{\pm}$ the kinetic lengths and $l_i$ the terrace widths. The coupling between crystal sublimation and step-step interactions induces non-linear, non-conservative terms in the equations of motion. Depending on the initial conditions, this leads to interrupted coarsening, anticoarsening of step bunches or periodic switching between step trains of different numbers of bunches.

preprint2012arXiv

Quantitative analyses of empirical fitness landscapes

The concept of a fitness landscape is a powerful metaphor that offers insight into various aspects of evolutionary processes and guidance for the study of evolution. Until recently, empirical evidence on the ruggedness of these landscapes was lacking, but since it became feasible to construct all possible genotypes containing combinations of a limited set of mutations, the number of studies has grown to a point where a classification of landscapes becomes possible. The aim of this review is to identify measures of epistasis that allow a meaningful comparison of fitness landscapes and then apply them to the empirical landscapes to discern factors that affect ruggedness. The various measures of epistasis that have been proposed in the literature appear to be equivalent. Our comparison shows that the ruggedness of the empirical landscape is affected by whether the included mutations are beneficial or deleterious and by whether intra- or intergenic epistasis is involved. Finally, the empirical landscapes are compared to landscapes generated with the Rough Mt.\ Fuji model. Despite the simplicity of this model, it captures the features of the experimental landscapes remarkably well.

preprint2011arXiv

Adaptive walks and extreme value theory

We study biological evolution in a high-dimensional genotype space in the regime of rare mutations and strong selection. The population performs an uphill walk which terminates at local fitness maxima. Assigning fitness randomly to genotypes, we show that the mean walk length is logarithmic in the number of initially available beneficial mutations, with a prefactor determined by the tail of the fitness distribution. This result is derived analytically in a simplified setting where the mutational neighborhood is fixed during the adaptive process, and confirmed by numerical simulations.

preprint2011arXiv

Correlations between record events in sequences of random variables with a linear trend

The statistics of records in sequences of independent, identically distributed random variables is a classic subject of study. One of the earliest results concerns the stochastic independence of record events. Recently, records statistics beyond the case of i.i.d. random variables have received much attention, but the question of independence of record events has not been addressed systematically. In this paper, we study this question in detail for the case of independent, non-identically distributed random variables, specifically, for random variables with a linearly moving mean. We find a rich pattern of positive and negative correlations, and show how their asymptotics is determined by the universality classes of extreme value statistics.

preprint2011arXiv

Evolutionary accessibility of mutational pathways

Functional effects of different mutations are known to combine to the total effect in highly nontrivial ways. For the trait under evolutionary selection (`fitness'), measured values over all possible combinations of a set of mutations yield a fitness landscape that determines which mutational states can be reached from a given initial genotype. Understanding the accessibility properties of fitness landscapes is conceptually important in answering questions about the predictability and repeatability of evolutionary adaptation. Here we theoretically investigate accessibility of the globally optimal state on a wide variety of model landscapes, including landscapes with tunable ruggedness as well as neutral `holey' landscapes. We define a mutational pathway to be accessible if it contains the minimal number of mutations required to reach the target genotype, and if fitness increases in each mutational step. Under this definition accessibility is high, in the sense that at least one accessible pathwayexists with a substantial probability that approaches unity as the dimensionality of the fitness landscape (set by the number of mutational loci) becomes large. At the same time the number of alternative accessible pathways grows without bound. We test the model predictions against an empirical 8-locus fitness landscape obtained for the filamentous fungus \textit{Aspergillus niger}. By analyzing subgraphs of the full landscape containing different subsets of mutations, we are able to probe the mutational distance scale in the empirical data. The predicted effect of high accessibility is supported by the empirical data and very robust, which we argue to reflect the generic topology of sequence spaces.

preprint2011arXiv

Evolutionary advantage of small populations on complex fitness landscapes

Background: Recent experimental and theoretical studies have shown that small asexual populations evolving on complex fitness landscapes may achieve a higher fitness than large ones due to the increased heterogeneity of adaptive trajectories. Here we introduce a class of haploid three-locus fitness landscapes that allows to investigate this scenario in a precise and quantitative way. Results: Our main result derived analytically shows how the probability of choosing the path of largest initial fitness increase grows with the population size. This makes large populations more likely to get trapped at local fitness peaks and implies an advantage of small populations at intermediate time scales. The range of population sizes where this effect is operative coincides with the onset of clonal interference. Additional studies using ensembles of random fitness landscapes show that the results achieved for a particular choice of three-locus landscape parameters are robust and also persist as the number of loci increases. Conclusions: Our study indicates that an advantage for small populations is likely whenever the fitness landscape contains local maxima. The advantage appears at intermediate time scales, which are long enough for trapping at local fitness maxima to have occurred but too short for peak escape by the creation of multiple mutants.

preprint2011arXiv

Rare events in population genetics: Stochastic tunneling in a two-locus model with recombination

We study the evolution of a population in a two-locus genotype space, in which the negative effects of two single mutations are overcompensated in a high fitness double mutant. We discuss how the interplay of finite population size, $N$, and sexual recombination at rate $r$ affects the escape times $t_\mathrm{esc}$ to the double mutant. For small populations demographic noise generates massive fluctuations in $t_\mathrm{esc}$. The mean escape time varies non-monotonically with $r$, and grows exponentially as $\ln t_{\mathrm{esc}} \sim N(r - r^\ast)^{3/2}$ beyond a critical value $r^\ast$.

preprint2011arXiv

Record statistics for biased random walks, with an application to financial data

We consider the occurrence of record-breaking events in random walks with asymmetric jump distributions. The statistics of records in symmetric random walks was previously analyzed by Majumdar and Ziff and is well understood. Unlike the case of symmetric jump distributions, in the asymmetric case the statistics of records depends on the choice of the jump distribution. We compute the record rate $P_n(c)$, defined as the probability for the $n$th value to be larger than all previous values, for a Gaussian jump distribution with standard deviation $σ$ that is shifted by a constant drift $c$. For small drift, in the sense of $c/σ\ll n^{-1/2}$, the correction to $P_n(c)$ grows proportional to arctan$(\sqrt{n})$ and saturates at the value $\frac{c}{\sqrt{2} σ}$. For large $n$ the record rate approaches a constant, which is approximately given by $1-(σ/\sqrt{2π}c)\textrm{exp}(-c^2/2σ^2)$ for $c/σ\gg 1$. These asymptotic results carry over to other continuous jump distributions with finite variance. As an application, we compare our analytical results to the record statistics of 366 daily stock prices from the Standard & Poors 500 index. The biased random walk accounts quantitatively for the increase in the number of upper records due to the overall trend in the stock prices, and after detrending the number of upper records is in good agreement with the symmetric random walk. However the number of lower records in the detrended data is significantly reduced by a mechanism that remains to be identified.

preprint2010arXiv

A pedestrian's view on interacting particle systems, KPZ universality, and random matrices

These notes are based on lectures delivered by the authors at a Langeoog seminar of SFB/TR12 "Symmetries and universality in mesoscopic systems" to a mixed audience of mathematicians and theoretical physicists. After a brief outline of the basic physical concepts of equilibrium and nonequilibrium states, the one-dimensional simple exclusion process is introduced as a paradigmatic nonequilibrium interacting particle system. The stationary measure on the ring is derived and the idea of the hydrodynamic limit is sketched. We then introduce the phenomenological Kardar-Parisi-Zhang (KPZ) equation and explain the associated universality conjecture for surface fluctuations in growth models. This is followed by a detailed exposition of a seminal paper of Johansson that relates the current fluctuations of the totally asymmetric simple exclusion process (TASEP) to the Tracy-Widom distribution of random matrix theory. The implications of this result are discussed within the framework of the KPZ conjecture.

preprint2010arXiv

Anti-Coarsening and Complex Dynamics of Step Bunches on Vicinal Surfaces during Sublimation

A sublimating vicinal crystal surface can undergo a step bunching instability when the attachment-detachment kinetics is asymmetric, in the sense of a normal Ehrlich-Schwoebel effect. Here we investigate this instability in a model that takes into account the subtle interplay between sublimation and step-step interactions, which breaks the volume-conserving character of the dynamics assumed in previous work. On the basis of a systematically derived continuum equation for the surface profile, we argue that the non-conservative terms pose a limitation on the size of emerging step bunches. This conclusion is supported by extensive simulations of the discrete step dynamics, which show breakup of large bunches into smaller ones as well as arrested coarsening and periodic oscillations between states with different numbers of bunches.

preprint2010arXiv

Bistability in two-locus models with selection, mutation, and recombination

The evolutionary effect of recombination depends crucially on the epistatic interactions between linked loci. A paradigmatic case where recombination is known to be strongly disadvantageous is a two-locus fitness landscape dis- playing reciprocal sign epistasis with two fitness peaks of unequal height. Focusing on the occurrence of bistability in the equilibrium solutions, we consider here the deterministic, haploid two-locus model with reversible mu- tations, selection and recombination. We find analytic formulae for the criti- cal recombination probability rc above which two stable stationary solutions appear which are localized on each of the two fitness peaks. We also derive the stationary genotype frequencies in various parameter regimes. When the recombination rate is close to rc and the fitness difference between the two peaks is small, we obtain a compact description in terms of a cubic polyno- mial which is analogous to the Landau theory of physical phase transitions.

preprint2010arXiv

Record-breaking temperatures reveal a warming climate

We present a mathematical analysis of records drawn from independent random variables with a drifting mean. To leading order the change in the record rate is proportional to the ratio of the drift velocity to the standard deviation of the underlying distribution. We apply the theory to time series of daily temperatures for given calendar days, obtained from historical climate recordings of European and American weather stations as well as re-analysis data. We conclude that the change in the mean temperature has increased the rate of record breaking events in a moderate but significant way: For the European station data covering the time period 1976-2005, we find that about 5 of the 17 high temperature records observed on average in 2005 can be attributed to the warming climate.

preprint2010arXiv

Records and sequences of records from random variables with a linear trend

We consider records and sequences of records drawn from discrete time series of the form $X_{n}=Y_{n}+cn$, where the $Y_{n}$ are independent and identically distributed random variables and $c$ is a constant drift. For very small and very large drift velocities, we investigate the asymptotic behavior of the probability $p_n(c)$ of a record occurring in the $n$th step and the probability $P_N(c)$ that all $N$ entries are records, i.e. that $X_1 < X_2 < ... < X_N$. Our work is motivated by the analysis of temperature time series in climatology, and by the study of mutational pathways in evolutionary biology.

preprint2009arXiv

Accurate rate coefficients for models of interstellar gas-grain chemistry

The methodology for modeling grain-surface chemistry has been greatly improved by taking into account the grain size and fluctuation effects. However, the reaction rate coefficients currently used in all practical models of gas-grain chemistry are inaccurate by a significant amount. We provide expressions for these crucial rate coefficients that are both accurate and easy to incorporate into gas-grain models. We use exact results obtained in earlier work, where the reaction rate coefficient was defined by a first-passage problem, which was solved using random walk theory. The approximate reaction rate coefficient presented here is easy to include in all models of interstellar gas-grain chemistry. In contrast to the commonly used expression, the results that it provides are in perfect agreement with detailed kinetic Monte Carlo simulations. We also show the rate coefficient for reactions involving multiple species.

preprint2009arXiv

Robustness and epistasis in mutation-selection models

We investigate the fitness advantage associated with the robustness of a phenotype against deleterious mutations using deterministic mutation-selection models of quasispecies type equipped with a mesa shaped fitness landscape. We obtain analytic results for the robustness effect which become exact in the limit of infinite sequence length. Thereby, we are able to clarify a seeming contradiction between recent rigorous work and an earlier heuristic treatment based on a mapping to a Schrödinger equation. We exploit the quantum mechanical analogy to calculate a correction term for finite sequence lengths and verify our analytic results by numerical studies. In addition, we investigate the occurrence of an error threshold for a general class of epistatic landscape and show that diminishing epistasis is a necessary but not sufficient condition for error threshold behavior.

preprint2009arXiv

The speed of evolution in large asexual populations

We consider an asexual biological population of constant size $N$ evolving in discrete time under the influence of selection and mutation. Beneficial mutations appear at rate $U$ and their selective effects $s$ are drawn from a distribution $g(s)$. After introducing the required models and concepts of mathematical population genetics, we review different approaches to computing the speed of logarithmic fitness increase as a function of $N$, $U$ and $g(s)$. We present an exact solution of the infinite population size limit and provide an estimate of the population size beyond which it is valid. We then discuss approximate approaches to the finite population problem, distinguishing between the case of a single selection coefficient, $g(s) = δ(s - s_b)$, and a continuous distribution of selection coefficients. Analytic estimates for the speed are compared to numerical simulations up to population sizes of order $10^{300}$.

preprint2007arXiv

Ehrlich-Schwoebel Effect for Organic Molecules: Direct Calculation of the Step Edge Barrier using Empirical Potentials

The step edge barrier of a prototypical organic semiconductor molecule, 3,4,9,10-perylene-tetracaboxylic-dianhydride (PTCDA) has been analysed by means of calculations based on emperical potentials. The minimum energy path (MEP) has been calculated for a single molecule on a substrate of three molecular layers between equivalent minimum energy positions within two neighboring unit cells. To determine the step edge barrier, we have calculated the MEP over a step to a fourth layer of molecules. We found energy barriers of E_D= 80 meV for in-layer diffusion and E_S = 750 meV for step crossing, indicating a strong Ehrlich-Schwoebel effect for PTCDA.

preprint2007arXiv

Morphological stability of electromigration-driven vacancy islands

The electromigration-induced shape evolution of two-dimensional vacancy islands on a crystal surface is studied using a continuum approach. We consider the regime where mass transport is restricted to terrace diffusion in the interior of the island. In the limit of fast attachment/detachment kinetics a circle translating at constant velocity is a stationary solution of the problem. In contrast to earlier work [O. Pierre-Louis and T.L. Einstein, Phys. Rev. B 62, 13697 (2000)] we show that the circular solution remains linearly stable for arbitrarily large driving forces. The numerical solution of the full nonlinear problem nevertheless reveals a fingering instability at the trailing end of the island, which develops from finite amplitude perturbations and eventually leads to pinch-off. Relaxing the condition of instantaneous attachment/detachment kinetics, we obtain non-circular elongated stationary shapes in an analytic approximation which compares favorably to the full numerical solution.

preprint2006arXiv

Kinetic Monte Carlo simulations of oscillatory shape evolution for electromigration-driven islands

The shape evolution of two-dimensional islands under electromigration-driven periphery diffusion is studied by kinetic Monte Carlo (KMC) simulations and continuum theory. The energetics of the KMC model is adapted to the Cu(100) surface, and the continuum model is matched to the KMC model by a suitably parametrized choice of the orientation-dependent step stiffness and step atom mobility. At 700 K shape oscillations predicted by continuum theory are quantitatively verified by the KMC simulations, while at 500 K qualitative differences between the two modeling approaches are found.

preprint2001arXiv

Coarsening of vortex ripples in sand

The coarsening of an array of vortex ripples prepared in an unstable state is discussed within the framework of a simple mass transfer model first introduced by K.H. Andersen et al. [Phys. Rev. E 63, 066308 (2001)]. Two scenarios for the selection of the final pattern are identified. When the initial state is homogeneous with uniform random perturbations, a unique final state is reached which depends only on the shape of the interaction function $f(λ)$. A potential formulation of the dynamics suggests that the final wavelength is determined by a Maxwell construction applied to $f(λ)$, but comparison with numerical simulations shows that this yields only an upper bound. In contrast, the evolution from a perfectly homogeneous state with a localized perturbation proceeds through the propagation of wavelength doubling fronts. The front speed can be predicted by standard marginal stability theory. In this case the final wavelength depends on the initial wavelength in a complicated manner which involves multiplication by factors of 2 and rational ratios such as 4/3.

preprint1999arXiv

Phase separation in disordered exclusion models

The effect of quenched disorder in the one-dimensional asymmetric exclusion process is reviewed. Both particlewise and sitewise disorder generically induces phase separation in a range of densities. In the particlewise case the existence of stationary product measures in the homogeneous phase implies that the critical density can be computed exactly, while for sitewise disorder only bounds are available. The coarsening of phase-separated domains starting from a homogeneous initial condition is addressed using scaling arguments and extremal statistics considerations. Some of these results have been obtained previously in the context of directed polymers subject to columnar disorder.