Catalog footprint

What is connected

56works
25topics
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

56 published item(s)

preprint2022arXiv

Multiresolution community analysis of international trade networks

The international trade network is a complex system where multiple trade blocs with varying sizes coexist and overlap with each other. However, the resulting structures of community detection in trade networks are often inconsistent and fails to capture the complex landscape of international trade. To address these problems, we propose a multiresolution framework that aggregates all the configuration information from a range of resolutions. This allows us to consider trade communities of different sizes and illuminate the underlying hierarchical structure of trade networks and its constituting blocks. Furthermore, by measuring membership inconsistency (MeI) of each country and conducting multiple regression analysis with various economic and political indicators, we demonstrate that there exists a positive correlation between the external instability of countries and their structural inconsistency in terms of network topology.

preprint2021arXiv

Different Environmental Conditions in Genetic Algorithm

We propose an extended genetic algorithm (GA) with different local environmental conditions. Genetic entities, or configurations, are put on nodes in a ring structure, and location-dependent environmental conditions are applied for each entity. Our GA is motivated by the geographic aspect of natural evolution: Geographic isolation reduces the diversity in a local group, but at the same time, can enhance intergroup diversity. Mating of genetic entities across different environments can make it possible to search for broad area of the fitness landscape. We validate our extended GA for finding the ground state of three-dimensional spin-glass system and find that the use of different environmental conditions makes it possible to find the lower-energy spin configurations at relatively shorter computation time. Our extension of GA belongs to a meta-optimization method and thus can be applied for a broad research area in which finding of the optimal state in a shorter computation time is the key problem.

preprint2021arXiv

Power-grid stability predictions using transferable machine learning

Complex network analyses have provided clues to improve power-grid stability with the help of numerical models. The high computational cost of numerical simulations, however, has inhibited the approach, especially when it deals with the dynamic properties of power grids such as frequency synchronization. In this study, we investigate machine learning techniques to estimate the stability of power-grid synchronization. We test three different machine learning algorithms -- random forest, support vector machine, and artificial neural network -- training them with two different types of synthetic power grids consisting of homogeneous and heterogeneous input-power distribution, respectively. We find that the three machine learning models better predict the synchronization stability of power-grid nodes when they are trained with the heterogeneous input-power distribution than the homogeneous one. With the real-world power grids of Great Britain, Spain, France, and Germany, we also demonstrate that the machine learning algorithms trained on synthetic power grids are transferable to the stability prediction of the real-world power grids, which implies the prospective applicability of machine learning techniques on power-grid studies.

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.

preprint2020arXiv

Extracting hierarchical backbones from bipartite networks

We propose a method for extracting hierarchical backbones from a bipartite network. Our method leverages the observation that a hierarchical relationship between two nodes in a bipartite network is often manifested as an asymmetry in the conditional probability of observing the connections to them from the other node set. Our method estimates both the importance and direction of the hierarchical relationship between a pair of nodes, thereby providing a flexible way to identify the essential part of the networks. Using semi-synthetic benchmarks, we show that our method outperforms existing methods at identifying planted hierarchy while offering more flexibility. Application of our method to empirical datasets---a bipartite network of skills and individuals as well as the network between gene products and Gene Ontology (GO) terms---demonstrates the possibility of automatically extracting or augmenting ontology from data.

preprint2019arXiv

Scaling relations and finite-size scaling in gravitationally correlated lattice percolation models

