Source author record

B. Kahng

B. Kahng 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

43works
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

43 published item(s)

preprint2026arXiv

AI-Driven Stabilization in Power Grids through Controlling Line Admittances

The global transition from traditional power plants to renewable energy sources introduces new challenges in grid stability, primarily because inverter-based technologies provide insufficient inertia. To address this, we introduce an artificial intelligence algorithm that autonomously stabilizes power grids by adaptively tuning admittance regulators in response to disturbances. This Adaptive Admittance Controller (AAC) algorithm not only stabilizes the system in real time but also identifies the best regulator locations, thereby unifying grid planning and real time control within a single framework. When tested on a real UK power grid, the AAC markedly reduces frequency deviations and rapidly restores nominal operation. In addition, the algorithm isolates a small number of key regulators and intervenes only on these, lowering both system complexity and cost. The AAC algorithm further reduces the nonlinearity effect, quickly stabilizing the frequency and power flow. This intelligent control scheme enables power grids to reliably return to stable operating conditions under a broad spectrum of fault scenarios. The proposed framework can also be used to mitigate cascading failures by adaptively controlling critical links in a variety of networked infrastructures, such as cascades of traffic congestion on road networks or fuse failures in energy-saving systems.

preprint2020arXiv

Competing synchronization on random networks

The synchronization pattern of a fully connected competing Kuramoto model with a uniform intrinsic frequency distribution $g(ω)$ was recently considered. This competing Kuramoto model assigns two coupling constants with opposite signs, $K_1 < 0$ and $K_2 > 0$, to the $1-p$ and $p$ fractions of nodes, respectively. This model has a rich phase diagram that includes incoherent, $π$, and traveling wave (TW) phases and a hybrid phase transition with abnormal properties that occurs through an intermediate metastable $π$ state. Here, we consider the competing Kuramoto model on Erdős--Rényi (ER) random networks. Numerical simulations and the mean-field solution based on the annealed network approximation reveal that in this case, when the mean degree of the random networks is large, the features of the phase diagram and transition types are consistent overall with those on completely connected networks. However, when the mean degree is small, the mean-field solution is not consistent with the numerical simulation results; specifically, the TW state does not occur, and thus the phase diagram is changed, owing to the strong heterogeneity of the local environment. By contrast, for the original Kuramoto oscillators, the annealed mean-field solution is consistent with the numerical simulation result for ER networks.

preprint2020arXiv

Covid-19 epidemic under the K-quarantine model: Network approach

The Covid-19 pandemic is ongoing worldwide, and the damage it has caused is unprecedented. For prevention, South Korea has adopted a local quarantine strategy rather than a global lockdown. This approach not only minimizes economic damage, but it also efficiently prevents the spread of the disease. In this work, the spread of COVID-19 under local quarantine measures is modeled using the Susceptible-Exposed-Infected-Recovered model on complex networks. In this network approach, the links connected to isolated people are disconnected and then reinstated when they are released. This link dynamics leads to time-dependent reproduction number. Numerical simulations are performed on networks with reaction rates estimated from empirical data. The temporal pattern of the cumulative number of confirmed cases is then reproduced. The results show that a large number of asymptomatic infected patients are detected as they are quarantined together with infected patients. Additionally, possible consequences of the breakdowns of local quarantine measures and social distancing are considered.

preprint2020arXiv

Effective potential approach to hybrid synchronization transitions

The Kuramoto model exhibits different types of synchronization transitions depending on the type of natural frequency distribution. To obtain these results, the Kuramoto self-consistency equation (SCE) approach has been used successfully. However, this approach affords only limited understanding of more detailed properties such as the stability and finite size effect. Here, we extend the SCE approach by introducing an effective potential, that is, an integral version of the SCE. We examine the landscape of this effective potential for second-order, first-order, and hybrid synchronization transitions in the thermodynamic limit. In particular, for the hybrid transition, we find that the minimum of effective potential displays a plateau across the region in which the order parameter jumps. This result suggests that the effective free energy can be used to determine a type of synchronization transition. For finite systems, the effective potential contains local minima at which the system can be trapped. Using numerical simulations, we determine the stability of the system as a function of system size and simulation time.

preprint2020arXiv

Hysteresis and criticality in hybrid percolation transitions

