Source author record

Seung Ki Baek

Seung Ki Baek 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

46works
14topics
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

46 published item(s)

preprint2022arXiv

A second-order stability analysis for the continuous model of indirect reciprocity

Reputation is one of key mechanisms to maintain human cooperation, but its analysis gets complicated if we consider the possibility that reputation does not reach consensus because of erroneous assessment. The difficulty is alleviated if we assume that reputation and cooperation do not take binary values but have continuous spectra so that disagreement over reputation can be analysed in a perturbative way. In this work, we carry out the analysis by expanding the dynamics of reputation to the second order of perturbation under the assumption that everyone initially cooperates with good reputation. The second-order theory clarifies the difference between Image Scoring and Simple Standing in that punishment for defection against a well-reputed player should be regarded as good for maintaining cooperation. Moreover, comparison among the leading eight shows that the stabilizing effect of justified punishment weakens if cooperation between two ill-reputed players is regarded as bad. Our analysis thus explains how Simple Standing achieves a high level of stability by permitting justified punishment and also by disregarding irrelevant information in assessing cooperation. This observation suggests which factors affect the stability of a social norm when reputation can be perturbed by noise.

preprint2022arXiv

Democracy and polarization in the National Assembly of the Republic of Korea

The median-voter hypothesis predicts convergence of party platforms across a one-dimensional political spectrum during majoritarian elections. Assuming that the convergence is reflected in legislative activity, we study the time evolution of political polarization in the National Assembly of the Republic of Korea for the past 70 years. By projecting the correlation of lawmakers onto the first principal axis, we observe a high degree of polarization from the early 1960's to the late 1980's before democratization. As predicted by the hypothesis, it showed a sharp decrease when party politics were revived in 1987. Since then, the political landscape has become more and more multi-dimensional under the action of party politics, which invalidates the assumption behind the hypothesis. For comparison, we also analyze co-sponsorship in the United States House of Representatives from 1979 to 2020, whose correlation matrix has been constantly high-dimensional throughout the observation period. Our analysis suggests a pattern of polarization evolving with democratic development, from which we can argue the power and the limitation of the median-voter hypothesis as an explanation of real politics.

preprint2022arXiv

Social norms in indirect reciprocity with ternary reputations

Indirect reciprocity is a key mechanism that promotes cooperation in social dilemmas by means of reputation. Although it has been a common practice to represent reputations by binary values, either `good' or `bad', such a dichotomy is a crude approximation considering the complexity of reality. In this work, we studied norms with three different reputations, i.e., `good', `neutral', and `bad'. Through massive supercomputing for handling more than thirty billion possibilities, we fully identified which norms achieve cooperation and possess evolutionary stability against behavioural mutants. By systematically categorizing all these norms according to their behaviours, we found similarities and dissimilarities to their binary-reputation counterpart, the leading eight. We obtained four rules that should be satisfied by the successful norms, and the behaviour of the leading eight can be understood as a special case of these rules. A couple of norms that show counter-intuitive behaviours are also presented. We believe the findings are also useful for designing successful norms with more general reputation systems.

preprint2021arXiv

Friendly-rivalry solution to the iterated $n$-person public-goods game