In some systems, the connecting probability (and thus the percolation process) between two sites depends on the geometric distance between them. To understand such process, we propose gravitationally correlated percolation models for link-adding networks on the two-dimensional lattice $G$ with two strategies $S_{\rm max}$ and $S_{\rm min}$, to add a link $l_{i,j}$ to connect site $i$ and site $j$ with mass $m_i$ and $m_j$, respectively; $m_i$ and $m_j$ are sizes of the clusters which contain site $i$ and site $j$, respectively. The probability to add the link $l_{i,j}$ is related to the generalized gravity $g_{ij} \equiv m_i m_j/r_{ij}^d$, where $r_{ij}$ is the geometric distance between $i$ and $j$, and $d$ is an adjustable decaying exponent. In the beginning of the simulation, all sites of $G$ are occupied and there is no link. In the simulation process, two inter-cluster links $l_{i,j}$ and $l_{k,n}$ are randomly chosen and the generalized gravities $g_{ij}$ and $g_{kn}$ are computed. In the strategy $S_{\rm max}$, the link with larger generalized gravity is added. In the strategy $S_{\rm min}$, the link with smaller generalized gravity is added, which include percolation on the Erd\H os-Rényi random graph and the Achlioptas process of explosive percolation as the limiting cases, $d \to \infty$ and $d \to 0$, respectively. Adjustable strategies facilitate or inhibit the network percolation in a generic view. We calculate percolation thresholds $T_c$ and critical exponents $β$ by numerical simulations. We also obtain various finite-size scaling functions for the node fractions in percolating clusters or arrival of saturation length with different intervening strategies.

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

Finite-Time-Finite-Size Scaling of the Kuramoto Oscillators

Phase transition in its strict sense can only be observed in an infinite system, for which equilibration takes an infinitely long time at criticality. In numerical simulations, we are often limited both by the finiteness of the system size and by the finiteness of the observation time scale. We propose that one can overcome this barrier by measuring the nonequilibrium temporal relaxation for finite systems and by applying the finite-time-finite-size scaling (FTFSS) which systematically uses two scaling variables, one temporal and the other spatial. The FTFSS method yields a smooth scaling surface, and the conventional finite-size scaling curves can be viewed as proper cross sections of the surface. The validity of our FTFSS method is tested for the synchronization transition of Kuramoto models in the globally-coupled structure and in the small-world network structure. Our FTFSS method is also applied to the Monte-Carlo dynamics of the globally-coupled q-state clock model.

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.

preprint2015arXiv

Winding number excitation detects phase transition in one-dimensional XY model with variable interaction range

We numerically study the critical behavior of the one-dimensional XY model of the size N with variable interaction range L. As expected, the standard local order parameter of the magnetization is shown to well detect the mean-field type transition which occurs at any nonzero value of L/N. The system is particularly interesting since the underlying one-dimensional structure allows us to study the topological excitation of the winding number across the whole system even though the system shares the mean-field transition with the globally-coupled system. We propose a novel nonlocal order parameter based on the width of the winding number distribution which exhibits a clear signature of the transition nature of the system.

preprint2014arXiv

International Transmission of Shocks and Fragility of a Bank Network

The weighted and directed network of countries based on the number of overseas banks is analyzed in terms of its fragility to the banking crisis of one country. We use two different models to describe transmission of shocks, one local and the other global. Depending on the original source of the crisis, the overall size of crisis impacts is found to differ country by country. For the two-step local spreading model, it is revealed that the scale of the first impact is determined by the out-strength, the total number of overseas branches of the country at the origin of the crisis, while the second impact becomes more serious if the in-strength at the origin is increased. For the global spreading model, some countries named "triggers" are found to play important roles in shock transmission, and the importance of the feed-forward-loop mechanism is pointed out. We also discuss practical policy implications of the present work.

preprint2014arXiv

Matchmaker, Matchmaker, Make Me a Match: Migration of Populations via Marriages in the Past

The study of human mobility is both of fundamental importance and of great potential value. For example, it can be leveraged to facilitate efficient city planning and improve prevention strategies when faced with epidemics. The newfound wealth of rich sources of data---including banknote flows, mobile phone records, and transportation data---has led to an explosion of attempts to characterize modern human mobility. Unfortunately, the dearth of comparable historical data makes it much more difficult to study human mobility patterns from the past. In this paper, we present an analysis of long-term human migration, which is important for processes such as urbanization and the spread of ideas. We demonstrate that the data record from Korean family books (called "jokbo") can be used to estimate migration patterns via marriages from the past 750 years. We apply two generative models of long-term human mobility to quantify the relevance of geographical information to human marriage records in the data, and we find that the wide variety in the geographical distributions of the clans poses interesting challenges for the direct application of these models. Using the different geographical distributions of clans, we quantify the "ergodicity" of clans in terms of how widely and uniformly they have spread across Korea, and we compare these results to those obtained using surname data from the Czech Republic. To examine population flow in more detail, we also construct and examine a population-flow network between regions. Based on the correlation between ergodicity and migration in Korea, we identify two different types of migration patterns: diffusive and convective. We expect the analysis of diffusive versus convective effects in population flows to be widely applicable to the study of mobility and migration patterns across different cultures.

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

