Showing 1–50 of 88 results for author: Andrade, J

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.

Comments: 29 pages, 10 figures

Fracture and size effect in mechanical metamaterials

Authors: J. Ulloa, M. P. Ariza, J. E. Andrade, M. Ortiz

Abstract: We resort to variational methods to evaluate the asymptotic behavior of fine metamaterials as a function of cell size. To zeroth order, the metamaterial behaves as a micropolar continuum with both displacement and rotation degrees of freedom, but exhibits linear-elastic fracture mechanics scaling and therefore no size effect. To higher order, the overall energetics of the metastructure can be char… ▽ More We resort to variational methods to evaluate the asymptotic behavior of fine metamaterials as a function of cell size. To zeroth order, the metamaterial behaves as a micropolar continuum with both displacement and rotation degrees of freedom, but exhibits linear-elastic fracture mechanics scaling and therefore no size effect. To higher order, the overall energetics of the metastructure can be characterized explicitly in terms of the solution of the zeroth-order continuum problem by the method of Γ-expansion. We present explicit expressions of the second-order correction for octet frames. As an application, we evaluate the compliance of double-cantilever octet specimens to second order and use the result to elucidate the dependence of the apparent toughness of the specimen on cell size. The analysis predicts the discreteness of the metamaterial lattice to effectively shield the crack-tip, a mechanism that we term lattice shielding. The theory specifically predicts anti-shielding, i. e., coarser is weaker, in agreement with recent experimental observations. △ Less

Submitted 26 June, 2024; originally announced July 2024.

Suspended accounts align with the Internet Research Agency misinformation campaign to influence the 2016 US election

Authors: Matteo Serafino, Zhenkun Zhou, Jose S. Andrade, Jr., Alexandre Bovet, Hernan A. Makse

Abstract: The ongoing debate surrounding the impact of the Internet Research Agency s (IRA) social media campaign during the 2016 U.S. presidential election has largely overshadowed the involvement of other actors. Our analysis brings to light a substantial group of suspended Twitter users, outnumbering the IRA user group by a factor of 60, who align with the ideologies of the IRA campaign. Our study demons… ▽ More The ongoing debate surrounding the impact of the Internet Research Agency s (IRA) social media campaign during the 2016 U.S. presidential election has largely overshadowed the involvement of other actors. Our analysis brings to light a substantial group of suspended Twitter users, outnumbering the IRA user group by a factor of 60, who align with the ideologies of the IRA campaign. Our study demonstrates that this group of suspended Twitter accounts significantly influenced individuals categorized as undecided or weak supporters, potentially with the aim of swaying their opinions, as indicated by Granger causality. △ Less

Submitted 13 March, 2024; v1 submitted 28 July, 2023; originally announced July 2023.

Comments: 7 Tables, 5 Figures

Correlation between the rheology of 2D-inks precursors and the droplet size generated from a capillary nozzle in dripping regime

Authors: Pedro C. Rijo, Josué M. O. Cremonezzi, Ricardo J. E. Andrade, Francisco J. Galindo-Rosales

Abstract: This study provides a complete rheological characterization of 2D nanomaterial dispersions, employed as 2D-inks precursors in printed electronics. Three different 2D nanomaterials (molybdenum disulfide (MoS2), graphene, and hexagonal boron nitride(hBN)) were dispersed in a Newtonian fluid (toluene) and a viscoelastic fluid (toluene + ethyl cellulose) with different polymer concentrations. The pres… ▽ More This study provides a complete rheological characterization of 2D nanomaterial dispersions, employed as 2D-inks precursors in printed electronics. Three different 2D nanomaterials (molybdenum disulfide (MoS2), graphene, and hexagonal boron nitride(hBN)) were dispersed in a Newtonian fluid (toluene) and a viscoelastic fluid (toluene + ethyl cellulose) with different polymer concentrations. The presence of nanoparticles does not change the shear rheology of the carrier fluid. Regarding the extensional rheology, the results showed that the pinch-off phenomenon is present in all Toluene suspensions; however, the presence of the ethyl cellulose introduces elasticity in the system, even leading to the formation of beads-on-a-string, and the relaxation times of the suspensions depends on the kind of nanoparticles present in the fluid. As controlling the droplet size when dispensing 2D-inks is of paramount importance for printed electronics, as well as for many other applications, here it is presented a correlation between the rheological properties of these 2D-inks precursors and their droplet size when generated from a capillary nozzle in dripping regime. △ Less

Submitted 29 June, 2023; originally announced June 2023.

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.

Comments: 12 pages, 6 figures and 2 tables. Version closer to the published paper

Report number: 174

Journal ref: Chaos, Solitons & Fractals, 2023

Continuum model for extraction and retention in porous media

Authors: André F. V. Matias, Daniel F. Valente-Matias, Nuno R. Neng, José M. F. Nogueira, José S. Andrade Jr., Rodrigo C. V. Coelho, Nuno A. M. Araújo

Abstract: Several natural and industrial processes involve the extraction or retention of a solute by a fluid invading a network of channels. Examples include aquifer contamination, chemical filtration, and coffee extraction. We propose a continuum equation to model these processes, parametrized by the Péclet number and the rate of mass transfer between the solid and the fluid. We study the time dependence… ▽ More Several natural and industrial processes involve the extraction or retention of a solute by a fluid invading a network of channels. Examples include aquifer contamination, chemical filtration, and coffee extraction. We propose a continuum equation to model these processes, parametrized by the Péclet number and the rate of mass transfer between the solid and the fluid. We study the time dependence of the extracted mass for different values of the parameter space. The continuum description is validated by combining extraction experiments with coffee and computational fluid dynamics. An analytical solution is derived for the limit of slow mass transfer, which is corroborated by numerical simulations. △ Less

Submitted 19 December, 2023; v1 submitted 6 April, 2023; originally announced April 2023.

Shortcuts to adiabatic population inversion via time-rescaling: stability and thermodynamic cost

Authors: Jair da Silva Andrade, Ângelo F. da Silva França, Bertúlio de Lima Bernardo

Abstract: A shortcut to adiabaticity (STA) is concerned with the fast and robust manipulation of the dynamics of a quantum system that reproduces the effect of an adiabatic process. A recently proposed method enables the generation of shortcuts from a prescribed slow dynamics by simply rescaling the time variable of the quantum evolution operator [B. L. Bernardo, Phys. Rev. Research 2, 013133 (2020)]. This… ▽ More A shortcut to adiabaticity (STA) is concerned with the fast and robust manipulation of the dynamics of a quantum system that reproduces the effect of an adiabatic process. A recently proposed method enables the generation of shortcuts from a prescribed slow dynamics by simply rescaling the time variable of the quantum evolution operator [B. L. Bernardo, Phys. Rev. Research 2, 013133 (2020)]. This time-rescaling method does not demand knowledge about the eigenvalues and eigenstates of the Hamiltonian, and in many cases no additional coupling fields. Here, we use this approach to study the problem of speeding up the population inversion of a two-level quantum system. The fidelity of the dynamics versus systematic errors in the control parameters are shown to be comparable with other STA schemes. From a quantum thermodynamic viewpoint, we also demonstrate that the main properties of the distribution of work required to drive the system along the shortcuts are unchanged with respect to the reference (slow) protocol. △ Less

Submitted 7 July, 2022; v1 submitted 29 April, 2022; originally announced May 2022.

Comments: 10 pages, 4 figures

Journal ref: Sci. Rep. 12, 11538 (2022)

Determination of the recombination coefficient in electrolytic solutions from impedance spectroscopy measurements

Authors: G. Barbero, N. G. Fytas, I. Lelidis, J. V. da Silva Andrade, F. C. M. Freire, A. J. Santana

Abstract: The dependence of the actual/effective ionic concentration on the nominal salt concentration is experimentally investigated for aqueous solutions of NaCl and KCl. The actual ionic bulk density is determined by means of the impedance spectroscopy technique using two approaches for the data analysis: A first approach scrutinizing the real part's plateau of the cell's electric impedance spectrum vers… ▽ More The dependence of the actual/effective ionic concentration on the nominal salt concentration is experimentally investigated for aqueous solutions of NaCl and KCl. The actual ionic bulk density is determined by means of the impedance spectroscopy technique using two approaches for the data analysis: A first approach scrutinizing the real part's plateau of the cell's electric impedance spectrum versus the bulk density of dopant when the impedance is written in the series representation. A second one via extraction of the ionic density from simultaneously fitting of the real and imaginary parts of the impedance spectrum by using the Poisson-Nernst-Planck model. The validity of this procedure reposes on the hypothesis that both the diffusion coefficient and the dielectric constant are independent from the salt concentration in the explored range of concentrations. Our experimental results indicate that a first-order reaction describes well the dissociation phenomenon. The dissociation-association ratio of the reaction has been calculated. Deviations from linearity appear at rather low concentrations of the order of mM. △ Less

Submitted 1 April, 2022; originally announced April 2022.