Repeated interaction promotes cooperation among rational individuals under the shadow of future, but it is hard to maintain cooperation when a large number of error-prone individuals are involved. One way to construct a cooperative Nash equilibrium is to find a `friendly-rivalry' strategy, which aims at full cooperation but never allows the co-players to be better off. Recently it has been shown that for the iterated Prisoner's Dilemma in the presence of error, a friendly rival can be designed with the following five rules: Cooperate if everyone did, accept punishment for your own mistake, punish defection, recover cooperation if you find a chance, and defect in all the other circumstances. In this work, we construct such a friendly-rivalry strategy for the iterated $n$-person public-goods game by generalizing those five rules. The resulting strategy makes a decision with referring to the previous $m=2n-1$ rounds. A friendly-rivalry strategy for $n=2$ inherently has evolutionary robustness in the sense that no mutant strategy has higher fixation probability in this population than that of a neutral mutant. Our evolutionary simulation indeed shows excellent performance of the proposed strategy in a broad range of environmental conditions when $n= 2$ and $3$.

preprint2020arXiv

Co-sponsorship analysis of party politics in the 20th National Assembly of Republic of Korea

We investigate co-sponsorship among lawmakers by applying the principal-component analysis to the bills introduced in the 20th National Assembly of Korea. The most relevant factor for co-sponsorship is their party membership, and we clearly observe a signal of a third-party system in action. To identify other factors than the party influence, we analyze how lawmakers are clustered inside each party, and the result shows significant similarity between their committee membership and co-sponsorship in case of the ruling party. In addition, by monitoring each lawmaker's similarity to the average behavior of his or her party, we have found that it begins to decrease approximately one month before the lawmaker actually changes the party membership.

preprint2016arXiv

Combination with anti-tit-for-tat remedies problems of tit-for-tat

One of the most important questions in game theory concerns how mutual cooperation can be achieved and maintained in a social dilemma. In Axelrod's tournaments of the iterated prisoner's dilemma, Tit-for-Tat (TFT) demonstrated the role of reciprocity in the emergence of cooperation. However, the stability of TFT does not hold in the presence of implementation error, and a TFT population is prone to neutral drift to unconditional cooperation, which eventually invites defectors. We argue that a combination of TFT and anti-TFT (ATFT) overcomes these difficulties in a noisy environment, provided that ATFT is defined as choosing the opposite to the opponent's last move. According to this TFT-ATFT strategy, a player normally uses TFT; turns to ATFT upon recognizing his or her own error; returns to TFT either when mutual cooperation is recovered or when the opponent unilaterally defects twice in a row. The proposed strategy provides simple and deterministic behavioral rules for correcting implementation error in a way that cannot be exploited by the opponent, and suppresses the neutral drift to unconditional cooperation.

preprint2016arXiv

Comparing reactive and memory-one strategies of direct reciprocity

Direct reciprocity is a mechanism for the evolution of cooperation based on repeated interactions. When individuals meet repeatedly, they can use conditional strategies to enforce cooperative outcomes that would not be feasible in one-shot social dilemmas. Direct reciprocity requires that individuals keep track of their past interactions and find the right response. However, there are natural bounds on strategic complexity: Humans find it difficult to remember past interactions accurately, especially over long timespans. Given these limitations, it is natural to ask how complex strategies need to be for cooperation to evolve. Here, we study stochastic evolutionary game dynamics in finite populations to systematically compare the evolutionary performance of reactive strategies, which only respond to the co-player's previous move, and memory-one strategies, which take into account the own and the co-player's previous move. In both cases, we compare deterministic strategy and stochastic strategy spaces. For reactive strategies and small costs, we find that stochasticity benefits cooperation, because it allows for generous-tit-for-tat. For memory one strategies and small costs, we find that stochasticity does not increase the propensity for cooperation, because the deterministic rule of win-stay, lose-shift works best. For memory one strategies and large costs, however, stochasticity can augment cooperation.

preprint2016arXiv

Particle in a box with a time-dependent $δ$-function potential

In quantum information processing, one often considers inserting a barrier into a box containing a particle to generate one bit of Shannon entropy. We formulate this problem as a one-dimensional Schrödinger equation with a time-dependent $δ$-function potential. It is a natural generalization of the particle in a box, a canonical example of quantum mechanics, and we present analytic and numerical investigations on this problem. After deriving an exact Volterra-type integral equation, composed of an infinite sum of modes, we show that approximate formulas with the lowest-frequency modes correctly capture the qualitative behavior of the wave function. If we take into account hundreds of modes, our numerical calculation shows that the quantum adiabatic theorem actually gives a very good approximation even if the barrier height diverges within finite time, as long as it is sufficiently longer than the characteristic time scale of the particle. In particular, if the barrier is slowly inserted at an asymmetric position, the particle is localized by the insertion itself, in accordance with a prediction of the adiabatic theorem. On the other hand, when the barrier is inserted quickly, the wave function becomes rugged after the insertion because of the energy transfer to the particle. Regardless of the position of the barrier, the fast insertion leaves the particle unlocalized so that we can obtain meaningful information by a which-side measurement. Our numerical procedure provides a precise way to calculate the wave function throughout the process, from which one can estimate the amount of this information for an arbitrary insertion protocol.

preprint2015arXiv

Anomalous response in the vicinity of spontaneous symmetry breaking

We propose a mechanism to induce negative AC permittivity in the vicinity of a ferroelectric phase transition involved with spontaneous symmetry breaking. This mechanism makes use of responses at low frequency, yielding a high gain and a large phase delay, when the system jumps over the free-energy barrier with the aid of external fields. We illustrate the mechanism by analytically studying spin models with the Glauber-typed dynamics under periodic perturbations. Then, we show that the scenario is supported by numerical simulations of mean-field as well as two-dimensional spin systems.

preprint2015arXiv

Evolution of popularity in given names

An individual's identity in a human society is specified by his or her name. Differently from family names, usually inherited from fathers, a given name for a child is often chosen at the parents' disposal. However, their decision cannot be made in a vacuum but affected by social conventions and trends. Furthermore, such social pressure changes in time, as new names gain popularity while some other names are gradually forgotten. In this paper, we investigate how popularity of given names has evolved over the last century by using datasets collected in Korea, the province of Quebec in Canada, and the United States. In each of these countries, the average popularity of given names exhibits typical patterns of rise and fall with a time scale of about one generation. We also observe that notable changes of diversity in given names signal major social changes.

preprint2015arXiv

Interrupted coarsening in the zero-temperature kinetic Ising chain driven by a periodic external field

If quenched to zero temperature, the one-dimensional Ising spin chain undergoes coarsening, whereby the density of domain walls decays algebraically in time. We show that this coarsening process can be interrupted by exerting a rapidly oscillating periodic field with enough strength to compete with the spin-spin interaction. By analyzing correlation functions and the distribution of domain lengths both analytically and numerically, we observe nontrivial correlation with more than one length scale at the threshold field strength.

preprint2015arXiv

Nash equilibrium and evolutionary dynamics in semifinalists' dilemma

We consider a tournament among four equally strong semifinalists. The players have to decide how much stamina to use in the semifinals, provided that the rest is available in the final and the third-place playoff. We investigate optimal strategies for allocating stamina to the successive matches when players' prizes (payoffs) are given according to the tournament results. From the basic assumption that the probability to win a match follows a nondecreasing function of stamina difference, we present symmetric Nash equilibria for general payoff structures. We find three different phases of the Nash equilibria in the payoff space. First, when the champion wins a much bigger payoff than the others, any pure strategy can constitute a Nash equilibrium as long as all four players adopt it in common. Second, when the first two places are much more valuable than the other two, the only Nash equilibrium is such that everyone uses a pure strategy investing all stamina in the semifinal. Third, when the payoff for last place is much smaller than the others, a Nash equilibrium is formed when every player adopts a mixed strategy of using all or none of its stamina in the semifinals. In a limiting case that only last place pays the penalty, this mixed-strategy profile can be proved to be a unique symmetric Nash equilibrium, at least when the winning probability follows a Heaviside step function. Moreover, by using this Heaviside step function, we study the tournament by using evolutionary replicator dynamics to obtain analytic solutions, which reproduces the corresponding Nash equilibria on the population level and gives information on dynamic aspects.

preprint2015arXiv

Quantum Isothermal Reversible Process of Particles in a Box with a Delta Potential

For an understanding of a heat engine working in the microscopic scale, it is often necessary to estimate the amount of reversible work extracted by isothermal expansion of the quantum gas used as its working substance. We consider an engine with a movable wall, modeled as an infinite square well with a delta peak inside. By solving the resulting one-dimensional Schrödinger equation, we obtain the energy levels and the thermodynamic potentials. Our result shows how quantum tunneling degrades the engine by decreasing the amount of reversible work during the isothermal expansion.

preprint2015arXiv

Symmetry restoration by pricing in a duopoly of perishable goods

Competition is a main tenet of economics, and the reason is that a perfectly competitive equilibrium is Pareto-efficient in the absence of externalities and public goods. Whether a product is selected in a market crucially relates to its competitiveness, but the selection in turn affects the landscape of competition. Such a feedback mechanism has been illustrated in a duopoly model by Lambert et al., in which a buyer's satisfaction is updated depending on the {\em freshness} of a purchased product. The probability for buyer $n$ to select seller $i$ is assumed to be $p_{n,i} \propto e^{ S_{n,i}/T}$, where $S_{n,i}$ is the buyer's satisfaction and $T$ is an effective temperature to introduce stochasticity. If $T$ decreases below a critical point $T_c$, the system undergoes a transition from a symmetric phase to an asymmetric one, in which only one of the two sellers is selected. In this work, we extend the model by incorporating a simple price system. By considering a greed factor $g$ to control how the satisfaction depends on the price, we argue the existence of an oscillatory phase in addition to the symmetric and asymmetric ones in the $(T,g)$ plane, and estimate the phase boundaries through mean-field approximations. The analytic results show that the market preserves the inherent symmetry between the sellers for lower $T$ in the presence of the price system, which is confirmed by our numerical simulations.

preprint2015arXiv

Theory of fads: Traveling-wave solution of evolutionary dynamics in a one-dimensional trait space

We consider an infinite-sized population where an infinite number of traits compete simultaneously. The replicator equation with a diffusive term describes time evolution of the probability distribution over the traits due to selection and mutation on a mean-field level. We argue that this dynamics can be expressed as a variant of the Fisher equation with high-order correction terms. The equation has a traveling-wave solution, and the phase-space method shows how the wave shape depends on the correction. We compare this solution with empirical time-series data of given names in Quebec, treating it as a descriptive model for the observed patterns. Our model explains the reason that many names exhibit a similar pattern of the rise and fall as time goes by. At the same time, we have found that their dissimilarities are also statistically significant.

preprint2014arXiv

Internal energy density of the critical three-state Potts model on the kagome lattice

It has been conjectured that the internal energy density of the Potts model on a semi-infinite strip with a width $L$ does not have any finite-size corrections at the critical point $K=K_c$. By factorizing the transfer matrix for the kagome lattice with larger widths, we have found that this conjecture is not correct in that the internal energy density slightly varies with $L$ at the critical point. From this size dependence of the internal energy density, we obtain an upper bound as $K_c < 1.0565615$, which is close to a recent estimate $K_c^{\rm JS} = 1.0565600(7)$ by Jacobsen and Scullard [arXiv:1204.0622]. We also obtain a lower bound as $K_c > 1.0560$ by calculating the correlation length along the strips.

preprint2014arXiv

Nonequilibrium steady state of the kinetic Glauber-Ising model under an alternating magnetic field

When periodically driven by an external magnetic field, a spin system can enter a phase of steady entrained oscillations with nonequilibrium probability distribution function. We consider an arbitrary magnetic field switching its direction with frequency comparable with the spin-flip rate and show that the resulting nonequilibrium probability distribution can be related to the system equilibrium distribution in the presence of a constant magnetic field of the same magnitude. We derive convenient approximate expressions for this exact relation and discuss their implications.

preprint2014arXiv

Stochastic resonance in the two-dimensional q-state clock models

We numerically study stochastic resonance in the two-dimensional q-state clock models from q = 2 to 7 under a weak oscillating magnetic field. As in the mean-field case, we observe double resonance peaks, but the detailed response strongly depends on the direction of the field modulation for q >= 5 where the quasiliquid phase emerges. We explain this behavior in terms of free-energy landscapes on the two-dimensional magnetization plane.

preprint2014arXiv

Surname statistics - Crossing the boundary between disciplines

Regarding "Surname distribution in population genetics and in statistical physics" by Rossi, we add comments on the following two points: One is a difference in his renormalization-group formulation from the traditional one in statistical physics despite apparent similarity. The other is the potential importance of the entropy maximization principle under constraints governing the evolution of surname statistics.

preprint2014arXiv

Universal statistics of the knockout tournament

We study statistics of the knockout tournament, where only the winner of a fixture progresses to the next. We assign a real number called competitiveness to each contestant and find that the resulting distribution of prize money follows a power law with an exponent close to unity if the competitiveness is a stable quantity and a decisive factor to win a match. Otherwise, the distribution is found narrow. The existing observation of power law distributions in various kinds of real sports tournaments therefore suggests that the rules of those games are constructed in such a way that it is possible to understand the games in terms of the contestants' inherent characteristics of competitiveness.

preprint2014arXiv

Zero-one-only process: a correlated random walk with a stochastic ratchet

The investigation of random walks is central to a variety of stochastic processes in physics, chemistry, and biology. To describe a transport phenomenon, we study a variant of the one-dimensional persistent random walk, which we call a zero-one-only process. It makes a step in the same direction as the previous step with probability $p$, and stops to change the direction with $1-p$. By using the generating-function method, we calculate its characteristic quantities such as the statistical moments and probability of the first return.

preprint2013arXiv

Faster is More Different: Mean-Field Dynamics of Innovation Diffusion

Based on a recent model of paradigm shifts by Bornholdt et al., we studied mean-field opinion dynamics in an infinite population where an infinite number of ideas compete simultaneously with their values publicly known. We found that a highly innovative society is not characterized by heavy concentration in highly valued ideas: Rather, ideas are more broadly distributed in a more innovative society with faster progress, provided that the rate of adoption is constant, which suggests a positive correlation between innovation and technological disparity. Furthermore, the distribution is generally skewed in such a way that the fraction of innovators is substantially smaller than has been believed in conventional innovation-diffusion theory based on normality. Thus, the typical adoption pattern is predicted to be asymmetric with slow saturation in the ideal situation, which is compared with empirical data sets.

preprint2013arXiv

Phase transition in a coevolving network of conformist and contrarian voters

In the coevolving voter model, each voter has one of two diametrically opposite opinions, and a voter encountering a neighbor with the opposite opinion may either adopt it or rewire the connection to another randomly chosen voter sharing the same opinion. As we smoothly change the relative frequency of rewiring compared to that of adoption, there occurs a phase transition between an active phase and a frozen phase. By performing extensive Monte Carlo calculations, we show that the phase transition is characterized by critical exponents β=0.54(1) and ν =1.5(1), which differ from the existing mean-field-type prediction. We furthermore extend the model by introducing a contrarian type that tries to have neighbors with the opposite opinion, and show that the critical behavior still belongs to the same universality class irrespective of such contrarians' fraction.

preprint2013arXiv

Residual discrete symmetry of the five-state clock model

It is well-known that the $q$-state clock model can exhibit a Kosterlitz-Thouless (KT) transition if $q$ is equal to or greater than a certain threshold, which has been believed to be five. However, recent numerical studies indicate that helicity modulus does not vanish in the high-temperature phase of the five-state clock model as predicted by the KT scenario. By performing Monte Carlo calculations under the fluctuating twist boundary condition, we show that it is because the five-state clock model does not have the fully continuous U(1) symmetry even in the high-temperature phase while the six-state clock model does. We suggest that the upper transition of the five-state clock model is actually a weaker cousin of the KT transition so that it is $q \ge 6$ that exhibits the genuine KT behavior.

preprint2012arXiv

Bounds of percolation thresholds on hyperbolic lattices

We analytically study bond percolation on hyperbolic lattices obtained by tiling a hyperbolic plane with constant negative Gaussian curvature. The quantity of our main concern is $p_{c2}$, the value of occupation probability where a unique unbounded cluster begins to emerge. By applying the substitution method to known bounds of the order-5 pentagonal tiling, we show that $p_{c2} \ge 0.382 508$ for the order-5 square tiling, $p_{c2} \ge 0.472 043$ for its dual, and $p_{c2} \ge 0.275 768$ for the order-5-4 rhombille tiling.

preprint2012arXiv

Cluster-size heterogeneity in the two-dimensional Ising model

We numerically investigate the heterogeneity in cluster sizes in the two-dimensional Ising model and verify its scaling form recently proposed in the context of percolation problems [Phys. Rev. E 84, 010101(R) (2011)]. The scaling exponents obtained via the finite-size scaling analysis are shown to be consistent with theoretical values of the fractal dimension $d_f$ and the Fisher exponent $τ$ for the cluster distribution. We also point out that strong finite-size effects exist due to the geometric nature of the cluster-size heterogeneity.

preprint2012arXiv

Critical condition of the water-retention model

We study how much water can be retained without leaking through boundaries when each unit square of a two-dimensional lattice is randomly assigned a block of unit bottom area but with different heights from zero to $n-1$. As more blocks are put into the system, there exists a phase transition beyond which the system retains a macroscopic volume of water. We locate the critical points and verify that the criticality belongs to the two-dimensional percolation universality class. If the height distribution can be approximated as continuous for large $n$, the system is always close to a critical point and the fraction of the area below the resulting water level is given by the percolation threshold. This provides a universal upper bound of areas that can be covered by water in a random landscape.

preprint2012arXiv

Double stochastic resonance in the mean-field $q$-state clock models

A magnetic system with a phase transition at temperature $T_c$ may exhibit double resonance peaks under a periodic external magnetic field because the time scale matches the external frequency at two different temperatures, one above $T_c$ and the other below $T_c$. We study the double resonance phenomena for the mean-field $q$-state clock model based on the heat-bath-typed master equation. We find double peaks as observed in the kinetic Ising case ($q=2$) for all $q\ge 4$, but for the three-state clock model ($q=3$), the existence of double peaks is possible only above a certain external frequency since it undergoes a discontinuous phase transition.

preprint2012arXiv

Dworkin's Paradox

How to distribute welfare in a society is a key issue in the subject of distributional justice, which is deeply involved with notions of fairness. Following a thought experiment by Dworkin, this work considers a society of individuals with different preferences on the welfare distribution and an official to mediate the coordination among them. Based on a simple assumption that an individual's welfare is proportional to how her preference is fulfilled by the actual distribution, we show that an egalitarian preference is a strict Nash equilibrium and can be favorable even in certain inhomogeneous situations. These suggest how communication can encourage and secure a notion of fairness.

preprint2012arXiv

Upper transition point for percolation on the enhanced binary tree: A sharpened lower bound

Hyperbolic structures are obtained by tiling a hyperbolic surface with negative Gaussian curvature. These structures generally exhibit two percolation transitions: a system-wide connection can be established at a certain occupation probability $p=p_{c1}$ and there emerges a unique giant cluster at $p_{c2} > p_{c1}$. There have been debates about locating the upper transition point of a prototypical hyperbolic structure called the enhanced binary tree (EBT), which is constructed by adding loops to a binary tree. This work presents its lower bound as $p_{c2} \gtrsim 0.55$ by using phenomenological renormalization-group methods and discusses some solvable models related to the EBT.

preprint2012arXiv

Vortex interaction on curved surfaces

The vortex-excitation energy on a sphere can be obtained by using the stereographic projection. By applying this method, we calculate the energy needed to create a vortex on a surface with a constant negative curvature. It is found that the energy is a linear function of the radius of the vortex. In accordance with this result, the interaction energy between a pair of vortices is also found to change linearly with the vortex separation distance. Explicit vortex configurations are obtained numerically with this interaction.

preprint2011arXiv

A Paradoxical Property of the Monkey Book

A "monkey book" is a book consisting of a random distribution of letters and blanks, where a group of letters surrounded by two blanks is defined as a word. We compare the statistics of the word distribution for a monkey book with the corresponding distribution for the general class of random books, where the latter are books for which the words are randomly distributed. It is shown that the word distribution statistics for the monkey book is different and quite distinct from a typical sampled book or real book. In particular the monkey book obeys Heaps' power law to an extraordinary good approximation, in contrast to the word distributions for sampled and real books, which deviate from Heaps' law in a characteristics way. The somewhat counter-intuitive conclusion is that a "monkey book" obeys Heaps' power law precisely because its word-frequency distribution is not a smooth power law, contrary to the expectation based on simple mathematical arguments that if one is a power law, so is the other.

preprint2011arXiv

Cluster Monte Carlo study of magnetic dipoles

We implement a cluster-update Monte Carlo algorithm to simulate magnetic dipoles of the XY-spin type confined in a two-dimensional plane. The long-range character and anisotropy in the dipole interaction are handled by using the Luijten-Blöte algorithm and the Dotsenko-Selke-Talapov algorithm, respectively. We have checked the performance of this cluster-update algorithm in comparison to the Metropolis algorithm and found that it equilibrated the system faster in terms of the number of flipped spins, although the overall computational complexity of the problem remained the same.

preprint2011arXiv

Critical temperatures of the three- and four-state Potts models on the kagome lattice

The value of the internal energy per spin is independent of the strip width for a certain class of spin systems on two dimensional infinite strips. It is verified that the Ising model on the kagome lattice belongs to this class through an exact transfer-matrix calculation of the internal energy for the two smallest widths. More generally, one can suggest an upper bound for the critical coupling strength $K_c(q)$ for the $q$-state Potts model from exact calculations of the internal energy for the two smallest strip widths. Combining this with the corresponding calculation for the dual lattice and using an exact duality relation enables us to conjecture the critical coupling strengths for the three- and four-state Potts models on the kagome lattice. The values are $K_c(q=3)=1.056~509~426~929~0$ and $K_c(q=4) = 1.149~360~587~229~2$, and the values can, in principle, be obtained to an arbitrary precision. We discuss the fact that these values are in the middle of earlier approximate results and furthermore differ from earlier conjectures for the exact values.

preprint2011arXiv

Hierarchical renormalization-group study on the planar bond-percolation problem

For certain hierarchical structures, one can study the percolation problem using the renormalization-group method in a very precise way. We show that the idea can be also applied to two-dimensional planar lattices by regarding them as hierarchical structures. Either a lower bound or an exact critical probability can be obtained with this method and the correlation-length critical exponent is approximately estimated as $ν\approx 1$.

preprint2011arXiv

Ising model on a hyperbolic plane with a boundary

A hyperbolic plane can be modeled by a structure called the enhanced binary tree. We study the ferromagnetic Ising model on top of the enhanced binary tree using the renormalization-group analysis in combination with transfer-matrix calculations. We find a reasonable agreement with Monte Carlo calculations on the transition point, and the resulting critical exponents suggest the mean-field surface critical behavior.

preprint2011arXiv

Kosterlitz-Thouless transition of magnetic dipoles on the two-dimensional plane

The universality class of a phase transition is often determined by factors like dimensionality and inherent symmetry. We study the magnetic dipole system in which the ground-state symmetry and the underlying lattice structure are coupled to each other in an intricate way. A two-dimensional (2D) square-lattice system of magnetic dipoles undergoes an order-disorder phase transition belonging to the 2D Ising universality class. According to Prakash and Henley [Phys. Rev. B {\bf 42}, 6572 (1990)], this can be related to the fourfold-symmetric ground states which suggests a similarity to the four-state clock model. Provided that this type of symmetry connection holds true, the magnetic dipoles on a honeycomb lattice, which possess sixfold-symmetric ground states, should exhibit a Kosterlitz-Thouless transition in accordance with the six-state clock model. This is verified through numerical simulations in the present investigation. However, it is pointed out that this symmetry argument does not always apply, which suggests that factors other than symmetry can be decisive for the universality class of the magnetic dipole system.

preprint2011arXiv

Quantum Monte Carlo study of the transverse-field quantum Ising model on infinite-dimensional structures

In a number of classical statistical-physical models, there exists a characteristic dimensionality called the upper critical dimension above which one observes the mean-field critical behavior. Instead of constructing high-dimensional lattices, however, one can also consider infinite-dimensional structures, and the question is whether this mean-field character extends to quantum-mechanical cases as well. We therefore investigate the transverse-field quantum Ising model on the globally coupled network and the Watts-Strogatz small-world network by means of quantum Monte Carlo simulations and the finite-size scaling analysis. We confirm that both the structures exhibit critical behavior consistent with the mean-field description. In particular, we show that the existing cumulant method has a difficulty in estimating the correct dynamic critical exponent and suggest that an order parameter based on the quantum-mechanical expectation value can be a practically useful numerical observable to determine critical behavior when there is no well-defined dimensionality.

preprint2011arXiv

The Ten Thousand Kims

In the Korean culture the family members are recorded in special family books. This makes it possible to follow the distribution of Korean family names far back in history. It is here shown that these name distributions are well described by a simple null model, the random group formation (RGF) model. This model makes it possible to predict how the name distributions change and these predictions are shown to be borne out. In particular, the RGF model predicts that, for married women entering a collection of family books in a certain year, the occurrence of the most common family name "Kim" should be directly proportional the total number of married women with the same proportionality constant for all the years. This prediction is also borne out to high degree. We speculate that it reflects some inherent social stability in the Korean culture. In addition, we obtain an estimate of the total population of the Korean culture down to year 500 AD, based on the RGF model and find about ten thousand Kims.

preprint2011arXiv

Zipf's law unzipped

Why does Zipf's law give a good description of data from seemingly completely unrelated phenomena? Here it is argued that the reason is that they can all be described as outcomes of a ubiquitous random group division: the elements can be citizens of a country and the groups family names, or the elements can be all the words making up a novel and the groups the unique words, or the elements could be inhabitants and the groups the cities in a country, and so on. A Random Group Formation (RGF) is presented from which a Bayesian estimate is obtained based on minimal information: it provides the best prediction for the number of groups with $k$ elements, given the total number of elements, groups, and the number of elements in the largest group. For each specification of these three values, the RGF predicts a unique group distribution $N(k)\propto \exp(-bk)/k^γ$, where the power-law index $γ$ is a unique function of the same three values. The universality of the result is made possible by the fact that no system specific assumptions are made about the mechanism responsible for the group division. The direct relation between $γ$ and the total number of elements, groups, and the number of elements in the largest group, is calculated. The predictive power of the RGF model is demonstrated by direct comparison with data from a variety of systems. It is shown that $γ$ usually takes values in the interval $1\leqγ\leq 2$ and that the value for a given phenomena depends in a systematic way on the total size of the data set. The results are put in the context of earlier discussions on Zipf's and Gibrat's laws, $N(k)\propto k^{-2}$ and the connection between growth models and RGF is elucidated.

preprint2010arXiv

Analytic results for the percolation transitions of the enhanced binary tree

Percolation for a planar lattice has a single percolation threshold, whereas percolation for a negatively curved lattice displays two separate thresholds. The enhanced binary tree (EBT) can be viewed as a prototype model displaying two separate percolation thresholds. We present an analytic result for the EBT model which gives two critical percolation threshold probabilities, $p_{c1}=1/2\sqrt{13}-3/2$ and $p_{c2}=1/2$, and yields a size-scaling exponent $Φ=\ln [\frac{p(1+p)}{1-p(1-p)}]/\ln 2$. It is inferred that the two threshold values give exact upper limits and that $p_{c1}$ is furthermore exact. In addition, we argue that $p_{c2}$ is also exact. The physics of the model and the results are described within the midpoint-percolation concept: Monte Carlo simulations are presented for the number of boundary points which are reached from the midpoint, and the results are compared to the number of routes from the midpoint to the boundary given by the analytic solution. These comparisons provide a more precise physical picture of what happens at the transitions. Finally, the results are compared to related works, in particular, the Cayley tree and Monte Carlo results for hyperbolic lattices as well as earlier results for the EBT model. It disproves a conjecture that the EBT has an exact relation to the thresholds of its dual lattice.

preprint2010arXiv

Comment on `Six-state clock model on the square lattice: Fisher zero approach with Wang-Landau sampling'