Structural Properties of Networks Grown via an Achlioptas Process

After the Achlioptas process (AP), which yields the so-called explosive percolation, was introduced, the number of papers on percolation phenomena has been literally exploding. Most of the existing studies, however, have focused only on the nature of phase transitions, not paying proper attention to the structural properties of the resulting networks, which compose the main theme of the present paper. We compare the resulting network structure of the AP with random networks and find, through observations of the distributions of the shortest-path length and the betweenness centrality in the giant cluster, that the AP makes the network less clustered and more fragile. Such structural characteristics are more directly seen by using snapshots of the network structures and are explained by the fact that the AP suppresses the formation of large clusters more strongly than the random process does. These structural differences between the two processes are shown to be less noticeable in growing networks than in static ones.

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

Allometric Exponent and Randomness

An allometric height-mass exponent $γ$ gives an approximative power-law relation $< M> \propto H^γ$ between the average mass $< M>$ and the height $H$, for a sample of individuals. The individuals in the present study are humans but could be any biological organism. The sampling can be for a specific age of the individuals or for an age-interval. The body-mass index (BMI) is often used for practical purposes when characterizing humans and it is based on the allometric exponent $γ=2$. It is here shown that the actual value of $γ$ is to large extent determined by the degree of correlation between mass and height within the sample studied: no correlation between mass and height means $γ=0$, whereas if there was a precise relation between mass and height such that all individuals had the same shape and density then $γ=3$. The connection is demonstrated by showing that the value of $γ$ can be obtained directly from three numbers characterizing the spreads of the relevant random Gaussian statistical distributions: the spread of the height and mass distributions together with the spread of the mass distribution for the average height. Possible implications for allometric relations in general are discussed.

preprint2013arXiv

Analysis of one-dimensional Yut-Nori game: winning strategy and avalanche-size distribution

In the Korean traditional board game Yut-Nori, teams compete by moving their pieces on the two-dimensional game board, and the team whose all pieces complete a round trip on the board wins. In every round, teams throw four wooden sticks of the shape of half-cut cylinder and the number of sticks that show belly sides, i.e., the flat sides, determines the number of steps the team's piece can advance on the board. It is possible to pile up one team's pieces if their sites are identical so that pieces as a group can move together afterwards (piling). If a piece of the opponent team is at the new site of one team's piece, the piece is caught and removed from the board, and the team is given one more chance to throw sticks and proceed (catching). For a simplicity, we simulate this game on one-dimensional board with the same number of sites as the original game, and show that catching is more advantageous strategy than piling to win. We also study the avalanche-size distribution in thermodynamic limit to find that it follows an exponential form.

preprint2013arXiv

Fractality of profit landscapes and validation of time series models for stock prices

We apply a simple trading strategy for various time series of real and artificial stock prices to understand the origin of fractality observed in the resulting profit landscapes. The strategy contains only two parameters $p$ and $q$, and the sell (buy) decision is made when the log return is larger (smaller) than $p$ ($-q$). We discretize the unit square $(p, q) \in [0, 1] \times [0, 1]$ into the $N \times N$ square grid and the profit $Π(p, q)$ is calculated at the center of each cell. We confirm the previous finding that local maxima in profit landscapes are scattered in a fractal-like fashion: The number M of local maxima follows the power-law form $M \sim N^{a}$, but the scaling exponent $a$ is found to differ for different time series. From comparisons of real and artificial stock prices, we find that the fat-tailed return distribution is closely related to the exponent $a \approx 1.6$ observed for real stock markets. We suggest that the fractality of profit landscape characterized by $a \approx 1.6$ can be a useful measure to validate time series model for stock prices.