Comments: 18 pages in preprint style, 12 figures, including appendix

Journal ref: J. Electroanal. Chem. 907, 116070 (2022)

The impact of one-way streets on the asymmetry of the shortest commuting routes

Authors: Hygor P. M. Melo, Diogo P. Mota, José S. Andrade Jr., Nuno A. M. Araújo

Abstract: On a daily commute, the shortest route from home to work rarely overlaps completely the shortest way back. We analyze this asymmetry for several cities and show that it exists even without traffic, due to a non-negligible fraction of one-way streets. For different pairs of origin-destination ($\rm OD$), we compute the log-ratio $r=\ln(\ell_{\rm D}/\ell_{\rm O})$, where $\ell_{\rm O}$ and… ▽ More On a daily commute, the shortest route from home to work rarely overlaps completely the shortest way back. We analyze this asymmetry for several cities and show that it exists even without traffic, due to a non-negligible fraction of one-way streets. For different pairs of origin-destination ($\rm OD$), we compute the log-ratio $r=\ln(\ell_{\rm D}/\ell_{\rm O})$, where $\ell_{\rm O}$ and $\ell_{\rm D}$ are the lengths of the shortest routes from $\rm O$ to $\rm D$ and from $\rm D$ to $\rm O$, respectively. While its average is zero, the amplitude of the fluctuations decays as a power law of the $\rm OD$ shortest path length, $r\sim \ell_{\rm O}^{-β}$. Similarly, the fraction of one-way streets in a shortest route also decays as $\ell_{\rm O}^{-α}$. Based on semi-analytic arguments, we show that $β=(1+α)/2$. Thus, the value of the exponent $β$ is related to correlations in the structure of the underlying street network. △ Less

Submitted 14 November, 2021; originally announced November 2021.

Comments: 6 pages, 4 figures

Mechanical basis and topological routes to cell elimination

Authors: Siavash Monfared, Guruswami Ravichandran, Jose E. Andrade, Amin Doostmohammadi

Abstract: Cell layers eliminate unwanted cells through the extrusion process, which underlines healthy versus flawed tissue behaviors. Although several biochemical pathways have been identified, the underlying mechanical basis including the forces involved in cellular extrusion remain largely unexplored. Utilizing a phase-field model of a three-dimensional cell layer, we study the interplay of cell extrusio… ▽ More Cell layers eliminate unwanted cells through the extrusion process, which underlines healthy versus flawed tissue behaviors. Although several biochemical pathways have been identified, the underlying mechanical basis including the forces involved in cellular extrusion remain largely unexplored. Utilizing a phase-field model of a three-dimensional cell layer, we study the interplay of cell extrusion with cell-cell and cell-substrate interactions, in a flat monolayer. Independent tuning of cell-cell versus cell-substrate adhesion forces reveals that extrusion events can be distinctly linked to defects in nematic and hexatic orders associated with cellular arrangements. Specifically, we show that by increasing relative cell-cell adhesion forces the cell monolayer can switch between the collective tendency towards five-fold, hexatic, disclinations relative to half-integer, nematic, defects for extruding a cell. We unify our findings by accessing three-dimensional mechanical stress fields to show that an extrusion event acts as a mechanism to relieve localized stress concentration. △ Less

Submitted 19 April, 2023; v1 submitted 17 August, 2021; originally announced August 2021.

Tracing contacts to evaluate the transmission of COVID-19 from highly exposed individuals in public transportation

Authors: Caio Ponte, Humberto A. Carmona, Erneson A. Oliveira, Carlos Caminha, Antonio S. Lima Neto, José S. Andrade Jr., Vasco Furtado

Abstract: We investigate, through a data-driven contact tracing model, the transmission of COVID-19 inside buses during distinct phases of the pandemic in a large Brazilian city. From this microscopic approach, we recover the networks of close contacts within consecutive time windows. A longitudinal comparison is then performed by upscaling the traced contacts with the transmission computed from a mean-fiel… ▽ More We investigate, through a data-driven contact tracing model, the transmission of COVID-19 inside buses during distinct phases of the pandemic in a large Brazilian city. From this microscopic approach, we recover the networks of close contacts within consecutive time windows. A longitudinal comparison is then performed by upscaling the traced contacts with the transmission computed from a mean-field compartmental model for the entire city. Our results show that the effective reproduction numbers inside the buses, $Re^{bus}$, and in the city, $Re^{city}$, followed a compatible behavior during the first wave of the local outbreak. Moreover, by distinguishing the close contacts of healthcare workers in the buses, we discovered that their transmission, $Re^{health}$, during the same period, was systematically higher than $Re^{bus}$. This result reinforces the need for special public transportation policies for highly exposed groups of people. △ Less

Submitted 9 June, 2021; originally announced June 2021.

Superspreading k-cores at the center of COVID-19 pandemic persistence

Authors: Matteo Serafino, Higor S. Monteiro, Shaojun Luo, Saulo D. S. Reis, Carles Igual, Antonio S. Lima Neto, Matias Travizano, José S. Andrade, Jr., Hernán A. Makse

Abstract: The spread of COVID-19 caused by the recently discovered SARS-CoV-2 virus has become a worldwide problem with devastating consequences. To slow down the spread of the pandemic, mass quarantines have been implemented globally, provoking further social and economic disruptions. Here, we implement a comprehensive contact tracing network analysis to find an optimized quarantine protocol to dismantle t… ▽ More The spread of COVID-19 caused by the recently discovered SARS-CoV-2 virus has become a worldwide problem with devastating consequences. To slow down the spread of the pandemic, mass quarantines have been implemented globally, provoking further social and economic disruptions. Here, we implement a comprehensive contact tracing network analysis to find an optimized quarantine protocol to dismantle the chain of transmission of coronavirus with minimal disruptions to society. We track billions of anonymized GPS human mobility datapoints from a compilation of hundreds of mobile apps deployed in Latin America to monitor the evolution of the contact network of disease transmission before and after the confinements. As a consequence of the lockdowns, people's mobility across the region decreases by $\sim$53\%, which results in a drastic disintegration of the transmission network by $\sim$90\%. However, this disintegration did not halt the spreading of the disease. Our analysis indicates that superspreading k-core structures persist in the transmission network to prolong the pandemic. Once the k-cores are identified, an optimized strategy to break the chain of transmission is to quarantine a minimal number of 'weak links' with high betweenness centrality connecting the large k-cores. As countries built contact tracing apps to fight the pandemic, our results could turn into a valuable resource to help deploy quarantine protocols with minimized disruptions. △ Less

Submitted 15 March, 2021; originally announced March 2021.

Comments: 7 pages, 4 figures

Journal ref: PLoS Comput Biol 18(4): e1009865 (2022)

Flow through time-evolving porous media: swelling and erosion

Authors: André F. V. Matias, Rodrigo C. V. Coelho, José S. Andrade Jr., Nuno A. M. Araújo

Abstract: The flow through a porous medium strongly depends on the boundary conditions, very often assumed to be static. Here, we consider changes in the medium due to swelling and erosion and extend existing Lattice-Boltzmann models to include both. We study two boundary conditions: a constant pressure drop and a constant flow rate. For a constant flow rate, the steady state depends solely on the erosion d… ▽ More The flow through a porous medium strongly depends on the boundary conditions, very often assumed to be static. Here, we consider changes in the medium due to swelling and erosion and extend existing Lattice-Boltzmann models to include both. We study two boundary conditions: a constant pressure drop and a constant flow rate. For a constant flow rate, the steady state depends solely on the erosion dynamics while for a constant pressure drop it depends also on the timescale of swelling. We analyze the competition between swelling and erosion and identify a transition between regimes where either swelling or erosion dominate. △ Less

Submitted 16 April, 2021; v1 submitted 31 December, 2020; originally announced December 2020.

Comments: 20 pages, 7 figures

Journal ref: Journal of Computational Science (2021) 101360

Frustrated bearings

Authors: R. S. Pires, A. A. Moreira, H. J. Herrmann, J. S. Andrade Jr

Abstract: In a bearing state, touching spheres (disks in two dimensions) roll on each other without slip. Here we frustrate a system of touching spheres by imposing two different bearing states on opposite sides and search for the configurations of lowest energy dissipation. If the dissipation between contacts of spheres is viscous (with random damping constants), the angular momentum continuously changes f… ▽ More In a bearing state, touching spheres (disks in two dimensions) roll on each other without slip. Here we frustrate a system of touching spheres by imposing two different bearing states on opposite sides and search for the configurations of lowest energy dissipation. If the dissipation between contacts of spheres is viscous (with random damping constants), the angular momentum continuously changes from one bearing state to the other. For Coulomb friction (with random friction coefficients) in two dimensions, a sharp line separates the two bearing states and we show that this line corresponds to the minimum cut. Astonishingly however, in three dimensions, intermediate bearing domains, that are not synchronized with either side, are energetically more favorable than the minimum-cut surface. Instead of a sharp cut, the steady state displays a fragmented structure. This novel type of state of minimum dissipation is characterized by a spanning network of slipless contacts that reaches every sphere. Such a situation becomes possible because in three dimensions bearing states have four degrees of freedom. △ Less