Hwang in [Phys. Rev. E 80, 042103 (2009)] suggested that the two transitions of the six-state clock model on the square lattice are not of the Kosterlitz-Thouless type. Here we show from simulations thatat the upper transition, the helicity modulus does make a discontinuous jump to zero. This gives strong evidence for a Kosterlitz-Thouless transition.

preprint2010arXiv

Equilibrium solution to the lowest unique positive integer game

We address the equilibrium concept of a reverse auction game so that no one can enhance the individual payoff by a unilateral change when all the others follow a certain strategy. In this approach the combinatorial possibilities to consider become very much involved even for a small number of players, which has hindered a precise analysis in previous works. We here present a systematic way to reach the solution for a general number of players, and show that this game is an example of conflict between the group and the individual interests.

preprint2010arXiv

Non-Kosterlitz-Thouless transitions for the $q$-state clock models

The $q$-state clock model with the cosine potential has a single phase transition for $q\leq4$ and two transitions for $q\geq5$. It is shown by Monte Carlo simulations that the helicity modulus for the five-state clock model ($q=5$) does not vanish at the high-temperature transition. This is in contrast to the clock models with $q\geq6$ for which the helicity modulus vanishes. This means that the transition for the five-state clock model differs from the Kosterlitz-Thouless (KT) transition. It is also shown that this change in the transition is caused by an interplay between the number of angular directions and the interaction potential: by slightly modifying the interaction potential, the KT transition for $q=6$ turns into the same non-KT transition. Likewise, the KT transition is recovered for $q=5$ when the Villain potential is used. Comparisons with other clock-model results are made and discussed.