preprint2013arXiv

Nonequilbrium work by charge control in a Josephson junction

We consider a single Josephson junction in the presence of time varying gate charge, and examine the nonequilibrium work done by the charge control in the framework of fluctuation theorems. We obtain the probability distribution functions of the works performed by forward protocol and by its time reversed protocol, which from the Crooks relation gives the estimation of the free energy changes ΔF =0. The reliability of ΔF estimated from the Jarzynksi equality is crucially dependent on protocol parameters, while Bennett's acceptance ratio method confirms consistently ΔF=0. The error of the Jarzynski estimator either grows or becomes saturated as the duration of the work protocol increases, which depends on the protocol rapidity determining the existence of the oscillatory motion of the phase difference across the junction. The average of the work also shows similar behaviors and its saturation value is given by the relative weight of the oscillatory trajectory with respect to running trajectories with constant acceleration. We discuss non-negativity of the work average and its relation to heat and entropy production associated with the circuit control.

preprint2013arXiv

Percolation properties of growing networks under an Achlioptas process

We study the percolation transition in growing networks under an Achlioptas process (AP). At each time step, a node is added in the network and, with the probability $δ$, a link is formed between two nodes chosen by an AP. We find that there occurs the percolation transition with varying $δ$ and the critical point $δ_c=0.5149(1)$ is determined from the power-law behavior of order parameter and the crossing of the fourth-order cumulant at the critical point, also confirmed by the movement of the peak positions of the second largest cluster size to the $δ_c$. Using the finite-size scaling analysis, we get $β/\barν=0.20(1)$ and $1/\barν=0.40(1)$, which implies $β\approx 1/2$ and $\barν \approx 5/2$. The Fisher exponent $τ= 2.24(1)$ for the cluster size distribution is obtained and shown to satisfy the hyperscaling relation.

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.

preprint2013arXiv

Thermodynamic arrow of time of quantum projective measurements

We investigate a thermodynamic arrow associated with quantum projective measurements in terms of the Jensen-Shannon divergence between the probability distribution of energy change caused by the measurements and its time reversal counterpart. Two physical quantities appear to govern the asymptotic values of the time asymmetry. For an initial equilibrium ensemble prepared at a high temperature, the energy fluctuations determine the convergence of the time asymmetry approaching zero. At low temperatures, finite survival probability of the ground state limits the time asymmetry to be less than $\ln 2$. We illustrate our results for a concrete system and discuss the fixed point of the time asymmetry in the limit of infinitely repeated projections.

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

Critical Points of Correlated Percolation in a Gravitational Link-adding Network Model

Motivated by the importance of geometric information in real systems, a new model for long-range correlated percolation in link-adding networks is proposed with the connecting probability decaying with a power-law of the distance on the two-dimensional(2D) plane. By overlapping it with Achlioptas process, it serves as a gravity model which can be tuned to facilitate or inhibit the network percolation in a generic view, cover a broad range of thresholds. Moreover, it yields a set of new scaling relations. In the present work, we develop an approach to determine critical points for them by simulating the temporal evolutions of type-I, type-II and type-III links(chosen from both inter-cluster links, an intra-cluster link compared with an inter-cluster one, and both intra-cluster ones, respectively) and corresponding average lengths. Numerical results have revealed objective competition between fractions, average lengths of three types of links, verified the balance happened at critical points. The variation of decay exponents $a$ or transmission radius $R$ always shifts the temporal pace of the evolution, while the steady average lengths and the fractions of links always keep unchanged just as the values in Achlioptas process. Strategy with maximum gravity can keep steady average length, while that with minimum one can surpass it. Without the confinement of transmission range, $\bar{l} \to \infty$ in thermodynamic limit, while $\bar{l}$ does not when with it. However, both mechanisms support critical points. In two-dimensional free space, the relevance of correlated percolation in link-adding process is verified by validation of new scaling relations with various exponent $a$, which violates the scaling law of Weinrib's.

preprint2012arXiv

