-
Unveiling a Family of Dimerized Quantum Magnets in Ternary Metal Borides
Authors:
Zhen Zhang,
Andrew P. Porter,
Yang Sun,
Kirill D. Belashchenko,
Gayatri Viswanathan,
Arka Sarkar,
Kirill Kovnir,
Kai-Ming Ho,
Vladimir Antropov
Abstract:
Dimerized quantum magnets are exotic crystalline materials where Bose-Einstein condensation of magnetic excitations can happen. However, known dimerized quantum magnets are limited to only a few oxides and halides. Here, we unveil 9 dimerized quantum magnets and 11 conventional antiferromagnets in ternary metal borides MTB$_4$ (M = Sc, Y, La, Ce, Lu, Mg, Ca, Al; T = V, Cr, Mn, Fe, Co, Ni). In this…
▽ More
Dimerized quantum magnets are exotic crystalline materials where Bose-Einstein condensation of magnetic excitations can happen. However, known dimerized quantum magnets are limited to only a few oxides and halides. Here, we unveil 9 dimerized quantum magnets and 11 conventional antiferromagnets in ternary metal borides MTB$_4$ (M = Sc, Y, La, Ce, Lu, Mg, Ca, Al; T = V, Cr, Mn, Fe, Co, Ni). In this type of structure, 3d transition-metal atoms T are arranged in dimers. Quantum magnetism in these compounds is dominated by strong antiferromagnetic interactions between Cr (both Cr and Mn for M = Mg and Ca) atoms within the structural dimers, with much weaker interactions between the dimers. These systems are proposed to be close to a quantum critical point between a disordered singlet spin-dimer phase, with a spin gap, and the ordered conventional Néel antiferromagnetic phase. This new family of dimerized quantum magnets greatly enriches the materials inventory that allows investigations of the spin-gap phase. All the quantum-, conventionally-, and non-magnetic systems identified, together with experimental synthesis methods of a phase suitable for characterization, provide a platform with abundant possibilities to tune the magnetic exchange coupling by doping and study this unconventional type of quantum phase transition. This work opens up new avenues for studying the quantum magnetism of spin dimers in borides and establishes a theoretical workflow for future searches for dimerized quantum magnets in other families or types of materials.
△ Less
Submitted 5 June, 2024;
originally announced June 2024.
-
Random minimum spanning tree and dense graph limits
Authors:
Jan Hladký,
Gopal Viswanathan
Abstract:
A theorem of Frieze from 1985 asserts that the total weight of the minimum spanning tree of the complete graph $K_n$ whose edges get independent weights from the distribution $UNIFORM[0,1]$ converges to Apéry's constant in probability, as $n\to\infty$. We generalize this result to sequences of graphs $G_n$ that converge to a graphon $W$. Further, we allow the weights of the edges to be drawn from…
▽ More
A theorem of Frieze from 1985 asserts that the total weight of the minimum spanning tree of the complete graph $K_n$ whose edges get independent weights from the distribution $UNIFORM[0,1]$ converges to Apéry's constant in probability, as $n\to\infty$. We generalize this result to sequences of graphs $G_n$ that converge to a graphon $W$. Further, we allow the weights of the edges to be drawn from different distributions (subject to moderate conditions). The limiting total weight $κ(W)$ of the minimum spanning tree is expressed in terms of a certain branching process defined on $W$, which was studied previously by Bollobás, Janson and Riordan in connection with the giant component in inhomogeneous random graphs.
△ Less
Submitted 18 October, 2023;
originally announced October 2023.
-
A fermionic path integral for exact enumeration of polygons on the simple cubic lattice
Authors:
G. M. Viswanathan
Abstract:
Enumerating polygons on regular lattices is a classic problem in rigorous statistical mechanics. The goal of enumerating polygons on the square lattice via fermionic path integration was achieved using a free-fermion quadratic action in the late 1970s. Given that polygon edges only link 2 vertices, it is considered plausible, if not natural, that an action of degree 2 in the Grassmann variables mi…
▽ More
Enumerating polygons on regular lattices is a classic problem in rigorous statistical mechanics. The goal of enumerating polygons on the square lattice via fermionic path integration was achieved using a free-fermion quadratic action in the late 1970s. Given that polygon edges only link 2 vertices, it is considered plausible, if not natural, that an action of degree 2 in the Grassmann variables might suffice to enumerate lattice polygons in any dimension. Nevertheless, on nonplanar lattices the problem has remained open for more than four decades. Here we derive the Grassmann action for exact enumeration of polygons on the simple cubic lattice. Moreover, we prove that this action is not quadratic but quartic -- corresponding to a model of interacting fermions.
△ Less
Submitted 7 June, 2023; v1 submitted 26 May, 2023;
originally announced May 2023.
-
Connectivity of inhomogeneous random graphs II
Authors:
Jan Hladký,
Gopal Viswanathan
Abstract:
Each graphon $W:Ω^2\rightarrow[0,1]$ yields an inhomogeneous random graph model $G(n,W)$. We show that $G(n,W)$ is asymptotically almost surely connected if and only if (i) $W$ is a connected graphon and (ii) the measure of elements of $Ω$ of $W$-degree less than $α$ is $o(α)$ as $α\rightarrow 0$. These two conditions encapsulate the absence of several linear-sized components, and of isolated vert…
▽ More
Each graphon $W:Ω^2\rightarrow[0,1]$ yields an inhomogeneous random graph model $G(n,W)$. We show that $G(n,W)$ is asymptotically almost surely connected if and only if (i) $W$ is a connected graphon and (ii) the measure of elements of $Ω$ of $W$-degree less than $α$ is $o(α)$ as $α\rightarrow 0$. These two conditions encapsulate the absence of several linear-sized components, and of isolated vertices, respectively.
We study in bigger detail the limit probability of the property that $G(n,W)$ contains an isolated vertex, and, more generally, the limit distribution of the minimum degree of $G(n,W)$.
△ Less
Submitted 5 May, 2023;
originally announced May 2023.
-
Information parity on cortical functional brain networks increases under psychedelic influences
Authors:
Aline Viol,
Gandhi M. Viswanathan,
Oleksandra Soldatkina,
Fernanda Palhano-Fontes,
Heloisa Onias,
Draulio de Araujo,
Philipp Hoevel
Abstract:
The physical basis of consciousness is one of the most intriguing open questions that contemporary science aims to solve. By approaching the brain as an interactive information system, complex network theory has greatly contributed to understand brain process in different states of mind. We study an non-ordinary state of mind by comparing resting-state functional brain networks of individuals in t…
▽ More
The physical basis of consciousness is one of the most intriguing open questions that contemporary science aims to solve. By approaching the brain as an interactive information system, complex network theory has greatly contributed to understand brain process in different states of mind. We study an non-ordinary state of mind by comparing resting-state functional brain networks of individuals in two different conditions: before and after the ingestion of the psychedelic brew Ayahuasca. In order to quantify the functional, statistical symmetries between brain region connectivity, we calculate the pairwise information parity of the functional brain networks. Unlike most of usual network quantification, which is done on a local or global scale, information parity quantifies the pairwise statistical similarities considering the entire network structure. We detect an increase in the average information parity on brain networks of individuals under psychedelic influences. Notably, the information parity between regions from the limbic system and frontal cortex is consistently higher for all the individuals while under the psychedelic influence.
△ Less
Submitted 28 July, 2022;
originally announced July 2022.
-
What does it take to solve the 3D Ising model? Minimal necessary conditions for a valid solution
Authors:
G. M. Viswanathan,
M. A. G. Portillo,
E. P. Raposo,
M. G. E. da Luz
Abstract:
Exact solution of the Ising model on the simple cubic lattice is one of the long-standing open problems in rigorous statistical mechanics. Indeed, it is generally believed that settling it would constitute a methodological breakthrough, fomenting great prospects for further application, similarly to what happened when Lars Onsager solved the two dimensional model eighty years ago. Hence, there hav…
▽ More
Exact solution of the Ising model on the simple cubic lattice is one of the long-standing open problems in rigorous statistical mechanics. Indeed, it is generally believed that settling it would constitute a methodological breakthrough, fomenting great prospects for further application, similarly to what happened when Lars Onsager solved the two dimensional model eighty years ago. Hence, there have been many attempts to find analytic expressions for the exact partition function $Z$, but all such attempts have failed due to unavoidable conceptual or mathematical obstructions. Given the importance of this simple yet paradigmatic model, here we set out clear-cut criteria for any claimed exact expression for $Z$ to be minimally plausible. Specifically, we present six necessary -- but not sufficient -- conditions that $Z$ must satisfy. These criteria will allow very quick plausibility checks of future claims. As illustrative examples, we discuss previous mistaken ``solutions,'' unveiling their shortcomings.
△ Less
Submitted 9 October, 2022; v1 submitted 24 May, 2022;
originally announced May 2022.
-
Numerical studies for an ab initio investigation into the Boltzmann prescription in statistical mechanics of large systems
Authors:
V. Dossetti,
G. M. Viswanathan,
V. M. Kenkre
Abstract:
We present numerical investigations into the question of the validity of the Boltzmann prescription in Statistical Mechanics for large systems, addressing the issue of whether extensivity of energy implies the extensivity of the Boltzmann entropy. The importance of the question stems from the fact that it is currently considered open by some investigators but quite settled by others. We report ab…
▽ More
We present numerical investigations into the question of the validity of the Boltzmann prescription in Statistical Mechanics for large systems, addressing the issue of whether extensivity of energy implies the extensivity of the Boltzmann entropy. The importance of the question stems from the fact that it is currently considered open by some investigators but quite settled by others. We report ab initio results for gas-like Hamiltonian systems with long-range as well as short-range interactions, based on simulations that explicitly consider more than $2^{30} \approx 10^9$ states of the full Hilbert space. The basis of the technique is Monte Carlo algorithms. Despite the largeness of the numbers used, careful inspection shows that the systems studied are still too small to settle uniquely the issues raised. Therefore, the new approach outlined represents a first step in addressing on first principles the question of non-extensive statistical mechanics. General theoretical comments are also supplied to supplement the numerical investigations.
△ Less
Submitted 6 January, 2022;
originally announced January 2022.
-
Spectrum of the tight-binding model on Cayley Trees and comparison with Bethe Lattices
Authors:
M. Ostilli,
Claudionor G. Bezerra,
G. M. Viswanathan
Abstract:
There are few exactly solvable lattice models and even fewer solvable quantum lattice models. Here we address the problem of finding the spectrum of the tight-binding model (equivalently, the spectrum of the adjacency matrix) on Cayley trees. Recent approaches to the problem have relied on the similarity between Cayley tree and the Bethe lattice. Here, we avoid to make any ansatz related to the Be…
▽ More
There are few exactly solvable lattice models and even fewer solvable quantum lattice models. Here we address the problem of finding the spectrum of the tight-binding model (equivalently, the spectrum of the adjacency matrix) on Cayley trees. Recent approaches to the problem have relied on the similarity between Cayley tree and the Bethe lattice. Here, we avoid to make any ansatz related to the Bethe lattice due to fundamental differences between the two lattices that persist even when taking the thermodynamic limit. Instead, we show that one can use a recursive procedure that starts from the boundary and then use the canonical basis to derive the complete spectrum of the tight-binding model on Cayley Trees. Our resulting algorithm is extremely efficient, as witnessed with remarkable large trees having hundred of shells. We also shows that, in the thermodynamic limit, the density of states is dramatically different from that of the Bethe lattice.
△ Less
Submitted 18 March, 2022; v1 submitted 12 June, 2021;
originally announced June 2021.
-
Threshold-free estimation of entropy from a Pearson matrix
Authors:
H. Felippe,
A. Viol,
D. B. de Araujo,
M. G. E. da Luz,
F. Palhano-Fontes,
H. Onias,
E. P. Raposo,
G. M. Viswanathan
Abstract:
There is demand in diverse fields for a reliable method of estimating the entropy associated with correlations. The estimation of a unique entropy directly from the Pearson correlation matrix has remained an open problem for more than half a century. All existing approaches lack generality insofar as they require thresholding choices that arbitrarily remove possibly important information. Here we…
▽ More
There is demand in diverse fields for a reliable method of estimating the entropy associated with correlations. The estimation of a unique entropy directly from the Pearson correlation matrix has remained an open problem for more than half a century. All existing approaches lack generality insofar as they require thresholding choices that arbitrarily remove possibly important information. Here we propose an objective procedure for directly estimating a unique entropy of a general Pearson matrix. We show that upon rescaling the Pearson matrix satisfies all necessary conditions for an analog of the von Neumann entropy to be well defined. No thresholding is required. We demonstrate the method by estimating the entropy from neuroimaging time series of the human brain under the influence of a psychedelic.
△ Less
Submitted 6 February, 2023; v1 submitted 9 June, 2021;
originally announced June 2021.
-
The double hypergeometric series for the partition function of the 2D anisotropic Ising model
Authors:
Gandhimohan M. Viswanathan
Abstract:
In 1944 Lars Onsager published the exact partition function of the ferromagnetic Ising model on the infinite square lattice in terms of a definite integral. Only in the literature of the last decade, however, has it been recast in terms of special functions. Until now all known formulas for the partition function in terms of special functions have been restricted to the important special case of t…
▽ More
In 1944 Lars Onsager published the exact partition function of the ferromagnetic Ising model on the infinite square lattice in terms of a definite integral. Only in the literature of the last decade, however, has it been recast in terms of special functions. Until now all known formulas for the partition function in terms of special functions have been restricted to the important special case of the isotropic Ising model with symmetric couplings. Indeed, the anisotropic model is more challenging because there are two couplings and hence two reduced temperatures, one for each of the two axes of the square lattice. Hence, standard special functions of one variable are inadequate to the task. Here, we reformulate the partition function of the anisotropic Ising model in terms of the Kampé de Fériet function, which is a double hypergeometric function in two variables that is more general than the Appell hypergeometric functions. Finally, we present hypergeometric formulas for the generating function of multipolygons of given length on the infinite square lattice, for isotropic as well as anisotropic edge weights. For the isotropic case, the results allow easy calculation, to arbitrary order, of the celebrated series found by Cyril Domb.
△ Less
Submitted 26 June, 2021; v1 submitted 7 April, 2021;
originally announced April 2021.
-
Comment on "Inverse Square Lévy Walks are not Optimal Search Strategies for d $\geq$ 2" [Phys. Rev. Lett. 124, 080601 (2020)]
Authors:
S. V. Buldyrev,
E. P. Raposo,
F. Bartumeus,
S. Havlin,
F. R. Rusch,
M. G. E. da Luz,
G. M. Viswanathan
Abstract:
It is widely accepted that inverse square Lévy walks are optimal search strategies because they maximize the encounter rate with sparse, randomly distributed, replenishable targets when the search restarts in the vicinity of the previously visited target, which becomes revisitable again with high probability, i.e., non-destructive foraging [Nature 401, 911 (1999)]. The precise conditions for the v…
▽ More
It is widely accepted that inverse square Lévy walks are optimal search strategies because they maximize the encounter rate with sparse, randomly distributed, replenishable targets when the search restarts in the vicinity of the previously visited target, which becomes revisitable again with high probability, i.e., non-destructive foraging [Nature 401, 911 (1999)]. The precise conditions for the validity of this Lévy flight foraging hypothesis (LFH) have been widely described in the literature [Phys. Life Rev. 14, 94 (2015)]. Nevertheless, three objecting claims to the LFH have been raised recently for $d \geq 2$: (i) the capture rate $η$ has linear dependence on the target density $ρ$ for all values of the Lévy index $α$; (ii) "the gain $η_{max}/η$ achieved by varying $α$ is bounded even in the limit $ρ\to 0 $" so that "tuning $α$ can only yield a marginal gain"; (iii) depending on the values of the radius of detection $a$, the restarting distance $l_c$ and the scale parameter $s$, the optimum is realized for a range of $α$ [Phys. Rev. Lett. 124, 080601 (2020)]. Here we answer each of these three criticisms in detail and show that claims (i)-(iii) do not actually invalidate the LFH. Our results and analyses restore the original result of the LFH for non-destructive foraging.
△ Less
Submitted 22 March, 2021; v1 submitted 19 March, 2021;
originally announced March 2021.
-
Eclipse timing variation of GK Vir: evidence of a possible Jupiter-like planet in a circumbinary orbit
Authors:
Leonardo A. Almeida,
Elielson S. Pereira,
Gislene M. Borges,
Augusto Damineli,
Tatiana A. Michtchenko,
Gandhi M. Viswanathan
Abstract:
Eclipse timing variation analysis has become a powerful method to discover planets around binary systems. We applied this technique to investigate the eclipse times of GK Vir. This system is a post-common envelope binary with an orbital period of 8.26 h. Here, we present 10 new eclipse times obtained between 2013 and 2020. We calculated the O-C diagram using a linear ephemeris and verified a clear…
▽ More
Eclipse timing variation analysis has become a powerful method to discover planets around binary systems. We applied this technique to investigate the eclipse times of GK Vir. This system is a post-common envelope binary with an orbital period of 8.26 h. Here, we present 10 new eclipse times obtained between 2013 and 2020. We calculated the O-C diagram using a linear ephemeris and verified a clear orbital period variation (OPV) with a cyclic behavior. We investigated if this variation could be explained by the Applegate mechanism, the apsidal motion, or the light travel time (LTT) effect. We found that the Applegate mechanism would hardly explain the OPV with its current theoretical description. We obtained using different approaches that the apsidal motion is a less likely explanation than the LTT effect. We showed that the LTT effect with one circumbinary body is the most likely cause for the OPV, which was reinforced by the orbital stability of the third body. The LTT best solution provided an orbital period of ~24 yr for the outer body. Under the assumption of coplanarity between the external body and the inner binary, we obtained a Jupiter-like planet around the GK Vir. In this scenario, the planet has one of the longest orbital periods, with a full observational baseline, discovered so far. However, as the observational baseline of GK Vir is smaller than twice the period found in the O-C diagram, the LTT solution must be taken as preliminary.
△ Less
Submitted 8 September, 2020;
originally announced September 2020.
-
The connection between Jackson and Hausdorff derivatives in the context of generalized statistical mechanics
Authors:
Andre A. Marinho,
G. M. Viswanathan,
Francisco A. Brito,
C. G. Bezerra
Abstract:
In literature one can find many generalizations of the usual Leibniz derivative, such as Jackson derivative, Tsallis derivative and Hausdorff derivative. In this article we present a connection between Jackson derivative and recently proposed Hausdorff derivative. On one hand, the Hausdorff derivative has been previously associated with non-extensivity in systems presenting fractal aspects. On the…
▽ More
In literature one can find many generalizations of the usual Leibniz derivative, such as Jackson derivative, Tsallis derivative and Hausdorff derivative. In this article we present a connection between Jackson derivative and recently proposed Hausdorff derivative. On one hand, the Hausdorff derivative has been previously associated with non-extensivity in systems presenting fractal aspects. On the other hand, the Jackson derivative has a solid mathematical basis because it is the $\overline{q}$-analog of the ordinary derivative and it also arises in quantum calculus. From a quantum deformed $\overline{q}$-algebra we obtain the Jackson derivative and then address the problem of $N$ non-interacting quantum oscillators. We perform an expansion in the quantum grand partition function from which we obtain a relationship between the parameter $\overline{q}$, related to Jackson derivative, and the parameters $ζ$ and $q$ related to Hausdorff derivative and Tsallis derivative, respectively.
△ Less
Submitted 30 May, 2020;
originally announced June 2020.
-
An analogue of Pillai's theorem for continued fraction normality and an application to subsequences
Authors:
Satyadev Nandakumar,
Subin Pulari,
Prateek Vishnoi,
Gopal Viswanathan
Abstract:
We show that two notions of continued fraction normality, one where overlapping occurrences of finite patterns are counted as distinct occurrences, and another where only disjoint occurrences are counted as distinct, are identical. This equivalence involves an analogue of a theorem due to S. S. Pillai in 1940 for base-$b$ expansions. The proof requires techniques which are fundamentally different,…
▽ More
We show that two notions of continued fraction normality, one where overlapping occurrences of finite patterns are counted as distinct occurrences, and another where only disjoint occurrences are counted as distinct, are identical. This equivalence involves an analogue of a theorem due to S. S. Pillai in 1940 for base-$b$ expansions. The proof requires techniques which are fundamentally different, since the continued fraction expansion utilizes a countably infinite alphabet, leading to a non-compact space. Utilizing the equivalence of these two notions, we provide a new proof of Heersink and Vandehey's recent result that selection of subsequences along arithmetic progressions does not preserve continued fraction normality.
△ Less
Submitted 10 September, 2019; v1 submitted 8 September, 2019;
originally announced September 2019.
-
Characterizing complex networks using Entropy-degree diagrams: unveiling changes in functional brain connectivity induced by Ayahuasca
Authors:
A. Viol,
Fernanda Palhano-Fontes,
Heloisa Onias,
Draulio B. de Araujo,
Philipp Hövel,
G. M. Viswanathan
Abstract:
Open problems abound in the theory of complex networks, which has found successful application to diverse fields of science. With the aim of further advancing the understanding of the brain's functional connectivity, we propose to evaluate a network metric which we term the geodesic entropy. This entropy, in a way that can be made precise, quantifies the Shannon entropy of the distance distributio…
▽ More
Open problems abound in the theory of complex networks, which has found successful application to diverse fields of science. With the aim of further advancing the understanding of the brain's functional connectivity, we propose to evaluate a network metric which we term the geodesic entropy. This entropy, in a way that can be made precise, quantifies the Shannon entropy of the distance distribution to a specific node from all other nodes. Measurements of geodesic entropy allow for the characterization of the structural information of a network that takes into account the distinct role of each node into the network topology. The measurement and characterization of this structural information has the potential to greatly improve our understanding of sustained activity and other emergent behaviors in networks, such as self-organized criticality sometimes seen in such contexts. We apply these concepts and methods to study the effects of how the psychedelic Ayahuasca affects the functional connectivity of the human brain. We show that the geodesic entropy is able to differentiate the functional networks of the human brain in two different states of consciousness in the resting state: (i) the ordinary waking state and (ii) a state altered by ingestion of the Ayahuasca. The entropy of the nodes of brain networks from subjects under the influence of Ayahuasca diverge significantly from those of the ordinary waking state. The functional brain networks from subjects in the altered state have, on average, a larger geodesic entropy compared to the ordinary state. We conclude that geodesic entropy is a useful tool for analyzing complex networks and discuss how and why it may bring even further valuable insights into the study of the human brain and other empirical networks.
△ Less
Submitted 26 September, 2018;
originally announced September 2018.
-
Correspondence between spanning trees and the Ising model on a square lattice
Authors:
G. M. Viswanathan
Abstract:
An important problem in statistical physics concerns the fascinating connections between partition functions of lattice models studied in equilibrium statistical mechanics on the one hand and graph theoretical enumeration problems on the other hand. We investigate the nature of the relationship between the number of spanning trees and the partition function of the Ising model on the square lattice…
▽ More
An important problem in statistical physics concerns the fascinating connections between partition functions of lattice models studied in equilibrium statistical mechanics on the one hand and graph theoretical enumeration problems on the other hand. We investigate the nature of the relationship between the number of spanning trees and the partition function of the Ising model on the square lattice. The spanning tree generating function $T(z)$ gives the spanning tree constant when evaluated at $z=1$, while giving he lattice green function when differentiated. It is known that for the infinite square lattice the partition function $Z(K)$ of the Ising model evaluated at the critical temperature $K=K_c$ is related to $T(1)$. Here we show that this idea in fact generalizes to all real temperatures. We prove that $ ( Z(K) {\rm sech~} 2K ~\!)^2 = k \exp\big[ T(k) \big] $, where $k= 2 \tanh(2K) {\rm sech}(2K)$. The identical Mahler measure connects the two seemingly disparate quantities $T(z)$ and $Z(K)$. In turn, the Mahler measure is determined by the random walk structure function. Finally, we show that the the above correspondence does not generalize in a straightforward manner to non-planar lattices.
△ Less
Submitted 6 May, 2020; v1 submitted 2 June, 2017;
originally announced June 2017.
-
The complex social network of surnames: A comparison between Brazil and Portugal
Authors:
G. D. Ferreira,
G. M. Viswanathan,
L. R. da Silva,
H. J. Herrmann
Abstract:
We present a study of social networks based on the analysis of Brazilian and Portuguese family names (surnames). We construct networks whose nodes are names of families and whose edges represent parental relations between two families. From these networks we extract the connectivity distribution, clustering coefficient, shortest path and centrality. We find that the connectivity distribution follo…
▽ More
We present a study of social networks based on the analysis of Brazilian and Portuguese family names (surnames). We construct networks whose nodes are names of families and whose edges represent parental relations between two families. From these networks we extract the connectivity distribution, clustering coefficient, shortest path and centrality. We find that the connectivity distribution follows an approximate power law. We associate the number of hubs, centrality and entropy to the degree of miscegenation in the societies in both countries. Our results show that Portuguese society has a higher miscegenation degree than Brazilian society. All networks analyzed lead to approximate inverse square power laws in the degree distribution. We conclude that the thermodynamic limit is reached for small networks (3 or 4 thousand nodes). The assortative mixing of all networks is negative, showing that the more connected vertices are connected to vertices with lower connectivity. Finally, the network of surnames presents some small world characteristics.
△ Less
Submitted 12 May, 2017;
originally announced May 2017.
-
Shannon entropy of brain functional complex networks under the influence of the psychedelic Ayahuasca
Authors:
A. Viol,
Fernanda Palhano-Fontes,
Heloisa Onias,
Draulio B. de Araujo,
G. M. Viswanathan
Abstract:
The entropic brain hypothesis holds that the key facts concerning psychedelics are partially explained in terms of increased entropy of the brain's functional connectivity. Ayahuasca is a psychedelic beverage of Amazonian indigenous origin with legal status in Brazil in religious and scientific settings. In this context, we use tools and concepts from the theory of complex networks to analyze rest…
▽ More
The entropic brain hypothesis holds that the key facts concerning psychedelics are partially explained in terms of increased entropy of the brain's functional connectivity. Ayahuasca is a psychedelic beverage of Amazonian indigenous origin with legal status in Brazil in religious and scientific settings. In this context, we use tools and concepts from the theory of complex networks to analyze resting state fMRI data of the brains of human subjects under two distinct conditions: (i) under ordinary waking state and (ii) in an altered state of consciousness induced by ingestion of Ayahuasca. We report an increase in the Shannon entropy of the degree distribution of the networks subsequent to Ayahuasca ingestion. We also find increased local and decreased global network integration. Our results are broadly consistent with the entropic brain hypothesis. Finally, we discuss our findings in the context of descriptions of "mind-expansion" frequently seen in self-reports of users of psychedelic drugs.
△ Less
Submitted 1 November, 2016;
originally announced November 2016.
-
Information entropy of classical versus explosive percolation
Authors:
T. M. Vieira,
G. M. Viswanathan,
L. R. da Silva
Abstract:
We study the Shannon entropy of the cluster size distribution in classical as well as explosive percolation, in order to estimate the uncertainty in the sizes of randomly chosen clusters. At the critical point the cluster size distribution is a power-law, i.e. there are clusters of all sizes, so one expects the information entropy to attain a maximum. As expected, our results show that the entropy…
▽ More
We study the Shannon entropy of the cluster size distribution in classical as well as explosive percolation, in order to estimate the uncertainty in the sizes of randomly chosen clusters. At the critical point the cluster size distribution is a power-law, i.e. there are clusters of all sizes, so one expects the information entropy to attain a maximum. As expected, our results show that the entropy attains a maximum at this point for classical percolation. Surprisingly, for explosive percolation the maximum entropy does not match the critical point. Moreover, we show that it is possible determine the critical point without using the conventional order parameter, just analysing the entropy's derivatives.
△ Less
Submitted 14 August, 2015; v1 submitted 21 November, 2014;
originally announced November 2014.
-
The hypergeometric series for the partition function of the 2-D Ising model
Authors:
G. M. Viswanathan
Abstract:
In 1944 Onsager published the formula for the partition function of the Ising model for the infinite square lattice. He was able to express the internal energy in terms of a special function, but he left the free energy as a definite integral. Seven decades later, the partition function and free energy have yet to be written in closed form, even with the aid of special functions. Here we evaluate…
▽ More
In 1944 Onsager published the formula for the partition function of the Ising model for the infinite square lattice. He was able to express the internal energy in terms of a special function, but he left the free energy as a definite integral. Seven decades later, the partition function and free energy have yet to be written in closed form, even with the aid of special functions. Here we evaluate the definite integral explicitly, using hypergeometric series. Let $β$ denote the reciprocal temperature, $J$ the coupling and $f$ the free energy per spin. We prove that $-βf = \ln(2 \cosh 2K) - κ^2\, {}_4F_3 [1,1,\tfrac{3}{2},\tfrac{3}{2};\ 2,2,2 ;\ 16 κ^2 ] $, where $_p F_q$ is the generalized hypergeometric function, $K=βJ$, and $2κ= {\rm tanh} 2K {\rm sech} 2K$.
△ Less
Submitted 24 July, 2015; v1 submitted 10 November, 2014;
originally announced November 2014.
-
How to efficiently destroy a network with limited information
Authors:
T. M. Vieira,
G. M. Viswanathan,
L. R. da Silva
Abstract:
We address the general problem of how best to attack and destroy a network by node removal, given limited or no prior information about the edges. We consider a family of strategies in which nodes are randomly chosen, but not removed. Instead, a random acquaintance (i.e., a first neighbour) of the chosen node is removed from the network. By assigning an informal cost to the information about the n…
▽ More
We address the general problem of how best to attack and destroy a network by node removal, given limited or no prior information about the edges. We consider a family of strategies in which nodes are randomly chosen, but not removed. Instead, a random acquaintance (i.e., a first neighbour) of the chosen node is removed from the network. By assigning an informal cost to the information about the network structure, we show using cost-benefit analysis that acquaintance removal is the optimal strategy to destroy networks efficiently.
△ Less
Submitted 10 November, 2014;
originally announced November 2014.
-
Activity, diffusion, and correlations in a two-dimensional conserved stochastic sandpile
Authors:
Sharon Dantas da Cunha,
Luciano Rodrigues da Silva,
Gandhimohan M. Viswanathan,
Ronald Dickman
Abstract:
We perform large-scale simulations of a two-dimensional restricted-height conserved stochastic sandpile, focusing on particle diffusion and mobility, and spatial correlations. Quasistationary (QS) simulations yield the critical particle density to high precision [$p_c = 0.7112687(2)$], and show that the diffusion constant scales in the same manner as the activity density, as found previously in th…
▽ More
We perform large-scale simulations of a two-dimensional restricted-height conserved stochastic sandpile, focusing on particle diffusion and mobility, and spatial correlations. Quasistationary (QS) simulations yield the critical particle density to high precision [$p_c = 0.7112687(2)$], and show that the diffusion constant scales in the same manner as the activity density, as found previously in the one-dimensional case. Short-time scaling is characterized by subdiffusive behavior (mean-square displacement $\sim t^γ$ with $γ< 1$), which is easily understood as a consequence of the initial decay of activity, $ρ(t) \sim t^{-δ}$, with $γ= 1- δ$. We verify that at criticality, the activity correlation function $C(r) \sim r^{-β/ν_\perp}$, as expected at an absorbing-state phase transition. Our results for critical exponents are consistent with, and somewhat more precise than, predictions derived from the Langevin equation for stochastic sandpiles in two dimensions.
△ Less
Submitted 5 May, 2014;
originally announced May 2014.
-
A Solvability criterion for Navier-Stokes equations in high dimensions
Authors:
T. M. Viswanathan,
G. M. Viswanathan
Abstract:
We define the Ladyzhenskaya-Lions exponent $α_{\rm {\tiny
\sc l}} (n)=({2+n})/4$ for Navier-Stokes equations with dissipation $-(-Δ)^α$ in ${\Bbb R}^n$, for all $n\geq 2$. We review the proof of strong global solvability when $α\geq α_{\rm
{\tiny \sc l}} (n)$, given smooth initial data. If the corresponding Euler equations for $n>2$ were to allow uncontrolled growth of the enstrophy…
▽ More
We define the Ladyzhenskaya-Lions exponent $α_{\rm {\tiny
\sc l}} (n)=({2+n})/4$ for Navier-Stokes equations with dissipation $-(-Δ)^α$ in ${\Bbb R}^n$, for all $n\geq 2$. We review the proof of strong global solvability when $α\geq α_{\rm
{\tiny \sc l}} (n)$, given smooth initial data. If the corresponding Euler equations for $n>2$ were to allow uncontrolled growth of the enstrophy ${1\over 2} \|\nabla u \|^2_{L^2}$, then no globally controlled coercive quantity is currently known to exist that can regularize solutions of the Navier-Stokes equations for $α<α_{\rm {\tiny \sc l}} (n)$. The energy is critical under scale transformations only for $α=α_{\rm {\tiny
\sc l}} (n)$.
△ Less
Submitted 16 October, 2009; v1 submitted 24 July, 2009;
originally announced July 2009.
-
Spontaneous symmetry breaking and finite time singularities in $d$-dimensional incompressible flow with fractional dissipation
Authors:
G. M. Viswanathan,
T. M. Viswanathan
Abstract:
We investigate the formation of singularities in the incompressible Navier-Stokes equations in $d\geq 2$ dimensions with a fractional Laplacian $|\nabla |^α$. We derive analytically a sufficient but not necessary condition for solutions to remain always smooth and show that finite time singularities cannot form for $α\geq α_c= 1+d/2$. Moreover, initial singularities become unstable for $α>α_c$.
We investigate the formation of singularities in the incompressible Navier-Stokes equations in $d\geq 2$ dimensions with a fractional Laplacian $|\nabla |^α$. We derive analytically a sufficient but not necessary condition for solutions to remain always smooth and show that finite time singularities cannot form for $α\geq α_c= 1+d/2$. Moreover, initial singularities become unstable for $α>α_c$.
△ Less
Submitted 9 July, 2008;
originally announced July 2008.
-
Universal photocurrent-voltage characteristics of dye sensitized nanocrystalline TiO$_2$ photoelectrochemical cells
Authors:
J. S. Agnaldo,
J. C. Cressoni,
G. M. Viswanathan
Abstract:
We propose a new linearizable model for the nonlinear photocurrent-voltage characteristics of nanocrystalline TiO$_2$ dye sensitized solar cells based on first principles and report predicted values for fill factors. Upon renormalization diverse experimental photocurrent-voltage data collapse onto a single universal function. These advances allow the estimation of the complete current-voltage cu…
▽ More
We propose a new linearizable model for the nonlinear photocurrent-voltage characteristics of nanocrystalline TiO$_2$ dye sensitized solar cells based on first principles and report predicted values for fill factors. Upon renormalization diverse experimental photocurrent-voltage data collapse onto a single universal function. These advances allow the estimation of the complete current-voltage curve and the fill factor from any three experimental data points, e.g., the open circuit voltage, the short circuit current and one intermediate measurement. The theoretical underpinning provides insight into the physical mechanisms responsible for the remarkably large fill factors as well as their known dependence on the open circuit voltage.
△ Less
Submitted 28 January, 2008;
originally announced January 2008.
-
Spontaneous symmetry breaking in amnestically induced persistence
Authors:
Marco Antonio Alves da Silva,
A. S. Ferreira,
G. M. Viswanathan,
J. C. Cressoni
Abstract:
We investigate a recently proposed non-Markovian random walk model characterized by loss of memories of the recent past and amnestically induced persistence. We report numerical and analytical results showing the complete phase diagram, consisting of 4 phases, for this system: (i) classical nonpersistence, (ii) classical persistence (iii) log-periodic nonpersistence and (iv) log-periodic persist…
▽ More
We investigate a recently proposed non-Markovian random walk model characterized by loss of memories of the recent past and amnestically induced persistence. We report numerical and analytical results showing the complete phase diagram, consisting of 4 phases, for this system: (i) classical nonpersistence, (ii) classical persistence (iii) log-periodic nonpersistence and (iv) log-periodic persistence driven by negative feedback. The first two phases possess continuous scale invariance symmetry, however log-periodicity breaks this symmetry. Instead, log-periodic motion satisfies discrete scale invariance symmetry, with complex rather than real fractal dimensions. We find for log-periodic persistence evidence not only of statistical but also of geometric self-similarity.
△ Less
Submitted 11 September, 2007; v1 submitted 22 August, 2007;
originally announced August 2007.
-
The origin of power-law distributions in deterministic walks: the influence of landscape geometry
Authors:
M. C. Santos,
D. Boyer,
O. Miramontes,
G. M. Viswanathan,
E. P. Raposo,
J. L. Mateos,
M. G. E. da Luz
Abstract:
We investigate the properties of a deterministic walk, whose locomotion rule is always to travel to the nearest site. Initially the sites are randomly distributed in a closed rectangular ($A/L \times L)$ landscape and, once reached, they become unavailable for future visits. As expected, the walker step lengths present characteristic scales in one ($L \to 0$) and two ($A/L \sim L$) dimensions. H…
▽ More
We investigate the properties of a deterministic walk, whose locomotion rule is always to travel to the nearest site. Initially the sites are randomly distributed in a closed rectangular ($A/L \times L)$ landscape and, once reached, they become unavailable for future visits. As expected, the walker step lengths present characteristic scales in one ($L \to 0$) and two ($A/L \sim L$) dimensions. However, we find scale invariance for an intermediate geometry, when the landscape is a thin strip-like region. This result is induced geometrically by a dynamical trapping mechanism, leading to a power law distribution for the step lengths. The relevance of our findings in broader contexts -- of both deterministic and random walks -- is also briefly discussed.
△ Less
Submitted 30 May, 2007;
originally announced May 2007.
-
Forecasting extreme events in collective dynamics: an analytic signal approach to detecting discrete scale invariance
Authors:
G. M. Viswanathan
Abstract:
A challenging problem in physics concerns the possibility of forecasting rare but extreme phenomena such as large earthquakes, financial market crashes, and material rupture. A promising line of research involves the early detection of precursory log-periodic oscillations to help forecast extreme events in collective phenomena where discrete scale invariance plays an important role. Here I inves…
▽ More
A challenging problem in physics concerns the possibility of forecasting rare but extreme phenomena such as large earthquakes, financial market crashes, and material rupture. A promising line of research involves the early detection of precursory log-periodic oscillations to help forecast extreme events in collective phenomena where discrete scale invariance plays an important role. Here I investigate two distinct approaches towards the general problem of how to detect log-periodic oscillations in arbitrary time series without prior knowledge of the location of the moveable singularity. I first show that the problem has a definite solution in Fourier space, however the technique involved requires an unrealistically large signal to noise ratio. I then show that the quadrature signal obtained via analytic continuation onto the imaginary axis, using the Hilbert transform, necessarily retains the log-periodicities found in the original signal. This finding allows the development of a new method of detecting log-periodic oscillations that relies on calculation of the instantaneous phase of the analytic signal. I illustrate the method by applying it to the well documented stock market crash of 1987. Finally, I discuss the relevance of these findings for parametric rather than nonparametric estimation of critical times.
△ Less
Submitted 4 December, 2006; v1 submitted 28 November, 2006;
originally announced November 2006.
-
Amnestically induced persistence in random walks
Authors:
J. C. Cressoni,
M. A. A. da Silva,
G. M. Viswanathan
Abstract:
We study how the Hurst exponent $α$ depends on the fraction $f$ of the total time $t$ remembered by non-Markovian random walkers that recall only the distant past. We find that otherwise nonpersistent random walkers switch to persistent behavior when inflicted with significant memory loss. Such memory losses induce the probability density function of the walker's position to undergo a transition…
▽ More
We study how the Hurst exponent $α$ depends on the fraction $f$ of the total time $t$ remembered by non-Markovian random walkers that recall only the distant past. We find that otherwise nonpersistent random walkers switch to persistent behavior when inflicted with significant memory loss. Such memory losses induce the probability density function of the walker's position to undergo a transition from Gaussian to non-Gaussian. We interpret these findings of persistence in terms of a breakdown of self-regulation mechanisms and discuss their possible relevance to some of the burdensome behavioral and psychological symptoms of Alzheimer's disease and other dementias.
△ Less
Submitted 6 February, 2007; v1 submitted 17 November, 2006;
originally announced November 2006.
-
Variance fluctuations in nonstationary time series: a comparative study of music genres
Authors:
Heather D. Jennings,
Plamen Ch. Ivanov,
A. M. Martins,
P. C. da Silva,
G. M. Viswanathan
Abstract:
An important problem in physics concerns the analysis of audio time series generated by transduced acoustic phenomena. Here, we develop a new method to quantify the scaling properties of the local variance of nonstationary time series. We apply this technique to analyze audio signals obtained from selected genres of music. We find quantitative differences in the correlation properties of high ar…
▽ More
An important problem in physics concerns the analysis of audio time series generated by transduced acoustic phenomena. Here, we develop a new method to quantify the scaling properties of the local variance of nonstationary time series. We apply this technique to analyze audio signals obtained from selected genres of music. We find quantitative differences in the correlation properties of high art music, popular music, and dance music. We discuss the relevance of these objective findings in relation to the subjective experience of music.
△ Less
Submitted 21 December, 2003; v1 submitted 15 December, 2003;
originally announced December 2003.
-
Critical behavior of an epidemic model of drug resistant diseases
Authors:
C. R. da Silva,
U. L. Fulco,
M. L. Lyra,
G. M. Viswanathan
Abstract:
In this work, we study the critical behavior of an epidemic propagation model that considers individuals that can develop drug resistance. In our lattice model, each site can be found in one of four states: empty, healthy, normally infected (not drug resistant) and strain infected (drug resistant) states. The most relevant parameters in our model are related to the mortality, cure and mutation r…
▽ More
In this work, we study the critical behavior of an epidemic propagation model that considers individuals that can develop drug resistance. In our lattice model, each site can be found in one of four states: empty, healthy, normally infected (not drug resistant) and strain infected (drug resistant) states. The most relevant parameters in our model are related to the mortality, cure and mutation rates. This model presents two distinct stationary active phases: a phase with co-existing normal and drug resistant infected individuals and an intermediate active phase with only drug resistant individuals. We employ a finite-size scaling analysis to compute the critical points the critical exponents ratio $β/ν$ governing the phase-transitions between these active states and the absorbing inactive state. Our results are consistent with the hypothesis that these transitions belong to the directed percolation universality class.
△ Less
Submitted 28 November, 2003;
originally announced December 2003.
-
Kinematics of stock prices
Authors:
M. Serva,
U. L. Fulco,
M. L. Lyra,
G. M. Viswanathan
Abstract:
We investigate the general problem of how to model the kinematics of stock prices without considering the dynamical causes of motion. We propose a stochastic process with long-range correlated absolute returns. We find that the model is able to reproduce the experimentally observed clustering, power law memory, fat tails and multifractality of real financial time series. We find that the distrib…
▽ More
We investigate the general problem of how to model the kinematics of stock prices without considering the dynamical causes of motion. We propose a stochastic process with long-range correlated absolute returns. We find that the model is able to reproduce the experimentally observed clustering, power law memory, fat tails and multifractality of real financial time series. We find that the distribution of stock returns is approximated by a Gaussian with log-normally distributed local variance and shows excellent agreement with the behavior of the NYSE index for a range of time scales.
△ Less
Submitted 4 September, 2002;
originally announced September 2002.
-
The origin of fat tailed distributions in financial time series
Authors:
G. M. Viswanathan,
U. L. Fulco,
M. L. Lyra,
M. Serva
Abstract:
A classic problem in physics is the origin of fat tailed distributions generated by complex systems. We study the distributions of stock returns measured over different time lags $τ.$ We find that destroying all correlations without changing the $τ= 1$ d distribution, by shuffling the order of the daily returns, causes the fat tails almost to vanish for $τ>1$ d. We argue that the fat tails are c…
▽ More
A classic problem in physics is the origin of fat tailed distributions generated by complex systems. We study the distributions of stock returns measured over different time lags $τ.$ We find that destroying all correlations without changing the $τ= 1$ d distribution, by shuffling the order of the daily returns, causes the fat tails almost to vanish for $τ>1$ d. We argue that the fat tails are caused by known long-range volatility correlations. Indeed, destroying only sign correlations, by shuffling the order of only the signs (but not the absolute values) of the daily returns, allows the fat tails to persist for $τ>1$ d.
△ Less
Submitted 6 November, 2002; v1 submitted 28 December, 2001;
originally announced December 2001.
-
Average Time Spent by Levy Flights and Walks on an Interval with Absorbing Boundaries
Authors:
S. V. Buldyrev,
S. Havlin,
A. Ya. Kazakov,
M. G. E. da Luz,
E. P. Raposo,
H. E. Stanley,
G. M. Viswanathan
Abstract:
We consider a Levy flyer of order alpha that starts from a point x0 on an interval [O,L] with absorbing boundaries. We find a closed-form expression for the average number of flights the flyer takes and the total length of the flights it travels before it is absorbed. These two quantities are equivalent to the mean first passage times for Levy flights and Levy walks, respectively. Using fraction…
▽ More
We consider a Levy flyer of order alpha that starts from a point x0 on an interval [O,L] with absorbing boundaries. We find a closed-form expression for the average number of flights the flyer takes and the total length of the flights it travels before it is absorbed. These two quantities are equivalent to the mean first passage times for Levy flights and Levy walks, respectively. Using fractional differential equations with a Riesz kernel, we find exact analytical expressions for both quantities in the continuous limit. We show that numerical solutions for the discrete Levy processes converge to the continuous approximations in all cases except the case of alpha approaching 2 and the cases of x0 near absorbing boundaries. For alpha larger than 2 when the second moment of the flight length distribution exists, our result is replaced by known results of classical diffusion. We show that if x0 is placed in the vicinity of absorbing boundaries, the average total length has a minimum at alpha=1, corresponding to the Cauchy distribution. We discuss the relevance of this result to the problem of foraging, which has received recent attention in the statistical physics literature.
△ Less
Submitted 20 May, 2001; v1 submitted 29 December, 2000;
originally announced December 2000.
-
The Strange Behavior of Critical Branched Polymers
Authors:
H. H. Aragao-Rego,
J. E. de Freitas,
Liacir S. Lucena,
G. M. Viswanathan
Abstract:
We find that 2-dimensional (2-D) critical branched polymers with no impurities conclusively belong to the same universality class as 2-D random percolation clusters, although pure critical 3-D branched polymers do not belong to the 3-D percolation universality class. We find, moreover, that the fractal dimension of critical branched polymers in 2-D is independent of the presence of a random envi…
▽ More
We find that 2-dimensional (2-D) critical branched polymers with no impurities conclusively belong to the same universality class as 2-D random percolation clusters, although pure critical 3-D branched polymers do not belong to the 3-D percolation universality class. We find, moreover, that the fractal dimension of critical branched polymers in 2-D is independent of the presence of a random environment, and explore possible explanations for these strange results. We also report that when there are no impurities the critical branching probability in 3-D is bc=3.34e-4 +- 0.16e-4.
△ Less
Submitted 10 August, 2000; v1 submitted 6 July, 2000;
originally announced July 2000.
-
Quantifying Nonstationary Radioactivity Concentration Fluctuations Near Chernobyl: A Complete Statistical Description
Authors:
G. M. Viswanathan,
S. V. Buldyrev,
E. K. Garger,
V. A. Kashpur,
L. S. Lucena,
A. Shlyakhter,
H. E. Stanley,
J. Tschiersch
Abstract:
We analyze nonstationary $^{137}$Cs atmospheric activity concentration fluctuations measured near Chernobyl after the 1986 disaster and find three new results: (i) the histogram of fluctuations is well described by a log-normal distribution, (ii) there is a pronounced spectral component with period T=1 y, and (iii) the fluctuations are long-range correlated. These findings allow us to quantify t…
▽ More
We analyze nonstationary $^{137}$Cs atmospheric activity concentration fluctuations measured near Chernobyl after the 1986 disaster and find three new results: (i) the histogram of fluctuations is well described by a log-normal distribution, (ii) there is a pronounced spectral component with period T=1 y, and (iii) the fluctuations are long-range correlated. These findings allow us to quantify two fundamental statistical properties of the data - the probability distribution and the correlation properties of the time series. We interpret our findings as evidence that the atmospheric radionuclide resuspension processes are tightly coupled to the surrounding ecosystems and to large time scale weather patterns.
△ Less
Submitted 19 May, 2000;
originally announced May 2000.