-
Mapping Inter-City Trade Networks to Maximum Entropy Models using Electronic Invoice Data
Authors:
Cesar I. N. Sampaio Filho,
Rilder S. Pires,
Humberto A. Carmona,
José S. Andrade Jr
Abstract:
We analyze the network of transactions among cities based on the electronic invoice database for the municipalities in the Ceará state, Brazil. It consists of approximately 3.7 billion records, registered during the period between the years 2016 to 2019. All the transactions are grouped in a unique dataset and represented as an asymmetrical adjacency matrix corresponding to a directed graph with c…
▽ More
We analyze the network of transactions among cities based on the electronic invoice database for the municipalities in the Ceará state, Brazil. It consists of approximately 3.7 billion records, registered during the period between the years 2016 to 2019. All the transactions are grouped in a unique dataset and represented as an asymmetrical adjacency matrix corresponding to a directed graph with connections weighted by the number of transactions among cities. Due to the large size of Ceará state, its unequal distribution of wealth, and spatially heterogeneous population, we initially determine communities of cities based on their mutual intensity of trades and verify to which extent their economic interests reflect a community cohesiveness. For the first task, we use the Infomap algorithm to detect the partition which provides the shortest description length and captures the optimal community structure of the network in terms of its associated flow dynamics. Surprisingly, the partition identified has five modules, whose two-dimensional geographical projections are all simply-connected domains. We proceed with the analysis of traded products by building bipartite structures represented in terms of adjacency matrices between municipalities and products, considering selling and buying. Using the revealed comparative advantage (RCA) concept, we define a non-monetary and binary activity index that can distinguish the RCA of a city in a class of goods or services as evidenced by trade flows. Finally, through the pairwise Maximum Entropy Model, we can associate to the largest communities their corresponding binary Ising-like Hamiltonian models. In an analogy with critical phenomena, our results reveal that each community operates at a "temperature" that is close to the corresponding "critical point", suggesting a high degree of "economic cohesiveness" in its trade network of cities.
△ Less
Submitted 18 July, 2024;
originally announced July 2024.
-
Tricritical behavior in epidemic dynamics with vaccination
Authors:
Marcelo A. Pires,
Cesar I. N. Sampaio Filho,
Hans J. Herrmann,
José S. Andrade Jr
Abstract:
We scrutinize the phenomenology arising from a minimal vaccination-epidemic (MVE) dynamics using three methods: mean-field approach, Monte Carlo simulations, and finite-size scaling analysis. The mean-field formulation reveals that the MVE model exhibits either a continuous or a discontinuous active-to-absorbing phase transition, accompanied by bistability and a tricritical point. However, on squa…
▽ More
We scrutinize the phenomenology arising from a minimal vaccination-epidemic (MVE) dynamics using three methods: mean-field approach, Monte Carlo simulations, and finite-size scaling analysis. The mean-field formulation reveals that the MVE model exhibits either a continuous or a discontinuous active-to-absorbing phase transition, accompanied by bistability and a tricritical point. However, on square lattices, we detect no signs of bistability, and we disclose that the active-to-absorbing state transition has a scaling invariance and critical exponents compatible with the continuous transition of the directed percolation universality class. Additionally, our findings indicate that the tricritical and crossover behaviors of the MVE dynamics belong to the universality class of mean-field tricritical directed percolation.
△ Less
Submitted 26 November, 2023; v1 submitted 13 June, 2023;
originally announced June 2023.
-
Random networks with q-exponential degree distribution
Authors:
Cesar I. N. Sampaio Filho,
Marcio M. Bastos,
Hans J. Herrmann,
André A. Moreira,
José S. Andrade Jr
Abstract:
We use the configuration model to generate networks having a degree distribution that follows a $q$-exponential, $P_q(k)=(2-q)λ[1-(1-q)λk]^{1/(q-1)}$, for arbitrary values of the parameters $q$ and $λ$. We study the assortativity and the shortest path of these networks finding that the more the distribution resembles a pure power law, the less well connected are the corresponding nodes. In fact, t…
▽ More
We use the configuration model to generate networks having a degree distribution that follows a $q$-exponential, $P_q(k)=(2-q)λ[1-(1-q)λk]^{1/(q-1)}$, for arbitrary values of the parameters $q$ and $λ$. We study the assortativity and the shortest path of these networks finding that the more the distribution resembles a pure power law, the less well connected are the corresponding nodes. In fact, the average degree of a nearest neighbor grows monotonically with $λ^{-1}$. Moreover, our results show that $q$-exponential networks are more robust against random failures and against malicious attacks than standard scale-free networks. Indeed, the critical fraction of removed nodes grows logarithmically with $λ^{-1}$ for malicious attacks. An analysis of the $k_s$-core decomposition shows that $q$-exponential networks have a highest $k_s$-core, that is bigger and has a larger $k_s$ than pure scale-free networks. Being at the same time well connected and robust, networks with $q$-exponential degree distribution exhibit scale-free and small-world properties, making them a particularly suitable model for application in several systems.
△ Less
Submitted 26 December, 2022;
originally announced December 2022.
-
The symbiotic contact process: phase transitions, hysteresis cycles, and bistability
Authors:
C. I. N. Sampaio Filho,
T. B. dos Santos,
N. A. M. Araújo,
H. A. Carmona,
A. A. Moreira,
J. S. Andrade Jr
Abstract:
We performed Monte Carlo simulations of the symbiotic contact process on different spatial dimensions ($d$). On the complete and random graphs (infinite dimension), we observe hysteresis cycles and bistable regions, what is consistent with the discontinuous absorbing-state phase transition predicted by mean-field theory. By contrast, on a regular square lattice, we find no signs of bistability or…
▽ More
We performed Monte Carlo simulations of the symbiotic contact process on different spatial dimensions ($d$). On the complete and random graphs (infinite dimension), we observe hysteresis cycles and bistable regions, what is consistent with the discontinuous absorbing-state phase transition predicted by mean-field theory. By contrast, on a regular square lattice, we find no signs of bistability or hysteretic behavior. This result suggests that the transition in two dimensions is rather continuous. Based on our numerical observations, we conjecture that the nature of the transition changes at the upper critical dimension ($d_c$), from continuous ($d<d_c$) to discontinuous ($d>d_c$).
△ Less
Submitted 28 June, 2018;
originally announced June 2018.
-
Elastic backbone defines a new transition in the percolation model
Authors:
Cesar I. N. Sampaio Filho,
José S. Andrade Jr.,
Hans J. Herrmann,
André A. Moreira
Abstract:
The elastic backbone is the set of all shortest paths. We found a new phase transition at $p_{eb}$ above the classical percolation threshold at which the elastic backbone becomes dense. At this transition in $2d$ its fractal dimension is $1.750\pm 0.003$, and one obtains a novel set of critical exponents $β_{eb} = 0.50\pm 0.02$, $γ_{eb} = 1.97\pm 0.05$, and $ν_{eb} = 2.00\pm 0.02$ fulfilling consi…
▽ More
The elastic backbone is the set of all shortest paths. We found a new phase transition at $p_{eb}$ above the classical percolation threshold at which the elastic backbone becomes dense. At this transition in $2d$ its fractal dimension is $1.750\pm 0.003$, and one obtains a novel set of critical exponents $β_{eb} = 0.50\pm 0.02$, $γ_{eb} = 1.97\pm 0.05$, and $ν_{eb} = 2.00\pm 0.02$ fulfilling consistent critical scaling laws. Interestingly, however, the hyperscaling relation is violated. Using Binder's cumulant, we determine, with high precision, the critical probabilities $p_{eb}$ for the triangular and tilted square lattice for site and bond percolation. This transition describes a sudden rigidification as a function of density when stretching a damaged tissue.
△ Less
Submitted 21 March, 2018;
originally announced March 2018.
-
Crossover from mean-field to $2d$ Directed Percolation in the contact process
Authors:
T. B. dos Santos,
C. I. N. Sampaio Filho,
N. A. M. Araújo,
C. L. N. Oliveira,
A. A. Moreira
Abstract:
We study the contact process on spatially embedded networks, consisting of a regular square lattice with long-range connections. To generate the networks, a long-range connection is randomly added to each node $i$ of a square lattice, following the probability, $P_{ij}\sim{r_{ij}^{-α}}$ , where $r_{ij}$ is the Manhattan distance between nodes $i$ and $j$, and the exponent $α$ is a tunable paramete…
▽ More
We study the contact process on spatially embedded networks, consisting of a regular square lattice with long-range connections. To generate the networks, a long-range connection is randomly added to each node $i$ of a square lattice, following the probability, $P_{ij}\sim{r_{ij}^{-α}}$ , where $r_{ij}$ is the Manhattan distance between nodes $i$ and $j$, and the exponent $α$ is a tunable parameter. Extensive Monte Carlo simulations and a finite-size scaling analysis for different values of $α$ reveal a crossover from the mean-field to $2d$ Directed Percolation universality class with increasing $α$, in the range $3<α<4$.
△ Less
Submitted 28 February, 2018; v1 submitted 28 February, 2018;
originally announced February 2018.
-
Noise induced phase transition in the $S$-state block voter model
Authors:
J. M. de Araújo,
C. I. N. Sampaio Filho,
F. G. B. Moreira
Abstract:
We use Monte Carlo simulations and finite-size scaling theory to investigate the phase transition and critical behavior of the $S$-state block voter model on square lattices. It is shown that the system exhibits an order-disorder phase transition at a given value of the noise parameter, which changes from a continuous transition for $S\le4$ to a discontinuous transition for $S\ge5$. Moreover, for…
▽ More
We use Monte Carlo simulations and finite-size scaling theory to investigate the phase transition and critical behavior of the $S$-state block voter model on square lattices. It is shown that the system exhibits an order-disorder phase transition at a given value of the noise parameter, which changes from a continuous transition for $S\le4$ to a discontinuous transition for $S\ge5$. Moreover, for the cases of continuous transition, the calculated critical exponents indicate that the present studied nonequilibrium model system is in the same universality class of its counterpart equilibrium two-dimensional S-state Potts model. We also provide a first estimation of the long-range exponents governing the dependence on the range of interaction of the magnetization, the susceptibility, and the derivative of Binder's cumulant.
△ Less
Submitted 29 May, 2018; v1 submitted 27 September, 2016;
originally announced September 2016.
-
Itinerant conductance in fuse-antifuse networks
Authors:
Cesar I. N. Sampaio Filho,
Andre A. Moreira,
Nuno A. N. Araujo,
Jose S. Andrade,
Hans J. Herrmann
Abstract:
We report on a novel dynamic phase in electrical networks, in which current channels perpetually change in time. This occurs when the individual elements of the network are fuse-antifuse devices, namely, become insulators within a certain finite interval of local applied voltages. As a consequence, the macroscopic current exhibits temporal fluctuations which increase with system size. We determine…
▽ More
We report on a novel dynamic phase in electrical networks, in which current channels perpetually change in time. This occurs when the individual elements of the network are fuse-antifuse devices, namely, become insulators within a certain finite interval of local applied voltages. As a consequence, the macroscopic current exhibits temporal fluctuations which increase with system size. We determine the conditions under which this exotic situation appears by establishing a phase diagram as a function of the applied field and the size of the insulating window. Besides its obvious application as a versatile electronic device, due to its rich variety of behaviors, this network model provides a possible description for particle-laden flow through porous media leading to dynamical clogging and reopening of the local channels in the pore space.
△ Less
Submitted 29 August, 2016;
originally announced August 2016.
-
Majority-vote model on spatially embedded networks: crossover from mean-field to Ising universality classes
Authors:
C. I. N. Sampaio Filho,
T. B. dos Santos,
A. A. Moreira,
F. G. B. Moreira,
J. S. Andrade Jr
Abstract:
We study through Monte Carlo simulations and finite-size scaling analysis the nonequilibrium phase transitions of the majority-vote model taking place on spatially embedded networks. These structures are built from an underlying regular lattice over which long-range connections are randomly added according to the probability, $P_{ij}\sim{r^{-α}}$, where $r_{ij}$ is the Manhattan distance between n…
▽ More
We study through Monte Carlo simulations and finite-size scaling analysis the nonequilibrium phase transitions of the majority-vote model taking place on spatially embedded networks. These structures are built from an underlying regular lattice over which long-range connections are randomly added according to the probability, $P_{ij}\sim{r^{-α}}$, where $r_{ij}$ is the Manhattan distance between nodes $i$ and $j$, and the exponent $α$ is a controlling parameter [J. M. Kleinberg, Nature 406, 845 (2000)]. Our results show that the collective behavior of this system exhibits a continuous order-disorder phase transition at a critical parameter, which is a decreasing function of the exponent $α$. Precisely, considering the scaling functions and the critical exponents calculated, we conclude that the system undergoes a crossover among distinct universality classes. For $α\le3$ the critical behavior is described by mean-field exponents, while for $α\ge4$ it belongs to the Ising universality class. Finally, in the region where the crossover occurs, $3<α<4$, the critical exponents are dependent on $α$.
△ Less
Submitted 29 February, 2016;
originally announced February 2016.
-
Characterizing the intrinsic correlations of scale-free networks
Authors:
J. B. de Brito,
C. I. N. Sampaio Filho,
A. A. Moreira,
J. S. Andrade Jr
Abstract:
Very often, when studying topological or dynamical properties of random scale-free networks, it is tacitly assumed that degree-degree correlations are not present. However, simple constraints, such as the absence of multiple edges and self-loops, can give rise to intrinsic correlations in these structures. In the same way that Fermionic correlations in thermodynamic systems are relevant only in th…
▽ More
Very often, when studying topological or dynamical properties of random scale-free networks, it is tacitly assumed that degree-degree correlations are not present. However, simple constraints, such as the absence of multiple edges and self-loops, can give rise to intrinsic correlations in these structures. In the same way that Fermionic correlations in thermodynamic systems are relevant only in the limit of low temperature, the intrinsic correlations in scale-free networks are relevant only when the extreme values for the degrees grow faster than the square-root of the network size. In this situation, these correlations can significantly affect the dependence of the average degree of the nearest neighbors of a given vertice on this vertices's degree. Here, we introduce an analytical approach that is capable to predict the functional form of this property. Moreover, our results indicate that random scale-free networks models are not self-averaging, that is, the second moment of their degree distribution may vary orders of magnitude among different realizations. Finally, we argue that the intrinsic correlations investigated here may have profound impact on the critical properties of random scale-free networks.
△ Less
Submitted 10 June, 2015;
originally announced June 2015.
-
Mandala Networks: ultra-robust, ultra-small-world and highly sparse graphs
Authors:
C. I. N. Sampaio Filho,
A. A. Moreira,
R. F. S. Andrade,
H. J. Herrmann,
J. S. Andrade Jr
Abstract:
The increasing demands in security and reliability of infrastructures call for the optimal design of their embedded complex networks topologies. The following question then arises: what is the optimal layout to fulfill best all the demands? Here we present a general solution for this problem with scale-free networks, like the Internet and airline networks. Precisely, we disclose a way to systemati…
▽ More
The increasing demands in security and reliability of infrastructures call for the optimal design of their embedded complex networks topologies. The following question then arises: what is the optimal layout to fulfill best all the demands? Here we present a general solution for this problem with scale-free networks, like the Internet and airline networks. Precisely, we disclose a way to systematically construct networks which are 100$\%$ robust against random failures as well as to malicious attacks. Furthermore, as the sizes of these networks increase, their shortest paths become asymptotically invariant and densities of links go to zero, making them ultra-small worlds and highly sparse, respectively. The first property is ideal for communication and navigation purposes, while the second is interesting economically.
△ Less
Submitted 27 June, 2014;
originally announced June 2014.