Phase transitions (PTs) are generally classified into second-order and first-order transitions, each exhibiting different intrinsic properties. For instance, a first-order transition exhibits latent heat and hysteresis when a control parameter is increased and then decreased across a transition point, whereas a second-order transition does not. Recently, hybrid percolation transitions (HPTs) are issued in diverse complex systems, in which the features of first-order and second-order PTs occur at the same transition point. Thus, the question whether hysteresis appears in an HPT arises. Herein, we investigate this fundamental question with a so-called restricted Erdős--Rényi random network model, in which a cluster fragmentation process is additionally proposed. The hysteresis curve of the order parameter was obtained. Depending on when the reverse process is initiated, the shapes of hysteresis curves change, and the critical behavior of the HPT is conserved throughout the forward and reverse processes.

preprint2020arXiv

Synchronization in leader-follower switching dynamics

The features of animal population dynamics, for instance, flocking and migration, are often synchronized for survival under large-scale climate change or perceived threats. These coherent phenomena have been explained using synchronization models. However, such models do not take into account asynchronous and adaptive updating of an individual's status at each time. Here, we modify the Kuramoto model slightly by classifying oscillators as leaders or followers, according to their angular velocity at each time, where individuals interact asymmetrically according to their leader/follower status. As the angular velocities of the oscillators are updated, the leader and follower status may also be reassigned. Owing to this adaptive dynamics, oscillators may cooperate by taking turns acting as a leader or follower. This may result in intriguing patterns of synchronization transitions, including hybrid phase transitions, and produce the leader-follower switching pattern observed in bird migration patterns.

preprint2020arXiv

Tricritical directed percolation with long-range interaction in one and two dimensions

Recently, the quantum contact process, in which branching and coagulation processes occur both coherently and incoherently, was theoretically and experimentally investigated in driven open quantum spin systems. In the semi-classical approach, the quantum coherence effect was regarded as a process in which two consecutive atoms are involved in the excitation of a neighboring atom from the inactive (ground) state to the active state (excited $s$ state). In this case, both second-order and first-order transitions occur. Therefore, a tricritical point exists at which the transition belongs to the tricritical directed percolation (TDP) class. On the other hand, when an atom is excited to the $d$ state, long-range interaction is induced. Here, to account for this long-range interaction, we extend the TDP model to one with long-range interaction in the form of $\sim 1/r^{d+σ}$ (denoted as LTDP), where $r$ is the separation, $d$ is the spatial dimension, and $σ$ is a control parameter. In particular, we investigate the properties of the LTDP class below the upper critical dimension $d_c=$ min$(3,\,1.5σ)$. We numerically obtain a set of critical exponents in the LTDP class and determine the interval of $σ$ for the LTDP class. Finally, we construct a diagram of universality classes in the space ($d$, $σ$).

preprint2019arXiv

Entropy production and fluctuation theorems on complex networks

Entropy production (EP) is known as a fundamental quantity for measuring the irreversibility of processes in thermal equilibrium and states far from equilibrium. In stochastic thermodynamics, the EP becomes more visible in terms of the probability density functions of the trajectories of a particle in the state space. Inspired by a previous result that complex networks can serve as state spaces, we consider a data packet transport problem on complex networks. Entropy is produced owing to the complexity of pathways as the packet travels back and forth between two nodes. The EPs are exactly enumerated along the shortest paths between every pair of nodes, and the functional form of the EP distribution is determined by extreme value analysis. The asymptote of the accumulated EP distribution is found to follow the Gumbel distribution.

preprint2019arXiv

Interevent time distribution, burst, and hybrid percolation transition

Critical phenomena of a second-order percolation transition are known to be independent of cluster merging or pruning process. However, those of a hybrid percolation transition (HPT), mixed properties of both first-order and second-order transitions, depend on the processes. The HPT induced by cluster merging is more intrigue and little understood than the other. Here, we construct a theoretical framework using the so-called restricted percolation model. In this model, clusters are ranked by size and partitioned into small- and large-cluster sets. As the cluster rankings are updated by cluster coalescence, clusters may move back and forth across the set boundary. The inter-event time (IET) between two consecutive crossing times have two distributions with power-law decays, which in turn characterize the criticality of the HPT. A burst of such crossing events occurs and signals the upcoming transition. We discuss a related phenomenon to this critical dynamics.

preprint2019arXiv

Simplicial SIS model in scale-free uniform hypergraph

The hypergraph offers a platform to study structural properties emerging from more complicated and higher-order than pairwise interactions among constituents and dynamical behavior such as the spread of information or disease. Recently, a simplicial contagion problem was introduced and considered using a simplicial susceptible-infected-susceptible (SIS) model. Although recent studies have investigated random hypergraphs with a Poisson-type facet degree distribution, hypergraphs in the real world can have a power-law type of facet degree distribution. Here, we consider the SIS contagion problem on scale-free uniform hypergraphs and find that a continuous or hybrid epidemic transition occurs when the hub effect is dominant or weak, respectively. We determine the critical exponents analytically and numerically. We discuss the underlying mechanism of the hybrid epidemic transition.