Submitted 2 September, 2020; originally announced September 2020.

Comments: 5 pages, 4 figures

Journal ref: Phys. Rev. Lett. 125, 104301 (2020)

Cracking urban mobility

Authors: H. A. Carmona, A. W. T. de Noronha, A. A. Moreira, N. A. M. Araujo, J. S. Andrade Jr

Abstract: Assessing the resilience of a road network is instrumental to improve existing infrastructures and design new ones. Here we apply the optimal path crack model (OPC) to investigate the mobility of road networks and propose a new proxy for resilience of urban mobility. In contrast to static approaches, the OPC accounts for the dynamics of rerouting as a response to traffic jams. Precisely, one simul… ▽ More Assessing the resilience of a road network is instrumental to improve existing infrastructures and design new ones. Here we apply the optimal path crack model (OPC) to investigate the mobility of road networks and propose a new proxy for resilience of urban mobility. In contrast to static approaches, the OPC accounts for the dynamics of rerouting as a response to traffic jams. Precisely, one simulates a sequence of failures (cracks) at the most vulnerable segments of the optimal origin-destination paths that are capable to collapse the system. Our results with synthetic and real road networks reveal that their levels of disorder, fractions of unidirectional segments and spatial correlations can drastically affect the vulnerability to traffic congestion. By applying the OPC to downtown Boston and Manhattan, we found that Boston is significantly more vulnerable than Manhattan. This is compatible with the fact that Boston heads the list of American metropolitan areas with the highest average time waste in traffic. Moreover, our analysis discloses that the origin of this difference comes from the intrinsic spatial correlations of each road network. Finally, we argue that, due to their global influence, the most important cracks identified with OPC can be used to pinpoint potential small rerouting and structural changes in road networks that are capable to substantially improve urban mobility. △ Less

Submitted 31 August, 2020; originally announced August 2020.

Comments: 7 pages, 9 figures

Journal ref: Phys. Rev. Research 2, 043132 (2020)

Circuits with broken fibration symmetries perform core logic computations in biological networks

Authors: Ian Leifer, Flaviano Morone, Saulo D. S. Reis, Jose S. Andrade Jr., Mariano Sigman, Hernan A. Makse

Abstract: We show that logic computational circuits in gene regulatory networks arise from a fibration symmetry breaking in the network structure. From this idea we implement a constructive procedure that reveals a hierarchy of genetic circuits, ubiquitous across species, that are surprising analogues to the emblematic circuits of solid-state electronics: starting from the transistor and progressing to ring… ▽ More We show that logic computational circuits in gene regulatory networks arise from a fibration symmetry breaking in the network structure. From this idea we implement a constructive procedure that reveals a hierarchy of genetic circuits, ubiquitous across species, that are surprising analogues to the emblematic circuits of solid-state electronics: starting from the transistor and progressing to ring oscillators, current-mirror circuits to toggle switches and flip-flops. These canonical variants serve fundamental operations of synchronization and clocks (in their symmetric states) and memory storage (in their broken symmetry states). These conclusions introduce a theoretically principled strategy to search for computational building blocks in biological networks, and present a systematic route to design synthetic biological circuits. △ Less

Submitted 23 June, 2020; originally announced June 2020.

Journal ref: PLoS Comput Biol 2020,16(6): e1007776

Flow through three-dimensional self-affine fractures

Authors: H. J. Seybold, H. A. Carmona, F. A. Leandro Filho, A. D. Araújo, F. Nepomuceno Filho, J. S. Andrade Jr

Abstract: We investigate through numerical simulations of the Navier-Stokes equations the influence of the surface roughness on the fluid flow through fracture joints. Using the Hurst exponent $H$ to characterize the roughness of the self-affine surfaces that constitute the fracture, our analysis reveal the important interplay between geometry and inertia on the flow. Precisely, for low values of Reynolds n… ▽ More We investigate through numerical simulations of the Navier-Stokes equations the influence of the surface roughness on the fluid flow through fracture joints. Using the Hurst exponent $H$ to characterize the roughness of the self-affine surfaces that constitute the fracture, our analysis reveal the important interplay between geometry and inertia on the flow. Precisely, for low values of Reynolds numbers Re, we use Darcy's law to quantify the hydraulic resistance $G$ of the fracture and show that its dependence on $H$ can be explained in terms of a simple geometrical model for the tortuosity $τ$ of the channel. At sufficiently high values of Re, when inertial effects become relevant, our results reveal that nonlinear corrections up to third-order to Darcy's law are aproximately proportional to $H$. These results imply that the resistance $G$ to the flow follows a universal behavior by simply rescaling it in terms of the fracture resistivity and using an effective Reynolds number, namely, Re/$H$. Our results also reveal the presence of quasi-one-dimensional channeling, even considering the absence of shear displacement between upper and lower surfaces of the self-affine fracture. △ Less

Submitted 8 June, 2020; originally announced June 2020.

Comments: 10 pages, 8 figures

Journal ref: Phys. Rev. Fluids 5, 104101 (2020)

Spatio-temporal characteristics of dengue outbreaks

Authors: Saulo D. S. Reis, Lucas Böttcher, João P. da C. Nogueira, Geziel S. Sousa, Antonio S. Lima Neto, Hans J. Herrmann, José S. Andrade Jr

Abstract: After their re-emergence in the last decades, dengue fever and other vector-borne diseases are a potential threat to the lives of millions of people. Based on a data set of dengue cases in the Brazilian city of Fortaleza, collected from 2011 to 2016, we study the spatio-temporal characteristics of dengue outbreaks to characterize epidemic and non-epidemic years. First, we identify regions that sho… ▽ More After their re-emergence in the last decades, dengue fever and other vector-borne diseases are a potential threat to the lives of millions of people. Based on a data set of dengue cases in the Brazilian city of Fortaleza, collected from 2011 to 2016, we study the spatio-temporal characteristics of dengue outbreaks to characterize epidemic and non-epidemic years. First, we identify regions that show a high prevalence of dengue cases and mosquito larvae in different years and also analyze their corresponding correlations. Our results show that the characteristic correlation length of the epidemic is of the order of the system size, suggesting that factors such as citizen mobility may play a major role as a drive for spatial spreading of vector-borne diseases. Inspired by this observation, we perform a mean-field estimation of the basic reproduction number and find that our estimated values agree well with the values reported for other regions, pointing towards similar underlying spreading mechanisms. These findings provide insights into the spreading characteristics of dengue in densely populated areas and should be of relevance for the design of improved disease containment strategies. △ Less

Submitted 4 June, 2020; originally announced June 2020.

Comments: 25 pages, 8 figures

Time evolution of the behaviour of Brazilian legislative Representatives using a complex network approach

Authors: Ludwing Marenco, Humberto A. Carmona, Felipe M. Cardoso, José S. Andrade Jr., Carlos Lenz Cesar

Abstract: The follow up of Representative behavior after elections is imperative for a democratic Representative system, at the very least to punish betrayal with no re-election. Our goal was to show how to follow Representatives' and how to show behavior in real situations and observe trends in political crises including the onset of game changing political instabilities. We used correlation and correlatio… ▽ More The follow up of Representative behavior after elections is imperative for a democratic Representative system, at the very least to punish betrayal with no re-election. Our goal was to show how to follow Representatives' and how to show behavior in real situations and observe trends in political crises including the onset of game changing political instabilities. We used correlation and correlation distance matrices of Brazilian Representative votes during four presidential terms. Re-ordering these matrices with Minimal Spanning Trees displays the dynamical formation of clusters for the sixteen year period, which includes one Presidential impeachment. The reordered matrices, colored by correlation strength and by the parties clearly show the origin of observed clusters and their evolution over time. When large clusters provide government support cluster breaks, political instability arises, which could lead to an impeachment, a trend we observed three years before the Brazilian President was impeached. We believe this method could be applied to foresee other political storms. △ Less

Submitted 6 February, 2020; v1 submitted 3 May, 2019; originally announced May 2019.

Comments: 11 pages, 4 Figures

Report number: 2020

Journal ref: PLoS ONE 15(2): e0226504 (2020)

Dynamics of racial segregation and gentrification in New York City

Authors: Felipe G. Operti, André A. Moreira, Andrea Gabrielli, Hernan Makse, José S. Andrade Jr

