BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Su-Chan Park

Su-Chan Park appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.stat-mechPopulations and Evolutioncond-mat.dis-nn

Catalog footprint

What is connected

22works
3topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

22 published item(s)

preprint2020arXiv

Branching annihilating random walks with long-range attraction in one dimension

We introduce and numerically study the branching annihilating random walks with long-range attraction (BAWL). The long-range attraction makes hopping biased in such a manner that particle's hopping along the direction to the nearest particle has larger transition rate than hopping against the direction. Still, unlike the Lévy flight, a particle only hops to one of its nearest-neighbor sites. The strength of bias takes the form $x^{-σ}$ with non-negative $σ$, where $x$ is the distance to the nearest particle from a particle to hop. By extensive Monte Carlo simulations, we show that the critical decay exponent $δ$ varies continuously with $σ$ up to $σ=1$ and $δ$ is the same as the critical decay exponent of the directed Ising (DI) universality class for $σ\ge 1$. Investigating the behavior of the density in the absorbing phase, we argue that $σ=1$ is indeed the threshold that separates the DI and non-DI critical behavior. We also show by Monte Carlo simulations that branching bias with symmetric hopping exhibits the same critical behavior as the BAWL.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Crossover behaviors in branching annihilating attracting walk

We introduce branching annihilating attracting walk (BAAW) in one dimension. The attracting walk is implemented by a biased hopping in such a way that a particle prefers hopping to a nearest neighbor located on the side where the nearest particle is found within the range of attraction. We study the BAAW with four offspring by extensive Monte Carlo simulation. At first, we find the critical exponents of the BAAW with infinite range of attraction, which are different from those of the directed Ising (DI) universality class. Our results are consistent with the recent observation [B. Daga and P. Ray, Phys. Rev. E {\bf 99}, 032104 (2019)]. Then, by studying crossover behaviors, we show that as far as the range of attraction is finite the BAAW belongs to the DI class. We conclude that the origin of non-DI critical behavior of the BAAW with infinite range of attraction is the long-range nature of the attraction.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Order-parameter critical exponent of absorbing phase transitions in one-dimensional systems with two symmetric absorbing states

Via extensive Monte Carlo simulations along with systematic analyses of corrections to scaling, we estimate the order parameter critical exponent $β$ of absorbing phase transitions in systems with two symmetric absorbing states. The value of $β$ was conjectured to be $\frac{13}{14}\approx 0.93$ and Monte Carlo simulation studies in the literature have repeatedly reproduced values consistent with the conjecture. In this paper, we systematically estimate $β$ by analyzing the effective exponent after finding how strong corrections to scaling are. We show that the widely accepted numerical value of $β$ is not correct. Rather, we obtain $β= 1.020(5)$ from different models with two symmetric absorbing states.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Order-disorder transition in the two dimensional interacting monomer-dimer model: Ising criticality

We study the order-disorder transition of the two dimensional interacting monomer-dimer model (IMD) which has two symmetric absorbing states. To be self-contained, we first estimate numerically the dynamic exponent $z$ of the two dimensional Ising model. From the relaxation dynamics of the magnetization at the critical point, we obtain $β/(νz) = 0.057~650(12)$, or $z = 2.168 26(45)$, where $β= \frac{1}{8}$ and $ν=1$ are exactly known exponents. We, then, compare the critical relaxation of the order parameter at the transition point of the IMD with that of the Ising model. We found that the critical relaxation exponent $β/(νz)$ is in good agreement with the Ising model, unlike the recent claim by Nam et al [JSTAT {\bf (2014)}, P08011]. We also claim that the Binder cumulant is not an efficient quantity to locate the order-disorder transition point of the model with two symmetric absorbing states.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Block renormalization study on the nonequilibrium chiral Ising model