preprint2016arXiv

Critical behavior of $k$-core percolation: Numerical studies

$k$-Core percolation has served as a paradigmatic model of discontinuous percolation for a long time. Recently it was revealed that the order parameter of $k$-core percolation of random networks additionally exhibits critical behavior. Thus $k$-core percolation exhibits a hybrid phase transition. Unlike the critical behaviors of ordinary percolation that are well understood, those of hybrid percolation transitions have not been thoroughly understood yet. Here, we investigate the critical behavior of $k$-core percolation of Erdős-Rényi networks. We find numerically that the fluctuations of the order parameter and the mean avalanche size diverge in different ways. Thus, we classify the critical exponents into two types: those associated with the order parameter and those with finite avalanches. The conventional scaling relations hold within each set, however, these two critical exponents are coupled. Finally we discuss some universal features of the critical behaviors of $k$-core percolation and the cascade failure model on multiplex networks.

preprint2016arXiv

Genuine Non-Self-Averaging and Ultra-Slow Convergence in Gelation

In irreversible aggregation processes droplets or polymers of microscopic size successively coalesce until a large cluster of macroscopic scale forms. This gelation transition is widely believed to be self-averaging, meaning that the order parameter (the relative size of the largest connected cluster) attains well-defined values upon ensemble averaging with no sample-to-sample fluctuations in the thermodynamic limit. Here, we report on anomalous gelation transition types. Depending on the growth rate of the largest clusters, the gelation transition can show very diverse patterns as a function of the control parameter, which includes multiple stochastic discontinuous transitions, genuine non-self-averaging and ultra-slow convergence of the transition point. Our framework may be helpful in understanding and controlling gelation.

preprint2016arXiv

Hybrid Percolation Transition in Cluster Merging Processes: Continuously Varying Exponents

Consider growing a network, in which every new connection is made between two disconnected nodes. At least one node is chosen randomly from a subset consisting of $g$ fraction of the entire population in the smallest clusters. Here we show that this simple strategy for improving connection exhibits a phase transition barely studied before, namely a hybrid percolation transition exhibiting the properties of both first-order and second-order phase transitions. The cluster size distribution of finite clusters at a transition point exhibits power-law behavior with a continuously varying exponent $τ$ in the range $2 < τ(g) \le 2.5$. This pattern reveals a necessary condition for a hybrid transition in cluster aggregation processes, which is comparable to the power-law behavior of the avalanche size distribution arising in models with link-deleting processes in interdependent networks.

preprint2016arXiv

Hybrid Phase Transition into an Absorbing State: Percolation and Avalanches

Interdependent networks are more fragile under random attacks than simplex networks, because interlayer dependencies lead to cascading failures and finally to a sudden collapse. This is a hybrid phase transition (HPT), meaning that at the transition point the order parameter has a jump but there are also critical phenomena related to it. Here we study these phenomena on the Erdős--Rényi and the two dimensional interdependent networks and show that the hybrid percolation transition exhibits two kinds of critical behaviors: divergence of the fluctuations of the order parameter and power-law size distribution of finite avalanches at a transition point. At the transition point, avalanches of infinite size occur thus the avalanche statistics also has the nature of a HPT. The exponent $β_m$ of the order parameter is $1/2$ under general conditions, while the value of the exponent $γ_m$ characterizing the fluctuations of the order parameter depends on the system. The critical behavior of the finite avalanches can be described by another set of exponents, $β_a$ and $γ_a$. These two critical behaviors are coupled by a scaling law: $1-β_m=γ_a$.

preprint2015arXiv

Efficient algorithm to compute mutually connected components in interdependent networks

Mutually connected components (MCCs) play an important role as a measure of resilience in the study of interdependent networks. Despite their importance, an efficient algorithm to obtain the statistics of all MCCs during the removal of links has thus far been absent. Here, using a well-known fully dynamic graph algorithm, we propose an efficient algorithm to accomplish this task. We show that the time complexity of this algorithm is approximately $O({N^{1.2} })$ for random graphs, which is more efficient than $O(N^{2})$ of the brute-force algorithm. We confirm the correctness of our algorithm by comparing the behavior of the order parameter as links are removed with existing results for three types of double-layer multiplex networks. We anticipate that this algorithm will be used for simulations of large-size systems that have been previously inaccessible.