Abstract: Racial residential segregation is interconnected with several other phenomena such as income inequalities, property values inequalities, and racial disparities in health and in education. Furthermore, recent literature suggests the phenomena of gentrification as a cause of perpetuation or increase of racial residential segregation in some American cities. In this paper, we analyze the dynamics of… ▽ More Racial residential segregation is interconnected with several other phenomena such as income inequalities, property values inequalities, and racial disparities in health and in education. Furthermore, recent literature suggests the phenomena of gentrification as a cause of perpetuation or increase of racial residential segregation in some American cities. In this paper, we analyze the dynamics of racial residential segregation for white, black, Asian, and Hispanic citizens in New York City in the years of 1990, 2000, and 2010. It was possible to observe that segregation between white and Hispanic citizens, and discrimination between white and Asian ones has grown, while segregation between white and black is quite stable. Furthermore, we analyzed the per capita income and the Gini coefficient in each segregated zone, showing that the highest inequalities occur in the zones where there is overlap of high-density zones of pair of races. Focusing on census tracts that have changed density of population during these twenty years, and, particularly, by analyzing white and black people's segregation, our analysis reveals that a positive flux of white (black) people is associated to a substantial increase (decrease) of the property values, as compared with the city mean. Furthermore, by clustering the region of high density of black citizens, we measured the variation of area and displacement of the four biggest clusters in the period from 1990 to 2010. The large displacements (~ 1.6 Km) observed for two of these clusters, namely, one in the neighborhood of Harlem and the other inside the borough of Brooklyn, led to the emergence of typically gentrified regions. △ Less

Submitted 17 April, 2019; v1 submitted 15 April, 2019; originally announced April 2019.

Comments: 19 pages, 12 figures

A universal approach for drainage basins

Authors: Erneson A. Oliveira, Rilder S. Pires, Rubens S. Oliveira, Vasco Furtado, Hans J. Herrmann, José S. Andrade Jr

Abstract: Drainage basins are essential to Geohydrology and Biodiversity. Defining those regions in a simple, robust and efficient way is a constant challenge in Earth Science. Here, we introduce a model to delineate multiple drainage basins through an extension of the Invasion Percolation-Based Algorithm (IPBA). In order to prove the potential of our approach, we apply it to real and artificial datasets. W… ▽ More Drainage basins are essential to Geohydrology and Biodiversity. Defining those regions in a simple, robust and efficient way is a constant challenge in Earth Science. Here, we introduce a model to delineate multiple drainage basins through an extension of the Invasion Percolation-Based Algorithm (IPBA). In order to prove the potential of our approach, we apply it to real and artificial datasets. We observe that the perimeter and area distributions of basins and anti-basins display long tails extending over several orders of magnitude and following approximately power-law behaviors. Moreover, the exponents of these power laws depend on spatial correlations and are invariant under the landscape orientation, not only for terrestrial, but lunar and martian landscapes. The terrestrial and martian results are statistically identical, which suggests that a hypothetical martian river would present similarity to the terrestrial rivers. Finally, we propose a theoretical value for the Hack's exponent based on the fractal dimension of watersheds, $γ=D/2$. We measure $γ=0.54 \pm 0.01$ for Earth, which is close to our estimation of $γ\approx 0.55$. Our study suggests that Hack's law can have its origin purely in the maximum and minimum lines of the landscapes. △ Less

Submitted 8 July, 2019; v1 submitted 20 December, 2018; originally announced December 2018.

Comments: 20 pages, 6 Figures, and 1 Table

Fundraising and vote distribution: a non-equilibrium statistical approach

Authors: H. P. M. Melo, N. A. M. Araujo, J. S. Andrade Jr

Abstract: The number of votes correlates strongly with the money spent in a campaign, but the relation between the two is not straightforward. Among other factors, the output of a ballot depends on the number of candidates, voters, and available resources. Here, we develop a conceptual framework based on Shannon entropy maximization and Superstatistics to establish a relation between the distributions of mo… ▽ More The number of votes correlates strongly with the money spent in a campaign, but the relation between the two is not straightforward. Among other factors, the output of a ballot depends on the number of candidates, voters, and available resources. Here, we develop a conceptual framework based on Shannon entropy maximization and Superstatistics to establish a relation between the distributions of money spent by candidates and their votes. By establishing such a relation, we provide a tool to predict the outcome of a ballot and to alert for possible misconduct either in the report of fundraising and spending of campaigns or on vote counting. As an example, we consider real data from a proportional election with $6323$ candidates, where a detailed data verification is virtually impossible, and show that the number of potential misconducting candidates to audit can be reduced to only nine. △ Less

Submitted 3 October, 2018; originally announced October 2018.

The hidden traits of endemic illiteracy in cities

Authors: Luiz G. A. Alves, Jose S. Andrade Jr., Quentin S. Hanley, Haroldo V. Ribeiro

Abstract: In spite of the considerable progress towards reducing illiteracy rates, many countries, including developed ones, have encountered difficulty achieving further reduction in these rates. This is worrying because illiteracy has been related to numerous health, social, and economic problems. Here, we show that the spatial patterns of illiteracy in urban systems have several features analogous to the… ▽ More In spite of the considerable progress towards reducing illiteracy rates, many countries, including developed ones, have encountered difficulty achieving further reduction in these rates. This is worrying because illiteracy has been related to numerous health, social, and economic problems. Here, we show that the spatial patterns of illiteracy in urban systems have several features analogous to the spread of diseases such as dengue and obesity. Our results reveal that illiteracy rates are spatially long-range correlated, displaying non-trivial clustering structures characterized by percolation-like transitions and fractality. These patterns can be described in the context of percolation theory of long-range correlated systems at criticality. Together, these results provide evidence that the illiteracy incidence can be related to a transmissible process, in which the lack of access to minimal education propagates in a population in a similar fashion to endemic diseases. △ Less

Submitted 26 September, 2018; originally announced September 2018.

Comments: 19 pages, 4 figures, 1 table, supplementary information; accepted for publication in Physica A

Journal ref: Physica A 515, 566-574 (2019)

CFD simulation of the wind field over a terrain with sand fences: Critical spacing for the protected soil area

Authors: Izael A. Lima, Eric J. R. Parteli, Yaping Shao, José S. Andrade Jr., Hans J. Herrmann, Ascânio D. Araújo

Abstract: Sand fences are often erected to reduce wind speed, prevent aeolian soil erosion, and induce sand deposition and dune formation in areas affected by sand encroachment and desertification. However, the search for the most efficient array of fences by means of field experiments alone poses a challenging task given that field experiments are affected by weather conditions. Here we apply Computational… ▽ More Sand fences are often erected to reduce wind speed, prevent aeolian soil erosion, and induce sand deposition and dune formation in areas affected by sand encroachment and desertification. However, the search for the most efficient array of fences by means of field experiments alone poses a challenging task given that field experiments are affected by weather conditions. Here we apply Computational Fluid Dynamic simulations to investigate the three-dimensional structure of the turbulent wind field over an array of fences of different sizes, porosity and spacing. Our goal is to perform a quantitative analysis of this structure in the absence of saltation or suspension transport, as first step toward the development of a continuum simulation of the aeolian soil in presence of the fences. We find that the area of soil protected against direct aerodynamic entrainment has two regimes, depending on the spacing $L_x$ between the fences. When $L_x$ is smaller than a critical value $L_{xc}$, the wake zones associated with each fence are inter-connected (regime A), while these wake zones appear separated from each other (regime B) when $L_x$ exceeds this critical value of spacing. The system undergoes a second order phase transition at $L_x = L_{xc}$, with the cross-wind width of the protected zone scaling with $\left[{1-{L}_x}/{L}_{xc}\right]^β$ in regime A, with $β\approx 0.32$. Our findings have implication for a better understanding of aeolian transport in the presence of sand fences, as well as to develop optimization strategies for measures to protect soils from wind erosion. △ Less

Submitted 1 April, 2019; v1 submitted 27 August, 2018; originally announced August 2018.

Comments: 35 pages, 13 figures

Density waves and jet emission asymmetry in Bose Fireworks

Authors: Han Fu, Lei Feng, Brandon M. Anderson, Logan W. Clark, Jiazhong Hu, Jeffery W. Andrade, Cheng Chin, K. Levin

Abstract: A Bose condensate subject to a periodic modulation of the two-body interactions was recently observed to emit matter-wave jets resembling "fireworks" [Nature 551, 356(2017)]. In this paper, combining experiment with numerical simulation, we demonstrate that these "Bose fireworks" represent a late stage in a complex time evolution of the driven condensate. We identify a "density wave" stage which p… ▽ More A Bose condensate subject to a periodic modulation of the two-body interactions was recently observed to emit matter-wave jets resembling "fireworks" [Nature 551, 356(2017)]. In this paper, combining experiment with numerical simulation, we demonstrate that these "Bose fireworks" represent a late stage in a complex time evolution of the driven condensate. We identify a "density wave" stage which precedes jet emission and results from interference of matterwaves. The density waves self-organize and self-amplify without the breaking of long range translational symmetry. Importantly, this density wave structure deterministically establishes the template for the subsequent patterns of the emitted jets. Our simulations, in good agreement with experiment, also address the apparent asymmetry in the jet pattern and show it is fully consistent with momentum conservation. △ Less

Submitted 23 July, 2018; originally announced July 2018.

Journal ref: Phys. Rev. Lett. 121, 243001 (2018)

Self-organization and Nonuniversal Anomalous Scaling in Non-Newtonian Turbulence