We present a numerical study on the ordering dynamics of a one-dimensional nonequilibrium Ising spin system with chirality. This system is characterized by a direction-dependent spin update rule. Pairs of $+-$ spins can flip to $++$ or $--$ with probability $(1-u)$ or to $-+$ with probability $u$ while $-+$ pairs are frozen. The system was found to evolve into the ferromagnetic ordered state at any $u<1$ exhibiting the power-law scaling of the characteristic length scale $ξ\sim t^{1/z}$ and the domain wall density $ρ\sim t^{-δ}$. The scaling exponents $z$ and $δ$ were found to vary continuously with the parameter $u$. In order to establish the anomalous power-law scaling firmly, we perform the block spin renormalization analysis proposed by Basu and Hinrichsen [U. Basu and H. Hinrichsen, J. Stat. Mech. (2011) P11023]. Domain walls of $b$ sites are coarse-grained into a block spin $σ^b$, and the relative frequencies of two-block patterns $σ^b_1 σ^b_2$ are measured in the $b\to\infty$ and $t\to\infty$ limit. These indices are expected to be universal. By performing extensive Monte Carlo simulations, we find that the indices also vary continuously with $u$ and that their values are consistent with the scaling exponents found in the previous study. This study serves as another evidence for the claim that the nonequilibrium chiral Ising model displays the power-law scaling behavior with continuously varying exponents.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Critical decay exponent of the pair contact process with diffusion

We investigate the one-dimensional pair contact process with diffusion (PCPD) by extensive Monte Carlo simulations, mainly focusing on the critical density decay exponent $δ$. To obtain an accurate estimate of $δ$, we first find the strength of corrections to scaling using the recently introduced method [S.-C. Park. J. Korean Phys. Soc. {\bf 62}, 469 (2013)]. For small diffusion rate ($d\le 0.5$), the leading corrections-to-scaling term is found to be $\sim t^{-0.15}$, whereas for large diffusion rate ($d=0.95$) it is found to be $\sim t^{-0.5}$. After finding the strength of corrections to scaling, effective exponents are systematically analyzed to conclude that the value of critical decay exponent $δ$ is $0.173(3)$ irrespective of $d$. This value should be compared with the critical decay exponent of the directed percolation, 0.1595. In addition, we study two types of crossover. At $d=0$, the phase boundary is discontinuous and the crossover from the pair contact process to the PCPD is found to be described by the crossover exponent $ϕ= 2.6(1)$. We claim that the discontinuity of the phase boundary cannot be consistent with the theoretical argument supporting the hypothesis that the PCPD should belong to the DP. At $d=1$, the crossover from the mean field PCPD to the PCPD is described by $ϕ= 2$ which is argued to be exact.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

High-precision Estimate of the Critical Exponents for the Directed Ising Universality Class

With extensive Monte Carlo simulations, we present high-precision estimates of the critical exponents of branching annihilating random walks with two offspring, a prototypical model of the directed Ising universality class in one dimension. To estimate the exponents accurately, we propose a systematic method to find corrections to scaling whose leading behavior is supposed to take the form $t^{-χ}$ in the long-time limit at the critical point. Our study shows that $χ\approx 0.75$ for the number of particles in defect simulations and $χ\approx 0.5$ for other measured quantities, which should be compared with the widely used value of $χ= 1$. Using $χ$ so obtained, we analyze the effective exponents to find that $β/ν_\| = 0.2872(2)$, $z = 1.7415(5)$, $η= 0.0000(2)$, and accordingly, $β/ν_\perp = 0.5000(6)$. Our numerical results for $β/ν_\|$ and $z$ are clearly different from the conjectured rational numbers $β/ν_\| = \frac{2}{7} \approx 0.2857$, $z = \frac{7}{4}= 1.75$ by Jensen [Phys. Rev. E, {\bf 50}, 3623 (1994)]. Our result for $β/ν_\perp$, however, is consistent with $\frac{1}{2}$, which is believed to be exact.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Coarsening Dynamics of Nonequilibrium Chiral Ising Models

We investigate a nonequilibrium coarsening dynamics of a one-dimensional Ising spin system with chirality. Only spins at domain boundaries are updated so that the model undergoes a coarsening to either of equivalent absorbing states with all spins + or -. Chirality is imposed by assigning different transition rates to events at down (+-) kinks and up (-+) kinks. The coarsening is characterized by power-law scalings of the kink density $ρ\sim t^{-δ}$ and the characteristic length scale $ξ\sim t^{1/z}$ with time $t$. Surprisingly the scaling exponents vary continuously with model parameters, which is not the case for systems without chirality. These results are obtained from extensive Monte Carlo simulations and spectral analyses of the time evolution operator. Our study uncovers the novel universality class of the coarsening dynamics with chirality.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Order-disorder transition in a model with two symmetric absorbing states