preprint2015arXiv

Forest Fire Model as a Supercritical Dynamic Model in Financial Systems

Recently, large-scale cascading failures in complex systems have garnered substantial attention. Such extreme events have been treated as an integral part of the self-organized criticality (SOC). Recent empirical work has suggested that some extreme events systematically deviate from the SOC paradigm, requiring a different theoretical framework. We shed additional theoretical light on this possibility by studying financial crisis. We build our model of financial crisis on the well-known forest fire model in scale-free networks. Our analysis shows a non-trivial scaling feature indicating supercritical behavior, which is independent of system size. Extreme events in the supercritical state result from bursting of a fat bubble, seeds of which are sown by a protracted period of a benign financial environment with few shocks. Our findings suggest that policymakers can control the magnitude of financial meltdowns by keeping the economy operating within reasonable duration of a benign environment.

preprint2015arXiv

Two Types of Discontinuous Percolation Transitions in Cluster Merging Processes

Percolation is a paradigmatic model in disordered systems and has been applied to various natural phenomena. The percolation transition is known as one of the most robust continuous transitions. However, recent extensive studies have revealed that a few models exhibit a discontinuous percolation transition (DPT) in cluster merging processes. Unlike the case of continuous transitions, understanding the nature of discontinuous phase transitions requires a detailed study of the system at hand, which has not been undertaken yet for DPTs. Here we examine the cluster size distribution immediately before an abrupt increase in the order parameter of DPT models and find that DPTs induced by cluster merging kinetics can be classified into two types. Moreover, the type of DPT can be determined by the key characteristic of whether the cluster kinetic rule is homogeneous with respect to the cluster sizes. We also establish the necessary conditions for each type of DPT, which can be used effectively when the discontinuity of the order parameter is ambiguous, as in the explosive percolation model.

preprint2014arXiv

Crossover behavior of conductivity in a discontinuous percolation model

When conducting bonds are occupied randomly in a two-dimensional square lattice, the conductivity of the system increases continuously as the density of those conducting bonds exceeds the percolation threshold. Such a behavior is well known in percolation theory; however, the conductivity behavior has not been studied yet when the percolation transition is discontinuous. Here we investigate the conductivity behavior through a discontinuous percolation model evolving under a suppressive external bias. Using effective medium theory, we analytically calculate the conductivity behavior as a function of the density of conducting bonds. The conductivity function exhibits a crossover behavior from a drastically to a smoothly increasing function beyond the percolation threshold in the thermodynamic limit. The analytic expression fits well our simulation data.

preprint2014arXiv

Fast Algorithm for Relaxation Processes in Big-data Systems

Relaxation processes driven by a Laplacian matrix can be found in many real-world big-data systems, for example, in search engines on the World-Wide-Web and the dynamic load balancing protocols in mesh networks. To numerically implement such processes, a fast-running algorithm for the calculation of the pseudo inverse of the Laplacian matrix is essential. Here we propose an algorithm which computes fast and efficiently the pseudo inverse of Markov chain generator matrices satisfying the detailed-balance condition, a general class of matrices including the Laplacian. The algorithm utilizes the renormalization of the Gaussian integral. In addition to its applicability to a wide range of problems, the algorithm outperforms other algorithms in its ability to compute within a manageable computing time arbitrary elements of the pseudo inverse of a matrix of size millions by millions. Therefore our algorithm can be used very widely in analyzing the relaxation processes occurring on large-scale networked systems.

preprint2014arXiv

Origin of the hub spectral dimension in scale-free networks

The return-to-origin probability and the first passage time distribution are essential quantities for understanding transport phenomena in diverse systems. The behaviors of these quantities typically depend on the spectral dimension $d_s$. However, it was recently revealed that in scale-free networks these quantities show a crossover between two power-law regimes characterized by $ d_s $ and the so-called hub spectral dimension $d_s^{\textrm{(hub)}}$ due to the heterogeneity of connectivities of each node. To understand the origin of $d_s^{\textrm{(hub)}}$ from a theoretical perspective, we study a random walk problem on hierarchical scale-free networks by using the renormalization group (RG) approach. Under the RG transformation, not only the system size but also the degree of each node changes due to the scale-free nature of the degree distribution. We show that the anomalous behavior of random walks involving the hub spectral dimension $d_s^{\textrm{(hub)}}$ is induced by the conservation of the power-law degree distribution under the RG transformation.