Authors: H. J. Seybold, H. A. Carmona, H. J. Herrmann, J. S. Andrade

Abstract: We investigate through Direct Numerical Simulations (DNS) the statistical properties of turbulent flows in the inertial subrange for non-Newtonian power-law fluids. The structural invariance found for the vortex size distribution is achieved through a self-organized mechanism at the microscopic scale of the turbulent motion that adjusts, according to the rheological properties of the fluid, the ra… ▽ More We investigate through Direct Numerical Simulations (DNS) the statistical properties of turbulent flows in the inertial subrange for non-Newtonian power-law fluids. The structural invariance found for the vortex size distribution is achieved through a self-organized mechanism at the microscopic scale of the turbulent motion that adjusts, according to the rheological properties of the fluid, the ratio between the viscous dissipations inside and outside the vortices. Moreover, the deviations from the K41 theory of the structure functions' exponents reveal that the anomalous scaling exhibits a systematic nonuniversal behavior with respect to the rheological properties of the fluids. △ Less

Submitted 20 June, 2018; originally announced June 2018.

Comments: 5 pages, 4 figures

Journal ref: Phys. Rev. Fluids 4, 064604 (2019)

On the decentralized navigation of multiple packages on transportation networks

Authors: Samuel M. da Silva, Saulo D. S. Reis, Ascânio D. Araújo, José S. Andrade, Jr

Abstract: We investigate by numerical simulation and finite-size analysis the impact of long-range shortcuts on a spatially embedded transportation network. Our networks are built from two-dimensional ($d=2$) square lattices to be improved by the addition of long-range shortcuts added with probability $P(r_{ij})\sim r_{ij}^{-α}$ [J. M. Kleinberg, Nature 406, 845 (2000)]. Considering those improved networks,… ▽ More We investigate by numerical simulation and finite-size analysis the impact of long-range shortcuts on a spatially embedded transportation network. Our networks are built from two-dimensional ($d=2$) square lattices to be improved by the addition of long-range shortcuts added with probability $P(r_{ij})\sim r_{ij}^{-α}$ [J. M. Kleinberg, Nature 406, 845 (2000)]. Considering those improved networks, we performed numerical simulation of multiple discrete package navigation and found a limit for the amount of packages flowing through the network. Such limit is characterized by a critical probability of creating packages $p_{c}$, where above this value a transition to a congested state occurs. Moreover, $p_{c}$ is found to follow a power-law, $p_{c}\sim L^{-γ}$, where $L$ is the network size. Our results indicate the presence of an optimal value of $α_{\rm min}\approx1.7$, where the parameter $γ$ reaches its minimum value and the networks are more resilient to congestion for larger system sizes. △ Less

Submitted 15 June, 2018; originally announced June 2018.

Comments: 7 pages, 5 figures

Journal ref: Phys. Rev. E 98, 032306 (2018)

All-optical attoclock for imaging tunnelling wavepackets

Authors: I. Babushkin, A. J. Galan, J. R. C. Andrade, A. Husakou, F. Morales, M. Kretschmar, T. Nagy, V. Vaičaitis, L. Shi, D. Zuber, L. Bergé, S. Skupin, I. A. Nikolaeva, N. A. Panov, D. E. Shipilo, O. G. Kosareva, A. N. Pfeiffer, A. Demircan, M. J. J. Vrakking, U. Morgner, M. Ivanov

Abstract: Recent experiments on measuring time-delays during tunnelling of cold atoms through an optically created potential barrier are reinvigorating the controversial debate regarding possible time-delays during light-induced tunnelling of an electron from an atom. Compelling theoretical and experimental arguments have been put forward to advocate opposite views, confirming or refuting the existence of f… ▽ More Recent experiments on measuring time-delays during tunnelling of cold atoms through an optically created potential barrier are reinvigorating the controversial debate regarding possible time-delays during light-induced tunnelling of an electron from an atom. Compelling theoretical and experimental arguments have been put forward to advocate opposite views, confirming or refuting the existence of finite tunnelling time delays. Yet, such a delay, whether present or not, is but a single quantity characterizing the tunnelling wavepacket; the underlying dynamics are richer. Here we propose to augment photo-electron detection in laser-induced tunnelling with detection of light emitted by the tunnelling electron -- the so-called Brunel radiation. Using a combination of single-color and two-color driving fields, we identify the all-optical signatures of the re-shaping of the tunnelling wavepacket as it emerges from the tunnelling barrier and moves away from the core. This reshaping includes not only an effective time-delay but also time-reversal asymmetry of the ionization process, which we describe theoretically and observe experimentally. We show how both delay and reshaping are mapped on the polarization properties of the Brunel radiation, with different harmonics behaving as different hands of a clock moving at different speeds. The all-optical detection paves the way to time-resolving optical tunnelling in condensed matter systems, e.g. tunnelling across bandgaps in solids, on the attosecond time-scale. △ Less

Submitted 20 September, 2023; v1 submitted 12 March, 2018; originally announced March 2018.

Comments: In the last version, only an acknowledgement has been changed in comparison to the pre-last one

A worldwide model for boundaries of urban settlements

Authors: Erneson A. Oliveira, Vasco Furtado, José S. Andrade Jr., Hernán A. Makse

Abstract: The shape of urban settlements plays a fundamental role in their sustainable planning. Properly defining the boundaries of cities is challenging and remains an open problem in the Science of Cities. Here, we propose a worldwide model to define urban settlements beyond their administrative boundaries through a bottom-up approach that takes into account geographical biases intrinsically associated w… ▽ More The shape of urban settlements plays a fundamental role in their sustainable planning. Properly defining the boundaries of cities is challenging and remains an open problem in the Science of Cities. Here, we propose a worldwide model to define urban settlements beyond their administrative boundaries through a bottom-up approach that takes into account geographical biases intrinsically associated with most societies around the world, and reflected in their different regional growing dynamics. The generality of the model allows to study the scaling laws of cities at all geographical levels: countries, continents, and the entire world. Our definition of cities is robust and holds to one of the most famous results in Social Sciences: Zipf's law. According to our results, the largest cities in the world are not in line with what was recently reported by the United Nations. For example, we find that the largest city in the world is an agglomeration of several small settlements close to each other, connecting three large settlements: Alexandria, Cairo, and Luxor. Our definition of cities opens the doors to the study of the economy of cities in a systematic way independently of arbitrary definitions that employ administrative boundaries. △ Less

Submitted 23 May, 2018; v1 submitted 16 November, 2017; originally announced November 2017.

The price of a vote: diseconomy in proportional elections

Authors: Hygor Piaget M. Melo, Saulo D. S. Reis, André A. Moreira, Hernán A. Makse, José S. Andrade Jr

Abstract: The increasing cost of electoral campaigns raises the need for effective campaign planning and a precise understanding of the return of such investment. Interestingly, despite the strong impact of elections on our daily lives, how this investment is translated into votes is still unknown. By performing data analysis and modeling, we show that top candidates spend more money \emph{per vote} than th… ▽ More The increasing cost of electoral campaigns raises the need for effective campaign planning and a precise understanding of the return of such investment. Interestingly, despite the strong impact of elections on our daily lives, how this investment is translated into votes is still unknown. By performing data analysis and modeling, we show that top candidates spend more money \emph{per vote} than the less successful and poorer candidates, a sublinearity that discloses a diseconomy of scale. We demonstrate that such electoral diseconomy arises from the competition between candidates due to inefficient campaign expenditure. Our approach succeeds in two important tests. First, it reveals that the statistical pattern in the vote distribution of candidates can be explained in terms of the independently conceived, but similarly skewed distribution of money campaign. Second, using a heuristic argument, we are able to predict a turnout percentage for a given election of approximately 63\%. This result is in good agreement with the average turnout rate obtained from real data. Due to its generality, we expect that our approach can be applied to a wide range of problems concerning the adoption process in marketing campaigns. △ Less

Submitted 15 November, 2017; originally announced November 2017.

Comments: 8 pages, 4 figures

Targeted Recovery as an Effective Strategy against Epidemic Spreading

Authors: Lucas Böttcher, José S. Andrade Jr., Hans J. Herrmann

Abstract: We propose a targeted intervention protocol where recovery is restricted to individuals that have the least number of infected neighbours. Our recovery strategy is highly efficient on any kind of network, since epidemic outbreaks are minimal when compared to the baseline scenario of spontaneous recovery. In the case of spatially embedded networks, we find that an epidemic stays strongly spatially… ▽ More We propose a targeted intervention protocol where recovery is restricted to individuals that have the least number of infected neighbours. Our recovery strategy is highly efficient on any kind of network, since epidemic outbreaks are minimal when compared to the baseline scenario of spontaneous recovery. In the case of spatially embedded networks, we find that an epidemic stays strongly spatially confined with a characteristic length scale undergoing a random walk. We demonstrate numerically and analytically that this dynamics leads to an epidemic spot with a flat surface structure and a radius that grows linearly with the spreading rate. △ Less