We study a model of two-dimensional interacting monomers which has two symmetric absorbing states and exhibits two kinds of phase transition; one is an order-disorder transition and the other is an absorbing phase transition. Our focus is around the order-disorder transition, and we investigate whether this transition is described by the critical exponents of the two-dimensional Ising model. By analyzing the relaxation dynamics of "staggered magnetization," the finite-size scaling, and the behavior of the magnetization in the presence of a symmetry-breaking field, we show that this model should belong to the Ising universality class. Our results along with the universality hypothesis support the idea that the order-disorder transition in two-dimensional models with two symmetric absorbing states is of the Ising universality class, contrary to the recent claim [K. Nam et al., J. Stat. Mech.: Theory Exp. (2011) L06001]. Furthermore, we illustrate that the Binder cumulant could be a misleading guide to the critical point in these systems.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Critical dynamics of nonconserved $N$-vector model with anisotropic nonequilibrium perturbations

We study dynamic field theories for nonconserving $N$-vector models that are subject to spatial-anisotropic bias perturbations. We first investigate the conditions under which these field theories can have a single length scale. When N=2 or $N \ge 4$, it turns out that there are no such field theories, and, hence, the corresponding models are pushed by the bias into the Ising class. We further construct nontrivial field theories for N=3 case with certain bias perturbations and analyze the renormalization-group flow equations. We find that the three-component systems can exhibit rich critical behavior belonging to two different universality classes.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Absence of the discontinuous transition in the one-dimensional triplet creation model

Although Hinrichsen in his unpublished work theoretically rebutted the possibility of the discontinuous transition in one-dimensional nonequilibrium systems unless there are additional conservation laws, long-range interactions, macroscopic currents, or special boundary conditions, we have recently observed the resurrection of the claim that the triplet creation model (TC) introduced by Dickman and Tomé [Phys. Rev. E {\bf 44}, 4833 (1991)] would show the discontinuous transition. By extensive simulations, however, we find that the one-dimensional TC does belong to the directed percolation universality class even for larger diffusion constant than the suggested tricritical point in the literature. Furthermore, we find that the phase boundary is well described by the crossover from the mean field to the directed percolation, which supports the claim that the one-dimensional TC does not exhibit a discontinuous transition.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Crossover from the parity-conserving pair contact process with diffusion to other universality classes

The pair contact process with diffusion (PCPD) with modulo 2 conservation (\pcpdt) [$2A\to 4A$, $2A\to 0$] is studied in one dimension, focused on the crossover to other well established universality classes: the directed Ising (DI) and the directed percolation (DP). First, we show that the \pcpdt shares the critical behaviors with the PCPD, both with and without directional bias. Second, the crossover from the \pcpdt to the DI is studied by including a parity-conserving single-particle process ($A \to 3A$). We find the crossover exponent $1/ϕ_1 = 0.57(3)$, which is argued to be identical to that of the PCPD-to-DP crossover by adding $A \to 2A$. This suggests that the PCPD universality class has a well defined fixed point distinct from the DP. Third, we study the crossover from a hybrid-type reaction-diffusion process belonging to the DP [$3A\to 5A$, $2A\to 0$] to the DI by adding $A \to 3A$. We find $1/ϕ_2 = 0.73(4)$ for the DP-to-DI crossover. The inequality of $ϕ_1$ and $ϕ_2$ further supports the non-DP nature of the PCPD scaling. Finally, we introduce a symmetry-breaking field in the dual spin language to study the crossover from the \pcpdt to the DP. We find $1/ϕ_3 = 1.23(10)$, which is associated with a new independent route from the PCPD to the DP.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Reentrant phase transition in a predator-prey model

We numerically investigate the six-species predator-prey game in complex networks as well as in $d$-dimensional hypercubic lattices with $d=1,2,..., 6$. The interaction topology of the six species contains two loops, each of which is composed of cyclically predating three species. As the mutation rate $P$ is lowered below the well-defined phase transition point, the $Z_2$ symmetry related with the interchange of the two loops is spontaneously broken, and it has been known that the system develops the defensive alliance in which three cyclically predating species defend each other against the invasion of other species. In the small-world network structure characterized by the rewiring probability $α$, the phase diagram shows the reentrant behavior as $α$ is varied, indicating a twofold role of the shortcuts. In $d$-dimensional regular hypercubic lattices, the system also exhibits the reentrant phase transition as $d$ is increased. We identify universality class of the phase transition and discuss the proper mean-field limit of the system.

0 citations0 reviewsNo code linkedOpen access
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}$.

0 citations0 reviewsNo code linkedOpen access