preprint2014arXiv

Origin of the mixed-order transition in multiplex networks: the Ashkin-Teller model

Recently, diverse phase transition (PT) types have been obtained in multiplex networks, such as discontinuous, continuous, and mixed-order PTs. However, they emerge from individual systems, and there is no theoretical understanding of such PTs in a single framework. Here, we study a spin model called the Ashkin-Teller (AT) model in a mono-layer scale-free network; this can be regarded as a model of two species of Ising spin placed on each layer of a double-layer network. The four-spin interaction in the AT model represents the inter-layer interaction in the multiplex network. Diverse PTs emerge depending on the inter-layer coupling strength and network structure. Especially, we find that mixed-order PTs occur at the critical end points. The origin of such behavior is explained in the framework of Landau-Ginzburg theory.

preprint2014arXiv

Recent advances and open challenges in percolation

Percolation is the paradigm for random connectivity and has been one of the most applied statistical models. With simple geometrical rules a transition is obtained which is related to magnetic models. This transition is, in all dimensions, one of the most robust continuous transitions known. We present a very brief overview of more than 60 years of work in this area and discuss several open questions for a variety of models, including classical, explosive, invasion, bootstrap, and correlated percolation.

preprint2013arXiv

Avoiding a Spanning Cluster in Percolation Models

When dynamics in a system proceeds under suppressive external bias, the system can undergo an abrupt phase transition, as it occurs for example in the epidemic spreading. Recently, an explosive percolation (EP) model was introduced in line with such phenomena. The order of the EP transition has not been clarified in a unified framework covering low dimensional systems and the mean-field limit. We introduce a stochastic model, in which a rule for dynamics is designed to avoid the formation of a spanning cluster through competitive selection in Euclidean space. We show by heuristic arguments that, in the thermodynamic limit and depending on a control parameter, the EP transition can be either continuous or discontinuous if $d < d_c$ and is always continuous if $d \geq d_c$, where $d$ is the spatial dimension and $d_c$ the upper critical dimension.

preprint2012arXiv

Branching process approach for Boolean bipartite networks of metabolic reactions

The branching process (BP) approach has been successful in explaining the avalanche dynamics in complex networks. However, its applications are mainly focused on unipartite networks, in which all nodes are of the same type. Here, motivated by a need to understand avalanche dynamics in metabolic networks, we extend the BP approach to a particular bipartite network composed of Boolean AND and OR logic gates. We reduce the bipartite network into a unipartite network by integrating out OR gates, and obtain the effective branching ratio for the remaining AND gates. Then the standard BP approach is applied to the reduced network, and the avalanche size distribution is obtained. We test the BP results with simulations on the model networks and two microbial metabolic networks, demonstrating the usefulness of the BP approach.

preprint2012arXiv

Discontinuous percolation in diffusion-limited cluster aggregation

Recently, the diffusion-limited cluster aggregation (DLCA) model was restudied as a real-world example of showing discontinuous percolation transitions (PTs). Because a larger cluster is less mobile in Brownian motion, it comes into contact with other clusters less frequently. Thus, the formation of a giant cluster is suppressed in the DLCA process. All clusters grow continuously with respect to time, but the largest cluster grows drastically with respect to the number of cluster merging events. Here, we study the discontinuous PT occurring in the DLCA model in more general dimensions such as two, three, and four dimensions. PTs are also studied for a generalized velocity, which scales with cluster size $s$ as $v_{s} \propto s^η$. For Brownian motion of hard spheres in three dimensions, the mean relative speed scales as $s^{-1/2}$ and the collision rate $σv_s$ scales as $\sim s^{1/6}$. We find numerically that the PT type changes from discontinuous to continuous as $η$ crosses over a tricritical point $η_{c} \approx 1.2$ (in two dimensions), $η_{c} \approx 0.8$ (in three dimensions), and $η_{c} \approx 0.4$ (in four dimensions). We illustrate the root of this crossover behavior from the perspective of the heterogeneity of cluster-size distribution. Finally, we study the reaction-limited cluster aggregation (RLCA) model in the Brownian process, in which cluster merging takes place with finite probability $r$. We find that the PTs in two and three dimensions are discontinuous even for small $r$ such as $r=10^{-3}$, but are continuous in four dimensions.

preprint2012arXiv

First passage time for random walks in heterogeneous networks