Submitted 16 October, 2017; originally announced October 2017.

Comments: 6 pages, 5 figures

Journal ref: Sci. Rep. 7, 14356 (2017)

Polymer/2D material nanocomposite manufacturing beyond laboratory frontiers

Authors: Pablo A. R. Munoz, Camila F. P. de Oliveira, Leice G. Amurin, Camila L. C. Rodriguez, Danilo A. Nagaoka, Maria Inês Bruno Tavares, Sergio H. Domingues, Ricardo J. E. Andrade, Guilhermino J. M. Fechine

Abstract: Polymer nanocomposites based on 2D materials as fillers are the target in the industrial sector, but the ability to manufacture them on a large scale is very limited, and there is a lack of tools to scale up the manufacturing process of these nanocomposites. Here, for the first time, a systematic and fundamental study showing how 2D materials are inserted into the polymeric matrix in order to obta… ▽ More Polymer nanocomposites based on 2D materials as fillers are the target in the industrial sector, but the ability to manufacture them on a large scale is very limited, and there is a lack of tools to scale up the manufacturing process of these nanocomposites. Here, for the first time, a systematic and fundamental study showing how 2D materials are inserted into the polymeric matrix in order to obtain nanocomposites using conventional and industrially scalable polymer processing machines leading to large-scale manufacturing are described. Two new strategies were used to insert pre-exfoliated 2D material into the polymer matrix, liquid-phase feeder, and solid-solid deposition. Characterizations were beyond micro and nanoscale, allowing the evaluation of the morphology for millimeter samples size. The methodologies described here are extendable to all thermoplastic polymers and 2D materials providing nanocomposites with suitable morphology to obtain singular properties and also triggering the start of the manufacturing process on a large scale. △ Less

Submitted 5 October, 2017; originally announced October 2017.

Comments: 19 pages, 13 figures, 7 tables

The light pollution as a surrogate for urban population of the US cities

Authors: Felipe G. Operti, Erneson A. Oliveira, Humberto A. Carmona, Javam C. Machado, José S. Andrade Jr

Abstract: We show that the definition of the city boundaries can have a dramatic influence on the scaling behavior of the night-time light (NTL) as a function of population (POP) in the US. Precisely, our results show that the arbitrary geopolitical definition based on the Metropolitan/Consolidated Metropolitan Statistical Areas (MSA/CMSA) leads to a sublinear power-law growth of NTL with POP. On the other… ▽ More We show that the definition of the city boundaries can have a dramatic influence on the scaling behavior of the night-time light (NTL) as a function of population (POP) in the US. Precisely, our results show that the arbitrary geopolitical definition based on the Metropolitan/Consolidated Metropolitan Statistical Areas (MSA/CMSA) leads to a sublinear power-law growth of NTL with POP. On the other hand, when cities are defined according to a more natural agglomeration criteria, namely, the City Clustering Algorithm (CCA), an isometric relation emerges between NTL and population. This discrepancy is compatible with results from previous works showing that the scaling behaviors of various urban indicators with population can be substantially different for distinct definitions of city boundaries. Moreover, considering the CCA definition as more adequate than the MSA/CMSA one because the former does not violate the expected extensivity between land population and area of their generated clusters, we conclude that, without loss of generality, the CCA measures of light pollution and population could be interchangeably utilized in future studies. △ Less

Submitted 15 June, 2017; originally announced June 2017.

Optimal array of sand fences

Authors: Izael A. Lima, Ascânio D. Araújo, Eric J. R. Parteli, José S. Andrade Jr., Hans J. Herrmann

Abstract: Sand fences are widely applied to prevent soil erosion by wind in areas affected by desertification. Sand fences also provide a way to reduce the emission rate of dust particles, which is triggered mainly by the impacts of wind-blown sand grains onto the soil and affects the Earth's climate. Many different types of fence have been designed and their effects on the sediment transport dynamics studi… ▽ More Sand fences are widely applied to prevent soil erosion by wind in areas affected by desertification. Sand fences also provide a way to reduce the emission rate of dust particles, which is triggered mainly by the impacts of wind-blown sand grains onto the soil and affects the Earth's climate. Many different types of fence have been designed and their effects on the sediment transport dynamics studied since many years. However, the search for the optimal array of fences has remained largely an empirical task. In order to achieve maximal soil protection using the minimal amount of fence material, a quantitative understanding of the flow profile over the relief encompassing the area to be protected including all employed fences is required. Here we use Computational Fluid Dynamics to calculate the average turbulent airflow through an array of fences as a function of the porosity, spacing and height of the fences. Specifically, we investigate the factors controlling the fraction of soil area over which the basal average wind shear velocity drops below the threshold for sand transport when the fences are applied. We introduce a cost function, given by the amount of material necessary to construct the fences. We find that, for typical sand-moving wind velocities, the optimal fence height (which minimizes this cost function) is around $50\,$cm, while using fences of height around $1.25\,$m leads to maximal cost. △ Less

Submitted 16 February, 2017; originally announced February 2017.

Comments: 10 pages

Surface Plasmon Assisted Gentle Ablation of Nanostructures by Femtosecond Oscillator

Authors: Liping Shi, Bianca Iwan, Quentin Ripault, Jose R. C. Andrade, Seunghwoi Han, Hyunwoong Kim, Willem Boutu, Dominik Franz, Rana Nicolas, Torsten Heidenblut, Carsten Reinhardt, Bert Bastiaens, Tamas Nagy, Ihar Babuskin, Uwe Morgner, Seung-Woo Kim, Gunter Steinmeyer, Hamed Merdji, Milutin Kovacev

Abstract: We experimentally demonstrate the use of subwavelength optical nanoantennae to assist the gentle ablation of nanostructures directly using ultralow fluence from a Ti: sapphire oscillator through the excitation of surface plasmon waves. We show that this ablation mechanism is the same for metal and dielectric. The analytical solutions of ablation threshold are in excellent agreement with the experi… ▽ More We experimentally demonstrate the use of subwavelength optical nanoantennae to assist the gentle ablation of nanostructures directly using ultralow fluence from a Ti: sapphire oscillator through the excitation of surface plasmon waves. We show that this ablation mechanism is the same for metal and dielectric. The analytical solutions of ablation threshold are in excellent agreement with the experiment estimations. Surface plasmon assisted locally enhanced ablation at nanoscale provides a method for nanomachining, manipulation and modification the nanostructures without collateral thermal damage to the materials. It is also shown that this ablation can deposit low-density high quality thin nano film. △ Less

Submitted 13 January, 2017; originally announced January 2017.

Comments: 4 pages, 4 figures

Human Mobility in Large Cities as a Proxy for Crime

Authors: Carlos Caminha, Vasco Furtado, Tarcisio H. C. Pequeno, Caio Ponte, Hygor P. M. Melo, Erneson A. Oliveira, José S. Andrade Jr

Abstract: We investigate at the subscale of the neighborhoods of a highly populated city the incidence of property crimes in terms of both the resident and the floating population. Our results show that a relevant allometric relation could only be observed between property crimes and floating population. More precisely, the evidence of a superlinear behavior indicates that a disproportional number of proper… ▽ More We investigate at the subscale of the neighborhoods of a highly populated city the incidence of property crimes in terms of both the resident and the floating population. Our results show that a relevant allometric relation could only be observed between property crimes and floating population. More precisely, the evidence of a superlinear behavior indicates that a disproportional number of property crimes occurs in regions where an increased flow of people takes place in the city. For comparison, we also found that the number of crimes of peace disturbance only correlates well, and in a superlinear fashion too, with the resident population. Our study raises the interesting possibility that the superlinearity observed in previous studies [Bettencourt et al., Proc. Natl. Acad. Sci. USA 104, 7301 (2007) and Melo et al., Sci. Rep. 4, 6239 (2014)] for homicides versus population at the city scale could have its origin in the fact that the floating population, and not the resident one, should be taken as the relevant variable determining the intrinsic microdynamical behavior of the system. △ Less

Submitted 27 October, 2016; originally announced October 2016.

Comments: 17 pages, 8 Figures

Crumpling Damaged Graphene

Authors: I. Giordanelli, M. Mendoza, J. S. Andrade, Jr., M. A. F. Gomes, H. J. Herrmann

Abstract: Through molecular mechanics we find that non-covalent interactions modify the fractality of crumpled damaged graphene. Pristine graphene membranes are damaged by adding random vacancies and carbon-hydrogen bonds. Crumpled membranes exhibit a fractal dimension of $ 2.71 \pm 0.02$ when all interactions between carbon atoms are considered, and $2.30 \pm 0.05$ when non-covalent interactions are suppre… ▽ More Through molecular mechanics we find that non-covalent interactions modify the fractality of crumpled damaged graphene. Pristine graphene membranes are damaged by adding random vacancies and carbon-hydrogen bonds. Crumpled membranes exhibit a fractal dimension of $ 2.71 \pm 0.02$ when all interactions between carbon atoms are considered, and $2.30 \pm 0.05$ when non-covalent interactions are suppressed. The transition between these two values, obtained by switching on/off the non-covalent interactions of equilibrium configurations, is shown to be reversible and independent on thermalisation. In order to explain this transition, we propose a theoretical model that is compatible with our numerical findings. Finally, we also compare damaged graphene membranes with other crumpled structures, as for instance, polymerised membranes and paper sheets, that share similar scaling properties. △ Less