preprint2010arXiv

Surface and bulk criticality in midpoint percolation

The concept of midpoint percolation has recently been applied to characterize the double percolation transitions in negatively curved structures. Regular $d$-dimensional hypercubic lattices are in the present work investigated using the same concept. Specifically, the site-percolation transitions at the critical thresholds are investigated for dimensions up to $d=10$ by means of the Leath algorithm. It is shown that the explicit inclusion of the boundaries provides a straightforward way to obtain critical indices, both for the bulk and surface parts. At and above the critical dimension $d=6$, it is found that the percolation cluster contains only a finite number of surface points in the infinite-size limit. This is in accordance with the expectation from studies of lattices with negative curvature. It is also found that the number of surface points, reached by the percolation cluster in the infinite limit, approaches 2d for large dimensions $d$. We also note that the size dependence in proliferation of percolating clusters for $d\ge 7$ can be obtained by solely counting surface points of the midpoint cluster.

preprint2009arXiv

Comment on `Monte Carlo simulation study of the two-stage percolation transition in enhanced binary trees'

The enhanced binary tree (EBT) is a nontransitive graph which has two percolation thresholds $p_{c1}$ and $p_{c2}$ with $p_{c1}<p_{c2}$. Our Monte Carlo study implies that the second threshold $p_{c2}$ is significantly lower than a recent claim by Nogawa and Hasegawa (J. Phys. A: Math. Theor. {\bf 42} (2009) 145001). This means that $p_{c2}$ for the EBT does not obey the duality relation for the thresholds of dual graphs $p_{c2}+\overline{p}_{c1}=1$ which is a property of a transitive, nonamenable, planar graph with one end. As in regular hyperbolic lattices, this relation instead becomes an inequality $p_{c2}+\overline{p}_{c1}<1$. We also find that the critical behavior is well described by the scaling form previously found for regular hyperbolic lattices.