The first passage time (FPT) for random walks is a key indicator of how fast information diffuses in a given system. Despite the role of FPT as a fundamental feature in transport phenomena, its behavior, particularly in heterogeneous networks, is not yet fully understood. Here, we study, both analytically and numerically, the scaling behavior of the FPT distribution to a given target node, averaged over all starting nodes. We find that random walks arrive quickly at a local hub, and therefore, the FPT distribution shows a crossover with respect to time from fast decay behavior (induced from the attractive effect to the hub) to slow decay behavior (caused by the exploring of the entire system). Moreover, the mean FPT is independent of the degree of the target node in the case of compact exploration. These theoretical results justify the necessity of using a random jump protocol (empirically used in search engines) and provide guidelines for designing an effective network to make information quickly accessible.

preprint2011arXiv

Discontinuous percolation transitions in real physical systems

We study discontinuous percolation transitions (PT) in the diffusion-limited cluster aggregation model of the sol-gel transition as an example of real physical systems, in which the number of aggregation events is regarded as the number of bonds occupied in the system. When particles are Brownian, in which cluster velocity depends on cluster size as $v_s \sim s^η$ with $η=-0.5$, a larger cluster has less probability to collide with other clusters because of its smaller mobility. Thus, the cluster is effectively more suppressed in growth of its size. Then the giant cluster size increases drastically by merging those suppressed clusters near the percolation threshold, exhibiting a discontinuous PT. We also study the tricritical behavior by controlling the parameter $η$, and the tricritical point is determined by introducing an asymmetric Smoluchowski equation.

preprint2011arXiv

Suppression effect on explosive percolations

When a group of people unknown to each other meet and familiarize among themselves, over time they form a community on a macroscopic scale. This phenomenon can be understood in the context of percolation transition (PT) of networks, which takes place continuously in the classical random graph model. Recently, a modified model was introduced in which the formation of the community was suppressed. Then the PT occurs explosively at a delayed transition time. Whether the explosive PT is indeed discontinuous or continuous becomes controversial. Here we show that type of PT depends on a detailed dynamic rule. Thus, when the dynamic rule is designed to suppress the growth of overall clusters, then the explosive PT could be discontinuous.

preprint2010arXiv

Cluster aggregation model for discontinuous percolation transition

The evolution of the Erdős-Rényi (ER) network by adding edges can be viewed as a cluster aggregation process. Such ER processes can be described by a rate equation for the evolution of the cluster-size distribution with the connection kernel $K_{ij}\sim ij$, where $ij$ is the product of the sizes of two merging clusters. Here, we study more general cases in which $K_{ij}$ is sub-linear as $K_{ij}\sim (ij)^ω$ with $0 \le ω< 1/2$; we find that the percolation transition (PT) is discontinuous. Moreover, PT is also discontinuous when the ER dynamics evolves from proper initial conditions. The rate equation approach for such discontinuous PTs enables us to uncover the mechanism underlying the explosive PT under the Achlioptas process.

preprint2010arXiv

Complete trails of co-authorship network evolution

The rise and fall of a research field is the cumulative outcome of its intrinsic scientific value and social coordination among scientists. The structure of the social component is quantifiable by the social network of researchers linked via co-authorship relations, which can be tracked through digital records. Here, we use such co-authorship data in theoretical physics and study their complete evolutionary trail since inception, with a particular emphasis on the early transient stages. We find that the co-authorship networks evolve through three common major processes in time: the nucleation of small isolated components, the formation of a tree-like giant component through cluster aggregation, and the entanglement of the network by large-scale loops. The giant component is constantly changing yet robust upon link degradations, forming the network's dynamic core. The observed patterns are successfully reproducible through a new network model.

preprint2010arXiv

Diversity and critical behavior in prisoner's dilemma game

The prisoner's dilemma (PD) game is a simple model for understanding cooperative patterns in complex systems consisting of selfish individuals. Here, we study a PD game problem in scale-free networks containing hierarchically organized modules and controllable shortcuts connecting separated hubs. We find that cooperator clusters exhibit a percolation transition in the parameter space (p,b), where p is the occupation probability of shortcuts and b is the temptation payoff in the PD game. The cluster size distribution follows a power law at the transition point. Such a critical behavior, resulting from the combined effect of stochastic processes in the PD game and the heterogeneous structure of complex networks, illustrates the diversity of social relationships and the self-organization of cooperator communities in real-world systems.

preprint2010arXiv

Finite-size scaling theory for explosive percolation transitions