Double resonance in the infinite-range quantum Ising model

We study quantum resonance behavior of the infinite-range kinetic Ising model at zero temperature. Numerical integration of the time-dependent Schrödinger equation in the presence of an external magnetic field in the $z$ direction is performed at various transverse field strengths $g$. It is revealed that two resonance peaks occur when the energy gap matches the external driving frequency at two distinct values of $g$, one below and the other above the quantum phase transition. From the similar observations already made in classical systems with phase transitions, we propose that the double resonance peaks should be a generic feature of continuous transitions, for both quantum and classical many-body systems.

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

Force correlations in molecular and stochastic dynamics

A molecular gas system in three dimensions is numerically studied by the energy conserving molecular dynamics (MD). The autocorrelation functions for the velocity and the force are computed and the friction coefficient is estimated. From the comparison with the stochastic dynamics (SD) of a Brownian particle, it is shown that the force correlation function in MD is different from the delta-function force correlation in SD in short time scale. However, as the measurement time scale is increased further, the ensemble equivalence between the microcanonical MD and the canonical SD is restored. We also discuss the practical implication of the result.

preprint2012arXiv

Fractal Profit Landscape of the Stock Market

We investigate the structure of the profit landscape obtained from the most basic, fluctuation based, trading strategy applied for the daily stock price data. The strategy is parameterized by only two variables, p and q. Stocks are sold and bought if the log return is bigger than p and less than -q, respectively. Repetition of this simple strategy for a long time gives the profit defined in the underlying two-dimensional parameter space of p and q. It is revealed that the local maxima in the profit landscape are spread in the form of a fractal structure. The fractal structure implies that successful strategies are not localized to any region of the profit landscape and are neither spaced evenly throughout the profit landscape, which makes the optimization notoriously hard and hypersensitive for partial or limited information. The concrete implication of this property is demonstrated by showing that optimization of one stock for future values or other stocks renders worse profit than a strategy that ignores fluctuations, i.e., a long-term buy-and-hold strategy.

preprint2012arXiv

Neutral theory of chemical reaction networks

To what extent do the characteristic features of a chemical reaction network reflect its purpose and function? In general, one argues that correlations between specific features and specific functions are key to understanding a complex structure. However, specific features may sometimes be neutral and uncorrelated with any system-specific purpose, function or causal chain. Such neutral features are caused by chance and randomness. Here we compare two classes of chemical networks: one that has been subjected to biological evolution (the chemical reaction network of metabolism in living cells) and one that has not (the atmospheric planetary chemical reaction networks). Their degree distributions are shown to share the very same neutral system-independent features. The shape of the broad distributions is to a large extent controlled by a single parameter, the network size. From this perspective, there is little difference between atmospheric and metabolic networks; they are just different sizes of the same random assembling network. In other words, the shape of the degree distribution is a neutral characteristic feature and has no functional or evolutionary implications in itself; it is not a matter of life and death.

preprint2012arXiv

Time Reversibility of Quantum Diffusion in Small-world Networks

We study the time-reversal dynamics of a tight-binding electron in the Watts-Strogatz (WS) small-world networks. The localized initial wave packet at time $t=0$ diffuses as time proceeds until the time-reversal operation, together with the momentum perturbation of the strength $η$, is made at the reversal time $T$. The time irreversibility is measured by $I \equiv |Π(t = 2T) - Π(t = 0)|$, where $Π$ is the participation ratio gauging the extendedness of the wavefunction and for convenience, $t$ is measured forward even after the time reversal . When $η= 0$, the time evolution after $T$ makes the wavefunction at $t=2T$ identical to the one at $t=0$, and we find I=0, implying a null irreversibility or a complete reversibility. On the other hand, as $η$ is increased from zero, the reversibility becomes weaker, and we observe enhancement of the irreversibility. We find that $I$ linearly increases with increasing $η$ in the weakly-perturbed region, and that the irreversibility is much stronger in the WS network than in the local regular network.

preprint2012arXiv

Traveling Baseball Players' Problem in Korea

