Trust snapshot

Quick read

Trust 21 - Emerging
20works
0followers
14topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

20 published item(s)

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

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

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.

preprint2011arXiv

A Paradoxical Property of the Monkey Book

A &#34;monkey book&#34; 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&#39; power law to an extraordinary good approximation, in contrast to the word distributions for sampled and real books, which deviate from Heaps&#39; law in a characteristics way. The somewhat counter-intuitive conclusion is that a &#34;monkey book&#34; obeys Heaps&#39; 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

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

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 &#34;Kim&#34; 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&#39;s law unzipped

Why does Zipf&#39;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&#39;s and Gibrat&#39;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&#39;

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

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

Selective pressure on metabolic network structures as measured from the random blind-watchmaker network

A random null model termed the Blind Watchmaker network (BW) has been shown to reproduce the degree distribution found in metabolic networks. This might suggest that natural selection has had little influence on this particular network property. We here investigate to what extent other structural network properties have evolved under selective pressure from the corresponding ones of the random null model: The clustering coefficient and the assortativity measures are chosen and it is found that these measures for the metabolic network structure are close enough to the BW-network so as to fit inside its reachable random phase space. It is furthermore shown that the use of this null model indicates an evolutionary pressure towards low assortativity and that this pressure is stronger for larger networks. It is also shown that selecting for BW networks with low assortativity causes the BW degree distribution to slightly deviate from its power-law shape in the same way as the metabolic networks. This implies that an equilibrium model with fluctuating degree distribution is more suitable as a null model, when identifying selective pressures, than a randomized counterpart with fixed degree sequence, since the overall degree sequence itself can change under selective pressure on other global network properties.

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&#39;

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.

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.

preprint2004arXiv

Correlations in Networks associated to Preferential Growth

Combinations of random and preferential growth for both on-growing and stationary networks are studied and a hierarchical topology is observed. Thus for real world scale-free networks which do not exhibit hierarchical features preferential growth is probably not the main ingredient in the growth process. An example of such real world networks includes the protein-protein interaction network in yeast, which exhibits pronounced anti-hierarchical features.

preprint1995arXiv

Verification of a New Non-Linear IV-exponent: Simulation of the 2D Coulomb Gas with Langevin Dynamics.

It has recently been suggested from scaling arguments that the non-linear IV-exponent a, for a two-dimensional superconductor is different from the exponent originally suggested by Ambegaokar et al. The relation between the new and the old exponent is a=a_AHNS-3. The new scaling behaviour is linked to the logarithmic vortex interaction and the long range time tail which this gives rise to. Consequently one may expect that the scaling behavior is generic for models which have these basic features. The simplest model of this type is the two-dimensional Coulomb gas model with Langevin dynamics. We here explicitly verify, through computer simulations, that the IV-characteristics of this model indeed scales according to the new scaling exponent a. Keywords: vortex, Coulomb gas, IV-exponent, Simulations, Langevin, 2D superconductor, thin films.