The finite-size scaling (FSS) theory for continuous phase transitions has been useful in determining the critical behavior from the size dependent behaviors of thermodynamic quantities. When the phase transition is discontinuous, however, FSS approach has not been well established yet. Here, we develop a FSS theory for the explosive percolation transition arising in the Erdős and Rényi model under the Achlioptas process. A scaling function is derived based on the observed fact that the derivative of the curve of the order parameter at the critical point $t_c$ diverges with system size in a power-law manner, which is different from the conventional one based on the divergence of the correlation length at $t_c$. We show that the susceptibility is also described in the same scaling form. Numerical simulation data for different system sizes are well collapsed on the respective scaling functions.

preprint2010arXiv

Reduction of high reset currents in unipolar resistance switching Pt/SrTiOx/Pt capacitors using acceptor doping

The high reset current, IR, in unipolar resistance switching is an important issue which should be resolved for practical applications in nonvolatile memories. We showed that,during the forming and set processes, the compliance current, Icomp, can work as a crucial parameter to reduce IR. Doping with Co or Mn can significantly reduce the leakage current in capacitors made using SrTiOx film, opening a larger operation window for Icomp. By decreasing Icomp with acceptor doping, we could reduce IR in SrTiOx films by a factor of approximately 20. Our work suggests that the decrease of Icomp by carrier doping could be a viable alternative for reducing IR in unipolar resistance switching.

preprint2010arXiv

Spectral dimensions of hierarchical scale-free networks with shortcuts

The spectral dimension has been widely used to understand transport properties on regular and fractal lattices. Nevertheless, it has been little studied for complex networks such as scale-free and small world networks. Here we study the spectral dimension and the return-to-origin probability of random walks on hierarchical scale-free networks, which can be either fractals or non-fractals depending on the weight of shortcuts. Applying the renormalization group (RG) approach to the Gaussian model, we obtain the spectral dimension exactly. While the spectral dimension varies between $1$ and $2$ for the fractal case, it remains at $2$, independent of the variation of network structure for the non-fractal case. The crossover behavior between the two cases is studied through the RG flow analysis. The analytic results are confirmed by simulation results and their implications for the architecture of complex systems are discussed.

preprint2009arXiv

Disassortativity of random critical branching trees

Random critical branching trees (CBTs) are generated by the multiplicative branching process, where the branching number is determined stochastically, independent of the degree of their ancestor. Here we show analytically that despite this stochastic independence, there exists the degree-degree correlation (DDC) in the CBT and it is disassortative. Moreover, the skeletons of fractal networks, the maximum spanning trees formed by the edge betweenness centrality, behave similarly to the CBT in the DDC. This analytic solution and observation support the argument that the fractal scaling in complex networks originates from the disassortativity in the DDC.

preprint2009arXiv

Large 1/f noise of unipolar resistance switching and its percolating nature

We investigated the 1/f noise of Pt/NiO/Pt capacitors that show unipolar resistance switching. When they were switched from the low to high resistance states, the power spectral density of the voltage fluctuation was increased by approximately five orders of magnitude. At 100 K, the relative resistance fluctuation, SR/R2, in the low resistance state displayed a power law dependence on the resistance R with exponent w = 1.6. This behavior can be explained by percolation theory; however, at higher temperatures or near the switching voltage, SR/R2 becomes enhanced further. This large 1/f noise can be therefore an important problem in the development of resistance random access memory devices.

preprint2009arXiv

Percolation Transitions in Scale-Free Networks under Achlioptas Process

It has been recently shown that the percolation transition is discontinuous in Erdős-Rényi networks and square lattices in two dimensions under the Achlioptas Process (AP). Here, we show that when the structure is highly heterogeneous as in scale-free networks, a discontinuous transition does not always occur: a continuous transition is also possible depending on the degree distribution of the scale-free network. This originates from the competition between the AP that discourages the formation of a giant component and the existence of hubs that encourages it. We also estimate the value of the characteristic degree exponent that separates the two transition types.

preprint2008arXiv

Jamming transition in traffic flow under the priority queuing protocol

Packet traffic in complex networks undergoes the jamming transition from free-flow to congested state as the number of packets in the system increases. Here we study such jamming transition when queues are operated by the priority queuing protocol and packets are guided by the dynamic routing protocol. We introduce a minimal model in which there are two types of packets distinguished by whether priority is assigned. Based on numerical simulations, we show that traffic is improved in the congested region under the priority queuing protocol, and it is worsened in the free-flow region. Also, we find that at the transition point, the waiting-time distribution follows a power law, and the power spectrum of traffic exhibits a crossover between two 1/f^a behaviors with exponent a ~ 1 and 1 < a < 2 in low and high frequency regime, respectively. This crossover is originated from a characteristic waiting time of packets in the queue.