We study the so-called the traveling tournament problem (TTP), to find an optimal tournament schedule. Differently from the original TTP, in which the total travel distance of all the participants is the objective function to minimize, we instead seek to maximize the fairness of the round robin tournament schedule of the Korean Baseball League. The standard deviation of the travel distances of teams is defined as the energy function, and the Metropolis Monte-Carlo method combined with the simulated annealing technique is applied to find the ground state configuration. The resulting tournament schedule is found to satisfy all the constraint rules set by the Korean Baseball Organization, but with drastically increased fairness in traveling distances.

preprint2011arXiv

Antiphase synchronization of two nonidentical pendulums

We numerically study the synchronization of two nonidentical pendulum motions, pivoting on a common movable frame in the point of view of the dynamic phase transition. When the difference in the pendulum lengths is not too large, it is shown that the system settles down into the dynamic state of the antiphase synchronization with the phase difference $π$. We observe that there is a bistable region where either the antiphase synchronized state or the desynchronized state can be stabilized. We also find that there exists a hysteresis effect around the dynamic phase transition as the length difference is adiabatically changed.

preprint2011arXiv

Comment on "Phase transition in a one-dimensional Ising ferromagnet at zero temperature using Glauber dynamics with a synchronous updating mode"

Sznajd-Weron in [Phys. Rev. E {\bf 82}, 031120 (2010)] suggested that the one-dimensional Ising model subject to the zero temperature synchronous Glauber dynamics exhibits a discontinuous phase transition. We show here instead that the phase transition is of a continuous nature and identify critical exponents: $β\approx 0$, $ν\approx 1$, and $z \approx 2$, via a systematic finite-size scaling analysis.

preprint2011arXiv

Continuity of the Explosive Percolation Transition

The explosive percolation problem on the complete graph is investigated via extensive numerical simulations. We obtain the cluster-size distribution at the moment when the cluster size heterogeneity becomes maximum. The distribution is found to be well described by the power-law form with the decay exponent $τ= 2.06(2)$, followed by a hump. We then use the finite-size scaling method to make all the distributions at various system sizes up to $N=2^{37}$ collapse perfectly onto a scaling curve characterized solely by the single exponent $τ$. We also observe that the instant of that collapse converges to a well-defined percolation threshold from below as $N\rightarrow\infty$. Based on these observations, we show that the explosive percolation transition in the model should be continuous, contrary to the widely-spread belief of its discontinuity.

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

Dissolution of traffic jam via additional local interactions

We use a cellular automata approach to numerically investigate traffic flow patterns on a single lane. The free-flow phase (F), the synchronized phase (S), and the jam phase (J) are observed and the transitions among them are studied as the vehicular density $ρ$ is slowly varied. If $ρ$ is decreased from well inside the J phase, the flux $Φ$ follows the lower branch of the hysteresis loop, implying that the adiabatic decrease of $ρ$ is not an efficient way to put the system back into S or F phases. We propose a simple way to help the system to escape out of J phase, which is based on the local information of the velocities of downstream vehicles.

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

Synchronization in interdependent networks

We explore the synchronization behavior in interdependent systems, where the one-dimensional (1D) network (the intranetwork coupling strength $J_{\rm I}$) is ferromagnetically intercoupled (the strength $J$) to the Watts-Strogatz (WS) small-world network (the intranetwork coupling strength $J_{\rm II}$). In the absence of the internetwork coupling ($J = 0$), the former network is well known not to exhibit the synchronized phase at any finite coupling strength, whereas the latter displays the mean-field transition. Through an analytic approach based on the mean-field approximation, it is found that for the weakly coupled 1D network ($J_{\rm I} \ll 1$) the increase of $J$ suppresses synchrony, because the nonsynchronized 1D network becomes a heavier burden for the synchronization process of the WS network. As the coupling in the 1D network becomes stronger, it is revealed by the renormalization group (RG) argument that the synchronization is enhanced as $J_{\rm I}$ is increased, implying that the more enhanced partial synchronization in the 1D network makes the burden lighter. Extensive numerical simulations confirm these expected behaviors, while exhibiting a reentrant behavior in the intermediate range of $J_{\rm I}$. The nonmonotonic change of the critical value of $J_{\rm II}$ is also compared with the result from the numerical RG calculation.

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.

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

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.