Submitted 23 May, 2016; originally announced May 2016.

Journal ref: Scientific Reports 6 (2016): 25891

Model for the growth of the World Airline Network

Authors: T. Verma, N. A. M. Araújo, J. Nagler, J. S. Andrade Jr., H. J. Herrmann

Abstract: We propose a probabilistic growth model for transport networks which employs a balance between popularity of nodes and the physical distance between nodes. By comparing the degree of each node in the model network and the WAN, we observe that the difference between the two is minimized for $α\approx 2$. Interestingly, this is the value obtained for the node-node correlation function in the WAN. Th… ▽ More We propose a probabilistic growth model for transport networks which employs a balance between popularity of nodes and the physical distance between nodes. By comparing the degree of each node in the model network and the WAN, we observe that the difference between the two is minimized for $α\approx 2$. Interestingly, this is the value obtained for the node-node correlation function in the WAN. This suggests that our model explains quite well the growth of airline networks. △ Less

Submitted 30 January, 2016; originally announced February 2016.

Comments: 4 pages, 3 figures

In Vitro Durability - Pivot bearing with Diamond Like Carbon for Ventricular Assist Devices

Authors: Rosa Corrêa Leoncio de Sá, Vladimir Jesus Trava Airoldi, Tarcísio Fernandes Leão, Evandro Drigo da Silva, Jeison Willian Gomes da Fonseca, Bruno Utiyama da Silva, Edir Branzoni Leal, João Roberto Moro, Aron José Pazin de Andrade, Eduardo Guy Perpétuo Bock

Abstract: Institute Dante Pazzanese of Cardiology (IDPC) develops Ventricular Assist Devices (VAD) that can stabilize the hemodynamics of patients with severe heart failure before, during and/or after the medical practice; can be temporary or permanent. The ADV's centrifugal basically consist of a rotor suspended for system pivoting bearing; the PIVOT is the axis with movement of rotational and the bearing… ▽ More Institute Dante Pazzanese of Cardiology (IDPC) develops Ventricular Assist Devices (VAD) that can stabilize the hemodynamics of patients with severe heart failure before, during and/or after the medical practice; can be temporary or permanent. The ADV's centrifugal basically consist of a rotor suspended for system pivoting bearing; the PIVOT is the axis with movement of rotational and the bearing is the bearing surface. As a whole system of an implantable VAD should be made of long-life biomaterial so that there is no degradation or deformation during application time; surface modification techniques have been widely studied and implemented to improve properties such as biocompatibility and durability of applicable materials. The Chemical Vapour Deposition technique allows substrates having melting point higher than 300 °C to be coated, encapsulated, with a diamond like carbon film (DLC); The test simulated the actual conditions in which the system of support remains while applying a ADV. The results have been prepared on a comparative basis; where you select the pivoting assembly which showed less deformation by abrasive wear. △ Less

Submitted 31 October, 2015; originally announced November 2015.

Comments: sent by Eduardo Bock

MSC Class: 00A79 ACM Class: B.2.2; E.0

Entropy Production and the Pressure-Volume Curve of the Lung

Authors: Cláudio L. N. Oliveira, Ascânio D. Araújo, Jason H. T. Bates, José S. Andrade Jr., Béla Suki

Abstract: We investigate analytically the production of entropy during a breathing cycle in healthy and diseased lungs. First, we calculate entropy production in healthy lungs by applying the laws of thermodynamics to the well-known transpulmonary pressure-volume ($P-V$) curves of the lung under the assumption that lung tissue behaves as an entropy spring-like rubber. The bulk modulus, $B$, of the lung is a… ▽ More We investigate analytically the production of entropy during a breathing cycle in healthy and diseased lungs. First, we calculate entropy production in healthy lungs by applying the laws of thermodynamics to the well-known transpulmonary pressure-volume ($P-V$) curves of the lung under the assumption that lung tissue behaves as an entropy spring-like rubber. The bulk modulus, $B$, of the lung is also derived from these calculations. Second, we extend this approach to elastic recoil disorders of the lung such as occur in pulmonary fibrosis and emphysema. These diseases are characterized by particular alterations in the $P-V$ relationship. For example, in fibrotic lungs $B$ increases monotonically with disease progression, while in emphysema the opposite occurs. These diseases can thus be mimicked simply by making appropriate adjustments to the parameters of the $P-V$ curve. Using Clausius's formalism, we show that entropy production, $ΔS$, is related to the hysteresis area, $ΔA$, enclosed by the $P-V$ curve during a breathing cycle, namely, $ΔS=ΔA/T$, where $T$ is the body temperature. Although $ΔA$ is highly dependent on the disease, such formula applies to healthy as well as diseased lungs, regardless of the disease stage. Finally, we use {\it ansatzs} to predict analytically the entropy produced by the fibrotic and emphysematous lungs. △ Less

Submitted 28 September, 2015; originally announced September 2015.

The cuttlefish Sepia officinalis (Sepiidae, Cephalopoda) constructs cuttlebone from a liquid-crystal precursor

Authors: Antonio G. Checa, Julyan H. E. Cartwright, Isabel Sánchez-Almazo, José P. Andrade, Francisco Ruiz-Raya

Abstract: Cuttlebone, the sophisticated buoyancy device of cuttlefish, is made of extensive superposed chambers that have a complex internal arrangement of calcified pillars and organic membranes. It has not been clear how this structure is assembled. We find that the membranes result from a myriad of minor membranes initially filling the whole chamber, made of nanofibres evenly oriented within each membran… ▽ More Cuttlebone, the sophisticated buoyancy device of cuttlefish, is made of extensive superposed chambers that have a complex internal arrangement of calcified pillars and organic membranes. It has not been clear how this structure is assembled. We find that the membranes result from a myriad of minor membranes initially filling the whole chamber, made of nanofibres evenly oriented within each membrane and slightly rotated with respect to those of adjacent membranes, producing a helical arrangement. We propose that the organism secretes a chitin-protein complex, which self-organizes layer-by-layer as a cholesteric liquid crystal, whereas the pillars are made by viscous fingering. The liquid crystallization mechanism permits us to homologize the elements of the cuttlebone with those of other coleoids and with the nacreous septa and the shells of nautiloids. These results challenge our view of this ultra-light natural material possessing desirable mechanical, structural and biological properties, suggesting that two self-organizing physical principles suffice to understand its formation. △ Less

Submitted 27 June, 2015; originally announced June 2015.

Journal ref: Scientific Reports 5, 11513, 2015

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.

Comments: 5 pages, 4 figures

Watersheds in disordered media

Authors: N. A. M. Araújo, K. J. Schrenk, H. J. Herrmann, J. S. Andrade Jr

Abstract: What is the best way to divide a rugged landscape? Since ancient times, watersheds separating adjacent water systems that flow, for example, toward different seas, have been used to delimit boundaries. Interestingly, serious and even tense border disputes between countries have relied on the subtle geometrical properties of these tortuous lines. For instance, slight and even anthropogenic modifica… ▽ More What is the best way to divide a rugged landscape? Since ancient times, watersheds separating adjacent water systems that flow, for example, toward different seas, have been used to delimit boundaries. Interestingly, serious and even tense border disputes between countries have relied on the subtle geometrical properties of these tortuous lines. For instance, slight and even anthropogenic modifications of landscapes can produce large changes in a watershed, and the effects can be highly nonlocal. Although the watershed concept arises naturally in geomorphology, where it plays a fundamental role in water management, landslide, and flood prevention, it also has important applications in seemingly unrelated fields such as image processing and medicine. Despite the far-reaching consequences of the scaling properties on watershed-related hydrological and political issues, it was only recently that a more profound and revealing connection has been disclosed between the concept of watershed and statistical physics of disordered systems. This review initially surveys the origin and definition of a watershed line in a geomorphological framework to subsequently introduce its basic geometrical and physical properties. Results on statistical properties of watersheds obtained from artificial model landscapes generated with long-range correlations are presented and shown to be in good qualitative and quantitative agreement with real landscapes. △ Less

Submitted 18 December, 2014; originally announced December 2014.

How does public opinion become extreme?

Authors: Marlon Ramos, Jia Shao, Saulo D. S. Reis, Celia Anteneodo, José S. Andrade Jr, Shlomo Havlin, Hernán A. Makse