preprint2008arXiv

Segregation in the annihilation of two-species reaction-diffusion processes on fractal scale-free networks

In the reaction-diffusion process $A+B \to \varnothing$ on random scale-free (SF) networks with the degree exponent $γ$, the particle density decays with time in a power law with an exponent $α$ when initial densities of each species are the same. The exponent $α$ is $α> 1$ for $2 < γ< 3$ and $α=1$ for $γ\ge 3$. Here, we examine the reaction process on fractal SF networks, finding that $α< 1$ even for $2 < γ< 3$. This slowly decaying behavior originates from the segregation effect: Fractal SF networks contain local hubs, which are repulsive to each other. Those hubs attract particles and accelerate the reaction, and then create domains containing the same species of particles. It follows that the reaction takes place at the non-hub boundaries between those domains and thus the particle density decays slowly. Since many real SF networks are fractal, the segregation effect has to be taken into account in the reaction kinetics among heterogeneous particles.

preprint2006arXiv

Internet data packet transport: from global topology to local queueing dynamics

We study structural feature and evolution of the Internet at the autonomous systems level. Extracting relevant parameters for the growth dynamics of the Internet topology, we construct a toy model for the Internet evolution, which includes the ingredients of multiplicative stochastic evolution of nodes and edges and adaptive rewiring of edges. The model reproduces successfully structural features of the Internet at a fundamental level. We also introduce a quantity called the load as the capacity of node needed for handling the communication traffic and study its time-dependent behavior at the hubs across years. The load at hub increases with network size $N$ as $\sim N^{1.8}$. Finally, we study data packet traffic in the microscopic scale. The average delay time of data packets in a queueing system is calculated, in particular, when the number of arrival channels is scale-free. We show that when the number of arriving data packets follows a power law distribution, $\sim n^{-λ}$, the queue length distribution decays as $n^{1-λ}$ and the average delay time at the hub diverges as $\sim N^{(3-λ)/(γ-1)}$ in the $N \to \infty$ limit when $2 < λ< 3$, $γ$ being the network degree exponent.

preprint2006arXiv

Synchronization transition of heterogeneously coupled oscillators on scale-free networks

We investigate the synchronization transition of the modified Kuramoto model where the oscillators form a scale-free network with degree exponent $λ$. An oscillator of degree $k_i$ is coupled to its neighboring oscillators with asymmetric and degree-dependent coupling in the form of $\couplingcoeff k_i^{η-1}$. By invoking the mean-field approach, we determine the synchronization transition point $J_c$, which is zero (finite) when $η> λ-2$ ($η< λ-2$). We find eight different synchronization transition behaviors depending on the values of $η$ and $λ$, and derive the critical exponents associated with the order parameter and the finite-size scaling in each case. The synchronization transition is also studied from the perspective of cluster formation of synchronized vertices. The cluster-size distribution and the largest cluster size as a function of the system size are derived for each case using the generating function technique. Our analytic results are confirmed by numerical simulations.

preprint2001arXiv

Robustness of the in-degree exponent for the world-wide web

We consider a stochastic model for directed scale-free networks following power-laws in the degree distributions in both incoming and outgoing directions. In our model, the number of vertices grow geometrically with time with growth rate p. At each time step, (i) each newly introduced vertex is connected to a constant number of already existing vertices with the probability linearly proportional to the in-degree of a selected vertex, and (ii) each existing vertex updates its outgoing edges through a stochastic multiplicative process with mean growth rate of outgoing edges g and variance $σ^2$. Using both analytic treatment and numerical simulations, we show that while the out-degree exponent $γ_{\rm out}$ depends on the parameters, the in-degree exponent $γ_{\rm in}$ has two distinct values, $γ_{\rm in}=2$ for $p > g$ and 1 for $p < g$, independent of different parameters values. The latter case has logarithmic correction to the power-law. Since the vertex growth rate p is larger than the degree growth rate g for the world-wide web (www) nowadays, the in-degree exponent appears robust as $γ_{\rm in}=2$ for the www.

preprint1997arXiv

Island Density in Homoepitaxial Growth:Improved Monte Carlo Results

We reexamine the density of two dimensional islands in the submonolayer regime of a homoepitaxially growing surface using the coarse grained Monte Carlo simulation with random sequential updating rather than parallel updating. It turns out that the power law dependence of the density of islands on the deposition rate agrees much better with the theoretical prediction than previous data obtained by other methods if random sequential instead of parallel updating is used.