Robustness of the avalanche dynamics in data packet transport on scale-free networks

We study the avalanche dynamics in the data packet transport on scale-free networks through a simple model. In the model, each vertex is assigned a capacity proportional to the load with a proportionality constant $1+a$. When the system is perturbed by a single vertex removal, the load of each vertex is redistributed, followed by subsequent failures of overloaded vertices. The avalanche size depends on the parameter $a$ as well as which vertex triggers it. We find that there exists a critical value $a_c$ at which the avalanche size distribution follows a power law. The critical exponent associated with it appears to be robust as long as the degree exponent is between 2 and 3, and is close in value to that of the distribution of the diameter changes by single vertex removal.Comment: 5 pages, 7 figures, final version published in PR

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^{-\lambda}$, the queue length distribution decays as $n^{1-\lambda}$ and the average delay time at the hub diverges as $\sim N^{(3-\lambda)/(\gamma-1)}$ in the $N \to \infty$ limit when $2 < \lambda < 3$, $\gamma$ being the network degree exponent.Comment: 5 pages, 4 figures, submitted to International Journal of Bifurcation and Chao

Evolution of the Protein Interaction Network of Budding Yeast: Role of the Protein Family Compatibility Constraint

Understanding of how protein interaction networks (PIN) of living organisms have evolved or are organized can be the first stepping stone in unveiling how life works on a fundamental ground. Here we introduce a hybrid network model composed of the yeast PIN and the protein family interaction network. The essential ingredient of the model includes the protein family identity and its robustness under evolution, as well as the three previously proposed ones: gene duplication, divergence, and mutation. We investigate diverse structural properties of our model with parameter values relevant to yeast, finding that the model successfully reproduces the empirical data.Comment: 5 pages, 5 figures, 1 table. Title changed. Final version published in JKP

Doppler radar wind field retrieval over the Po Valley

Although methods of using multiple Doppler radars to study wind fields have long been proposed, and many research studies have been made, very few operational radar operators adopt methods which require the use of specific scanning strategies to allow the extraction of wind information. Here we report a collaborative study on dual-Doppler radars based on two Doppler radars in the Po valley, Italy. Unusually, the radars are only about 90 km apart, though operated by the same authority. The wind field syntheses are carried out on a 30 km by 30 km region where the two radars have overlapping scan coverage. An iterative method based on the linear wind model and the equation of mass continuity is used to construct the wind fields. The methodology has been validated by two different methods. The first method is to reconstruct the radial wind observed by each radar, and the second method is calculating and comparing the along-track component with that derived from the observations. Both two comparisons show good agreement with the original data