preprint2010arXiv

The Sensitivity of Respondent-driven Sampling Method

Researchers in many scientific fields make inferences from individuals to larger groups. For many groups however, there is no list of members from which to take a random sample. Respondent-driven sampling (RDS) is a relatively new sampling methodology that circumvents this difficulty by using the social networks of the groups under study. The RDS method has been shown to provide unbiased estimates of population proportions given certain conditions. The method is now widely used in the study of HIV-related high-risk populations globally. In this paper, we test the RDS methodology by simulating RDS studies on the social networks of a large LGBT web community. The robustness of the RDS method is tested by violating, one by one, the conditions under which the method provides unbiased estimates. Results reveal that the risk of bias is large if networks are directed, or respondents choose to invite persons based on characteristics that are correlated with the study outcomes. If these two problems are absent, the RDS method shows strong resistance to low response rates and certain errors in the participants' reporting of their network sizes. Other issues that might affect the RDS estimates, such as the method for choosing initial participants, the maximum number of recruitments per participant, sampling with or without replacement and variations in network structures, are also simulated and discussed.

preprint2010arXiv

Voter model on a directed network: Role of bidirectional opinion exchanges

The voter model with the node update rule is numerically investigated on a directed network. We start from a directed hierarchical tree, and split and rewire each incoming arc at the probability $p$. In order to discriminate the better and worse opinions, we break the $Z_2$ symmetry ($σ= \pm 1$) by giving a little more preference to the opinion $σ= 1$. It is found that as $p$ becomes larger, introducing more complicated pattern of information flow channels, and as the network size $N$ becomes larger, the system eventually evolves to the state in which more voters agree on the better opinion, even though the voter at the top of the hierarchy keeps the worse opinion. We also find that the pure hierarchical tree makes opinion agreement very fast, while the final absorbing state can easily be influenced by voters at the higher ranks. On the other hand, although the ordering occurs much slower, the existence of complicated pattern of bidirectional information flow allows the system to agree on the better opinion.

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.

preprint2009arXiv

Heterogeneous attachment strategies optimize the topology of dynamic wireless networks

In optimizing the topology of wireless networks built of a dynamic set of spatially embedded agents, there are many trade-offs to be dealt with. The network should preferably be as small (in the sense that the average, or maximal, pathlength is short) as possible, it should be robust to failures, not consume too much power, and so on. In this paper, we investigate simple models of how agents can choose their neighbors in such an environment. In our model of attachment, we can tune from one situation where agents prefer to attach to others in closest proximity, to a situation where distance is ignored (and thus attachments can be made to agents further away). We evaluate this scenario with several performance measures and find that the optimal topologies, for most of the quantities, is obtained for strategies resulting in a mix of most local and a few random connections.

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.

preprint2007arXiv

A High Robustness and Low Cost Model for Cascading Failures

We study numerically the cascading failure problem by using artificially created scale-free networks and the real network structure of the power grid. The capacity for a vertex is assigned as a monotonically increasing function of the load (or the betweenness centrality). Through the use of a simple functional form with two free parameters, revealed is that it is indeed possible to make networks more robust while spending less cost. We suggest that our method to prevent cascade by protecting less vertices is particularly important for the design of more robust real-world networks to cascading failures.

preprint2007arXiv

Scale-freeness for networks as a degenerate ground state: A Hamiltonian formulation

The origin of scale-free degree distributions in the context of networks is addressed through an analogous non-network model in which the node degree corresponds to the number of balls in a box and the rewiring of links to balls moving between the boxes. A statistical mechanical formulation is presented and the corresponding Hamiltonian is derived. The energy, the entropy, as well as the degree distribution and its fluctuations are investigated at various temperatures. The scale-free distribution is shown to correspond to the degenerate ground state, which has small fluctuations in the degree distribution and yet a large entropy. We suggest an implication of our results from the viewpoint of the stability in evolution of networks.