Abstract: We investigate the emergence of extreme opinion trends in society by employing statistical physics modeling and analysis on polls that inquire about a wide range of issues such as religion, economics, politics, abortion, extramarital sex, books, movies, and electoral vote. The surveys lay out a clear indicator of the rise of extreme views. The precursor is a nonlinear relation between the fraction… ▽ More We investigate the emergence of extreme opinion trends in society by employing statistical physics modeling and analysis on polls that inquire about a wide range of issues such as religion, economics, politics, abortion, extramarital sex, books, movies, and electoral vote. The surveys lay out a clear indicator of the rise of extreme views. The precursor is a nonlinear relation between the fraction of individuals holding a certain extreme view and the fraction of individuals that includes also moderates, e.g., in politics, those who are "very conservative" versus "moderate to very conservative" ones. We propose an activation model of opinion dynamics with interaction rules based on the existence of individual "stubbornness" that mimics empirical observations. According to our modeling, the onset of nonlinearity can be associated to an abrupt bootstrap-percolation transition with cascades of extreme views through society. Therefore, it represents an early-warning signal to forecast the transition from moderate to extreme views. Moreover, by means of a phase diagram we can classify societies according to the percolative regime they belong to, in terms of critical fractions of extremists and people's ties. △ Less

Submitted 15 December, 2014; originally announced December 2014.

Comments: 28 pages, 5 figures

Journal ref: Scientific Reports 5, Article number: 10032 (2015)

Modelling the influence of photospheric turbulence on solar flare statistics

Authors: M. Mendoza, A. Kaydul, L. de Arcangelis, J. S. Andrade Jr., H. J. Herrmann

Abstract: Solar flares stem from the reconnection of twisted magnetic field lines in the solar photosphere. The energy and waiting time distributions of these events follow complex patterns that have been carefully considered in the past and that bear some resemblance with earthquakes and stockmarkets. Here we explore in detail the tangling motion of interacting flux tubes anchored in the plasma and the ene… ▽ More Solar flares stem from the reconnection of twisted magnetic field lines in the solar photosphere. The energy and waiting time distributions of these events follow complex patterns that have been carefully considered in the past and that bear some resemblance with earthquakes and stockmarkets. Here we explore in detail the tangling motion of interacting flux tubes anchored in the plasma and the energy ejections resulting when they recombine. The mechanism for energy accumulation and release in the flow is reminiscent of self-organized criticality. From this model we suggest the origin for two important and widely studied properties of solar flare statistics, including the time-energy correlations. We first propose that the scale-free energy distribution of solar flares is largely due to the twist exerted by the vorticity of the turbulent photosphere. Second, the long-range temporal and time-energy correlations appear to arise from the tube-tube interactions. The agreement with satellite measurements is encouraging. △ Less

Submitted 16 October, 2014; originally announced October 2014.

Comments: 24 pages, 8 figures

Journal ref: Nat. Commun. 5:5035 (2014)

Avoiding catastrophic failure in correlated networks of networks

Authors: Saulo D. S. Reis, Yanqing Hu, Andrés Babino, José S. Andrade Jr., Santiago Canals, Mariano Sigman, Hernán A. Makse

Abstract: Networks in nature do not act in isolation but instead exchange information, and depend on each other to function properly. An incipient theory of Networks of Networks have shown that connected random networks may very easily result in abrupt failures. This theoretical finding bares an intrinsic paradox: If natural systems organize in interconnected networks, how can they be so stable? Here we pro… ▽ More Networks in nature do not act in isolation but instead exchange information, and depend on each other to function properly. An incipient theory of Networks of Networks have shown that connected random networks may very easily result in abrupt failures. This theoretical finding bares an intrinsic paradox: If natural systems organize in interconnected networks, how can they be so stable? Here we provide a solution to this conundrum, showing that the stability of a system of networks relies on the relation between the internal structure of a network and its pattern of connections to other networks. Specifically, we demonstrate that if network inter-connections are provided by hubs of the network and if there is a moderate degree of convergence of inter-network connection the systems of network are stable and robust to failure. We test this theoretical prediction in two independent experiments of functional brain networks (in task- and resting states) which show that brain networks are connected with a topology that maximizes stability according to the theory. △ Less

Submitted 27 April, 2015; v1 submitted 18 September, 2014; originally announced September 2014.

Comments: 40 pages, 7 figures

Journal ref: Nature Physics, Vol 10, 762-767 (2014)

Explosive metallic breakdown

Authors: Cláudio L. N. Oliveira, Nuno A. M. Araújo, José S. Andrade Jr., Hans J. Herrmann

Abstract: We investigate the metallic breakdown of a substrate on which highly conducting particles are adsorbed and desorbed with a probability that depends on the local electric field. We find that, by tuning the relative strength $q$ of this dependence, the breakdown can change from continuous to explosive. Precisely, in the limit in which the adsorption probability is the same for any finite voltage dro… ▽ More We investigate the metallic breakdown of a substrate on which highly conducting particles are adsorbed and desorbed with a probability that depends on the local electric field. We find that, by tuning the relative strength $q$ of this dependence, the breakdown can change from continuous to explosive. Precisely, in the limit in which the adsorption probability is the same for any finite voltage drop, we can map our model exactly onto the $q$-state Potts model and thus the transition to a jump occurs at $q=4$. In another limit, where the adsorption probability becomes independent of the local field strength, the traditional bond percolation model is recovered. Our model is thus an example of a possible experimental realization exhibiting a truly discontinuous percolation transition. △ Less

Submitted 10 July, 2014; originally announced July 2014.

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.

Comments: 8 pages, 5 figures

Searching for superspreaders of information in real-world social media

Authors: Sen Pei, Lev Muchnik, Jose S. Andrade Jr., Zhiming Zheng, Hernan A. Makse

Abstract: A number of predictors have been suggested to detect the most influential spreaders of information in online social media across various domains such as Twitter or Facebook. In particular, degree, PageRank, k-core and other centralities have been adopted to rank the spreading capability of users in information dissemination media. So far, validation of the proposed predictors has been done by simu… ▽ More A number of predictors have been suggested to detect the most influential spreaders of information in online social media across various domains such as Twitter or Facebook. In particular, degree, PageRank, k-core and other centralities have been adopted to rank the spreading capability of users in information dissemination media. So far, validation of the proposed predictors has been done by simulating the spreading dynamics rather than following real information flow in social networks. Consequently, only model-dependent contradictory results have been achieved so far for the best predictor. Here, we address this issue directly. We search for influential spreaders by following the real spreading dynamics in a wide range of networks. We find that the widely-used degree and PageRank fail in ranking users' influence. We find that the best spreaders are consistently located in the k-core across dissimilar social platforms such as Twitter, Facebook, Livejournal and scientific publishing in the American Physical Society. Furthermore, when the complete global network structure is unavailable, we find that the sum of the nearest neighbors' degree is a reliable local proxy for user's influence. Our analysis provides practical instructions for optimal design of strategies for "viral" information dissemination in relevant applications. △ Less

Submitted 10 August, 2014; v1 submitted 7 May, 2014; originally announced May 2014.

Comments: 12 pages, 7 figures

Journal ref: Scientific Reports, 4, 5547 (2014)

Statistical Signs of Social Influence on Suicides

Authors: Hygor Piaget M. Melo, André A. Moreira, Hernán A. Makse, José S. Andrade Jr

Abstract: Certain currents in sociology consider society as being composed of autonomous individuals with independent psychologies. Others, however, deem our actions as strongly influenced by the accepted standards of social behavior. The later view was central to the positivist conception of society when in 1887 Émile Durkheim published his monograph Suicide (Durkheim, 1897). By treating the suicide as a s… ▽ More Certain currents in sociology consider society as being composed of autonomous individuals with independent psychologies. Others, however, deem our actions as strongly influenced by the accepted standards of social behavior. The later view was central to the positivist conception of society when in 1887 Émile Durkheim published his monograph Suicide (Durkheim, 1897). By treating the suicide as a social fact, Durkheim envisaged that suicide rates should be determined by the connections (or the lack of them) between people and society. Under the same framework, Durkheim considered that crime is bound up with the fundamental conditions of all social life and serves a social function. In this sense, and regardless of its extremely deviant nature, crime events are somehow capable to release certain social tensions and so have a purging effect in society. The social effect on the occurrence of homicides has been previously substantiated (Bettencourt et al., 2007; Alves et al., 2013), and confirmed here, in terms of a superlinear scaling relation: by doubling the population of a Brazilian city results in an average increment of 135 % in the number of homicides, rather than the expected isometric increase of 100 %, as found, for example, for the mortality due to car crashes. Here we present statistical signs of the social influence on the suicide occurrence in cities. Differently from homicides (superlinear) and fatal events in car crashes (isometric), we find sublinear scaling behavior between the number of suicides and city population, with allometric power-law exponents, $β= 0.836 \pm 0.009$ and $0.870 \pm 0.002$, for all cities in Brazil and US, respectively. The fact that the frequency of suicides is disproportionately small for larger cities reveals a surprisingly beneficial aspect of living and interacting in larger and more complex social networks. △ Less

Submitted 11 February, 2014; originally announced February 2014.

Comments: 7 pages, 4 figures