LIVRO DE RESUMOS
Encontro Nacional
ENFE 15 de Física Estatística 1-4 de Novembro de 2015 Vitória-ES
COMITÊ NACIONAL Roberto Andrade (UFBA) Ronald Dickman (UFMG) Dora Izzo (UFRJ) Tarcísio Marciano (UnB) João Plascak (UFPB) Carmen Prado (USP) Silvio Salinas (USP) José Soares Andrade Jr. (UFC) Constantino Tsallis (CBPF)
COMITÊ ORGANIZADOR Celia Anteneodo (PUC-Rio) Jeferson Arenzon (UFRGS) Marcia Barbosa (UFRGS) Marcos G. E. da Luz (UFPR) Gandhi Viswanathan (UFRN)
Local do evento: Sheraton Vitoria Hotel, Av. Saturnino de Brito 217, Praia do Canto, Vitória, ES Mais informações: http://enfe.fis.puc-rio.br Contato: [email protected]
UFES-CCE
1 Plenary talks (8:30 - 10:00) Posters (10:00-11:00 / 16:05-17:05)
PLENARY SESSION COMPLEXITY & INTERDISCIPLI- NARY TOPICS [02/11/2015 - 09:00 - Room Vit´oria] A random walk picture of basketball scoring and [02/11/2015 - P001] lead-change dynamics, Aaron Clauset, Marina Spatial organization and mobility effects in collec- tive hunting and defense strategies of predator- Kogan, Sid Redner, Santa Fe Institute By analy- ´ zing recently available play-by-play data from all regular- prey systems, Annette Cazaubiel, Ecole Normale season games from multiple seasons of the National Bas- Sup´erieure - International Center of Fundamental Phy- ketball Association (NBA), we present evidence that, sics, Paris - France, Alessandra F. Lutz¨ , Jeferson basketball scoring during a game is well described by a J. Arenzon, Instituto de F´ısica - UFRGS , Porto Ale- continuous-time anti-persistent random walk. The time gre RS - Brasil There is a myriad of strategies that intervals between successive scoring events follow an ex- predators utilize to increase their rate of success. Among ponential distribution, with essentially no correlations them, preys may be attacked in a cooperative, coordina- between different scoring intervals. We will also argue ted way, these actions being correlated in space and time. that the heterogeneity of team strengths plays a minor The number of known examples of coordinated hunting, role in understanding the statistical properties of basket- whether intra or interspecies, has increased in the last ye- ball scoring. ars and examples include hawks, crocodiles, spiders, etc. As intriguing applications of this random-walk picture, Although there are some additional costs, hunting or de- we show that: (i) the distribution of times when the last fending in group may bring several benefits for predators lead change occurs, (ii) the distribution of times when and preys, respectively, what have been widely studied. the score difference is maximal, and (iii) the distribution Despite these mounting evidences, much less attention for the fraction of game time that one team is leading are has been dedicated to model such behavior. all given by the celebrated arcsine law–a beautiful and This problem has been recently considered within a game surprising property of random walks. We also use the theoretical framework in which the abundances of preys random-walk picture to construct the criterion for when and predators were assumed constant and only the frac- a lead of a specified size is ”safe”as a function of the time tion of those populations using either an individual or remaining in the game. This prediction generally agrees collective strategy evolves. Lett et al (2004 Theor. Pop. with comprehensive data on more than 1.25 million sco- Biol. 65 263) considered a mean field approach in which ring events in roughly 40,000 games across four profes- these densities are described by Lotka-Volterra-like equa- sional or semiprofessional team sports, and are more ac- tions, taking into account some of the advantages and curate than popular heuristics that are currently used in disadvantages for both preys and predators choosing a sports analytics. grouping strategy. More specifically, it is assumed that grouping lowers the risk of predation at the cost of incre- [02/11/2015 - 09:30 - Room Vit´oria] asing the competition for resources, while predators have CLUSTER APPROACH TO GELS AND a greater probability of success at the expense of having GLASSES, A. Coniglio, CNR-SPIN, Department of to share the prey with others. Physics, University of Naples “Federico II”, Via Cinthia, We present a spatial version of this model that locates 80126 Napoli, Italy A percolation theory is presented individuals or groups on a lattice and study it in the li- to describe the dynamics of the sol-gel transition[1]. mits of both low and high population viscosity (with or The same approach at mean field level is shown to without diffusion, respectively), and compare these re- describe also the dynamical critical behavior predicted sults with the mean field predictions. Of particular inte- by mode coupling theory (MCT) for the continuous glass rest is the coexistence region with both grouped and indi- transition[2]. A similar approach is extended to MCT vidual predators and prey persist within the population. for the discontinuous glass transition, more appropriate When compared with the mean field case, fundamental to describe the standard molecular glass transition. It is differences appear and are strongly affected by finite size shown that the relevant model now is given by Bootstrap effects. Percolation. This approach will provide a geometrical and physical interpretation of the critical exponents, [02/11/2015 - P002] Density classification performance of the elucidating the scaling laws and the universal aspect of Gacs-Kurdyumov-Levin four-states cellular MCT. automaton model IV and related automata, J. Ricardo G. Mendonc¸a, Rolf E. O. Simoes,˜ [1] A. Fierro, T. Abete and A. Coniglio, J. Chem. Phys. EACH/USP Almost four decades ago (in 1978), Gacs, 131, (2009) 194906. Kurdyumov, and Levin (GKL) introduced three different [2] J..J. Arenzon, A. Coniglio A. Fierro, and M. Sellitto. cellular automata (CA), which they called models II, Phys. Rev. 90, (2014) 020301(R) IV, and VI, to investigate whether nonequilibrium interacting particle systems are capable of displaying phase transitions. Their objective was to examine the “positive probabilities conjecture,” according to which one-dimensional particle systems with short-range inte- ractions and positive transition probabilities are always ergodic. This conjecture has been disproved—much 2 Abstracts - ENFE - 02/11/2015 to the awe of the practicing community—many times is to make a change of variable t for a new log-time since then, with the introduction of several models that τ ln(tc t), then to study the power spectrum of the have become archetypal models in theoretical computer new≡ series− thus generated. A consequence of this method science and nonequilibrium statistical mechanics. is the non-uniformity of the sample data, i.e. unequal spa- As a by-product of their investigations, GKL introduced cing between data points. The FFT-based techniques are the density classification problem in the cellular auto- not applicable, but one solution is to use the Lomb perio- mata literature. The density classification task consists dogram of Scargle, which is suitable for unevenly sampled in classifying arrays of symbols according to their ini- points. tial density of symbols using local rules, and is comple- We applying this method to study the Brazilian finan- ted successfully if a correct verdict as to which was the cial market, with the aim of detecting discrete scale inva- initial majority state is obtained in time at most linear riance in the Bovespa (Bolsa de Valores de S˜ao Paulo) in the size of the input array. Density classification is stock market index. Some historical price series have a nontrivial task for CA in which cells interact over fi- been selected for the periods in 1999, 2001 and 2008. We nite neighbourhoods, because then the cells have to achi- report evidence of detection of possible log-periodicity eve a global consensus cooperating locally only. Ultima- before breaks. tely, that means that information should flow through [02/11/2015 - P004] the entire system, be processed by the cells, and be not Study of String-Like Excitations in Artificial Spin destroyed or become incoherent in the process—entropy Ice Through Linear Chains of Magnetic Dipoles, must loose to work in the task, a relevant property in the Denis da Mata Oliveira, Lucas Alvares da Silva theoretical analysis of data processing and storage under Mol,´ UFMG In this work we study a system of classi- noise. For one-dimensional locally interacting systems of cal magnetic Ising-like dipoles (spins) on the square lat- autonomous and memoryless cells, emergence of collec- tice interacting exclusively by magnetic dipolar interac- tive behavior is required in these cases. In this context, tion. The spins are positioned in sequence, forming linear GKL model II has been extensively scrutinized as a model chains on the lattices links (strings), parallel to its plane. system related with the concepts of emergence, communi- These strings are constructed by Self-Avoiding Walks cation, efficiency, and connectivity. The other two GKL (SAW’s). For each string length, we generate all reacha- models, however, did not receive much attention. ble microstate the system in order to make the accurate Here we characterize the density classification perfor- analysis of the system through the canonical ensemble. mance of Gacs, Kurdyumov, and Levin’s “model IV,” a Our purpose is to better understand the behavior of ex- four-states cellular automaton with three absorbing sta- citation type strings in artificial spin ice, since the studied tes, by Monte Carlo simulations. We show that model IV chains are similar to those excitations, which present as compares well with its sibling model II in the density clas- remarkable feature the presence of magnetic monopole sification task, the additional states being barely relevant quasi-particles associated to its extremes. The studied for its performance. We also investigate the performance system is similar to the homopolymers in network-based of model IV under the influence of noise and show that SAW’s. We observed signs of phase transition which are it cannot perform the density classification task reliably analyzed by such quantities as the end-to-end distance at any nonzero level of noise, an indication that, most and the radius of gyration. Among the phase transi- probably, it becomes ergodic in this case. tion signs observed, the compact-extended transition of strings is well identified, such as the θ-transition of ho- [02/11/2015 - P003] Discrete Scale Invariance in Self-Organized Criti- mopolymers. We also present some results for the he- cality Systems, Andre´ Luis Brito Querino, Uni- xagonal lattice. We observed that the low temperature properties of the square lattice are determined by confi- versidade de S˜ao Paulo USP Recently, studies have shown evidence of log-periodic behavior in non- gurations that satisfy a rule of alternation in the direction hierarchical systems. A known case of log-periodicity of the walk, while in the hexagonal lattice the minimum or discrete scale invariance are systems that have a ge- end-to-end distance is the key factor. ometric hierarchy, for example the model Potts on the Thanks financial support: CNPq e FAPEMIG. diamond structure. The usual solutions of the renorma- [02/11/2015 - P005] lization group show that such systems have power laws Role of dimensionality of complex networks with 1 xb with complex exponents b C when near a critical metrics: Connection with nonextensive statistical point. An interesting fact is the∈ emergence of such pro- mechanics, S.G.A. Brito, L.R. da Silva, Depar- perties in real systems, for instance rupture and break- tamento de F´ısica Te´orica e Experimental - UFRN, C. down of complex materials and financial crashes. These Tsallis, Centro Brasileiro de Pesquisas F´ısicas - CBPF may be examples of complex systems with self-organized and Santa Fe Institute The study of networks is per- criticality (SOC). ceived in several fields of the science since many real sys- The detection of discrete scale invariance, or log- tems can be modeled as networks. The networks are periodicity in non-hierarchical systems presents nume- everywhere from social science to physics, biology, eco- rous difficulties. Parametric estimates using log-periodic nomics and other areas. Over the last decade a large functions can be flawed due to large fluctuations in va- number of empirical studies has been identifying pecu- lues, beyond the problem of degeneracy and multiple lo- liar properties in very different networks, for instance, cal minima in parametric regression estimation. For these the Internet and Online Social Networks (e.g. Facebook) reasons most research focuses on the use of nonparame- to citation networks and networks of neurons. Over the tric methods in detecting discrete scale invariance. A past years the concept of nonextensive statistical mecha- method widely used for the study of log-periodic data nics has been extremely successful in applications of the ENFE - 02/11/2015 - Abstracts 3 complex systems, in particular, on the complex networks. spatial scenarios in models under certain special condi- Deep connections are known to exist between scale-free tions. We follow another route, where the social inte- networks and non-Gibbsian statistics. For example, the ractions between any two agents is given by the descent typical degree distributions at the thermodynamical limit along the gradient of a cost function deduced from a Baye- k/κ are of the form P (k) eq− , where the q-exponential sian learning formalism. The cost functions depends on z ∝ 1 a hyperparameter that estimates the trust of one agent form e [1 + (1 q)z] 1 q optimizes the nonadditive q ≡ − − on the information provided by the other. If the expec- entropy Sq, basis of nonextensive statistical mechanics (which recovers Boltzmann-Gibbs statistical mechanics ted value of the total cost function is relevant informa- at the q 1 limit). We introduce and study here d- tion, Maximum Entropy permits characterizing the state dimensional→ growing networks with preferential attach- of the society. Furthermore we introduce a dynamics on αA the trust parameters, which increases when agents concur ment involving Euclidean distances through r− (αA ij ≥ and decreases otherwise. We study the resulting phase di- 0). Reinforcing the connection with q-statistics we nume- agram in the case of large number of interacting agents rically verify that the q-exponential degree distributions on a complete social graph, hence under sympatric con- exhibit, for both q and κ, universal dependences on the ditions. Simulations show that there is evolution of as- ratio α /d. Moreover, the q = 1 limit is exponentially A sortative distrust in rich cultural environments measured achieved by increasing α /d to infinity. A by the diversity of the set of issues under discussions. [02/11/2015 - P006] High distrust leads to antilearning which leads to mul- Contribution to the study of complex networks: tiple groups which hold different opinions on the set of Affinity model with metrics, S.G.A. Brito, L.R. issues. We simulate conditions of political pressure and da Silva, Departamento de F´ısica Te´orica e Experimen- interaction that describe the House of Congress of Bra- tal - UFRN Currently the interest in large-scale sys- zil and are able to qualitatively replicate voting patterns tems with a high degree of complexity has been much dis- through four presidential cycles during the years of 1994 cussed in the scientific community in various areas of kno- to 2010. wledge. As an example, the Internet, protein interaction, collaboration of film actors, among others. To better un- [02/11/2015 - P008] Spectral analysis of complex network with tuna- derstand the behavior of interconnected systems, several ble degree distribution and clustering, Roberta models in the area of complex networks have been propo- Pires Lins Machado, Josue´ Xavier de Carvalho, sed. Barab´asi and Albert proposed a model in which the Universidade Federal do Rio de Janeiro - Campus Xer´em connection between the constituents of the system could Network models provide a natural way to describe real dynamically and which favors older sites, reproducing a data set in diverse fields as for example biology, sociology, characteristic behavior in some real systems: connecti- ecology, internet, global economy and many others. Real vity distribution of scale invariant. However, this model networks display topological features, such as heavy tail neglects two factors, among others, observed in real sys- degree distribution, high clustering coefficient and assor- tems: homophily and metrics. Given the importance of tativity or disassortativity, that can not be modeled by a these two terms in the global behavior of networks, we totally regular or random graphs. Many random network propose in this dissertation study a dynamic model of pre- models has been proposed with the aim to capture featu- ferential binding to three essential factors that are respon- res regularly found in empirical networks. Two of most sible for competition for links: (i) connectivity (the more studied classes of such models are the Scale-free network connected sites are privileged in the choice of links) (ii) model (SF) and the Small-world network model (SW). homophily (similar connections between sites are more The original SF model display power law degree distri- attractive), (iii) metric (the link is favored by the proxi- bution and vanishing clustering coefficient in the limit mity of the sites). Within this proposal, we analyze the of large networks. The SW model interpolates between behavior of the distribution of connectivity and dynamic regular and random network by means a parameter p, evolution of the network are affected by the metric by that characterizes the disorder in the system. For small α parameter that controls the importance of distance A values of p the model exhibits high clustering coefficient. in the preferential binding) and homophily by η (cha- Since each graph/network can be represented by a matrix racteristic intrinsic site). We realized that the increased (adjacency, laplacian or normalized laplacian) its spectra importance as the distance in the preferred connection, (eigenvalues and eigenvectors) encodes information about the connections between sites and become local connec- its topology. The fluctuations of networks spectra can tivity distribution is characterized by a typical range. In be analysed through the framework of Random Matrices parallel, we adjust the curves of connectivity distribution, k/ηq Theory (RMT). Results for SW model have revealed that for different values of αA, the equation P (k)= P0e−q as the disorder level increases (decreasing clustering co- from the statistical non-extensive Tsallis. efficient) the spectral fluctuations follows the description [02/11/2015 - P007] of a transition between two different ensemble of RMT Sympatric multiculturalism or how distrust po- (Poisson-GOE transition). The main focus of this work larizes societies of Bayesian agents into groups, is to clarify the effects of clustering on the spectral pro- Felippe Alves, Nestor Caticha, IF Universidade de perties of complex networks. The clustering quantifies S˜ao Paulo While social interactions tend to decrease the likelihood that two neighbors of a vertex are neigh- differences in opinions, multiplicity of groups and indivi- bors themselves and it is related to number of triangles in dual opinion differences persist in human societies. Axel- the network (the vast majority of real networks display a rod identified homophily and social conformity seeking high density of triangles). The spectra of totally random as basic interactions that can lead to multiculturalism in networks (no correlation between vertices) is well descri- 4 Abstracts - ENFE - 02/11/2015 bed by GOE ensemble of RMT. According to RMT it is lit`ecnica de Catalunya The heterogeneous topology of expected that network models with high clustering (cor- a complex network can have a very relevant impact on the relation) deviates from GOE prediction. Here we analyse properties of dynamical systems running on top of it. Al- under the RMT framework the spectral properties of dif- ready classical studies in network science have thus shown ferent classes of random networks where both the degree- that a heterogeneous connectivity pattern can lead to a dependent clustering coefficient and the degree distribu- null percolation threshold, set a strong resilience against tion are tunable with special attention to the Scale-free random failures, as well as to induce a vanishing epide- networks. mic threshold for disease propagation. Similar and ad- [02/11/2015 - P009] ditional remarkable effects have been observed in a wide A resource based model for the competition of two variety of dynamical processes. Such dynamical effects, species with different metabolic pathways, Andre´ originally reported for static networks, in which nodes Amado†, Jorge Lenin†, Weini Huang‡, Paulo R. and edges are fixed and do not change over time, can A. Campos†, Fernando F. Ferreira§, †Evolutionary take a different, more complex turn when one considers Dynamics Lab, Department of Physics, Federal Univer- the intrinsic time-varying, temporal nature of many real sity of Pernambuco, 50670-901, Recife-PE, Brazil networks. Indeed, networked systems are often not static, ‡Department Evolutionary Theory, Max Planck Institute but show connections which appear and disappear with for Evolutio The evolution of cooperation is one of some characteristic time scales. Social networks represent the most intriguing conundrum in evolutionary biology. the prototypical example of this behavior, being defined Natural selection favours traits that increase individuals in terms of a sequence of social contacts that are continu- ability to reproduce or survive. One individual is said to ously established and broken. This mixing of time scales be cooperative if it provides a benefit to another indivi- can induce new phenomenology on dynamics on temporal dual or to a group at the expense of its own relative fit- networks, in stark contrast with what is observed in static ness. According to evolutionary theory, such cooperative networks. Moreover, the bursty nature of the time evolu- behaviour is detrimental to the individual and should be tion of temporal network contacts, characterized by long counter-selected. Instead, the selfish behaviour is expec- stretches of inactivity, interspersed by bursts of intense ted to rise and invade in a group of cooperators. However, activity, can complicate the picture, inducing for example cooperation is observed everywhere. Despite the greater a dynamical slowing down in dynamical processes as va- attention the issue has received in the last decades our ried as epidemic spreading, diffusion or synchronization. understanding about the underlying mechanisms that ex- The random walk is one of the simplest dynamical pro- plain the emergence and maintenance of cooperation is cesses, although still underlying many practical realistic still poor. Here we combine the frameworks of resource- applications such as diffusion, searching, community de- based modeling and evolutionary game theory to study tection and spreading dynamics. Even in this simplest of the conditions under which cooperative strains can thrive cases, a time-varying substrate can induce very noticea- in the context of metabolic pathways. ble differences with respect to the behavior expected in Heterotrophic organisms can produce adenosine static networks. In this work, we investigate the dynamic triphosphate (ATP) from two different mechanisms: relaxation of random walks on temporal networks by fo- inefficient (high rate) in which high amount of resources cusing in the recently proposed activity driven model [N. is consumed or efficient (high yield) with less resource Perra, et. al. Sci.Rep. 2, 469 (2012)]. For realistic acti- consumption. Despite cells with a higher rate but lower vity distributions with a power-law form, we observe the yield of ATP production may have a selective advantage presence of a very slow relaxation dynamics compatible compared with high yield cells, the latter is abundant with aging effects. A theoretical description of this pro- in nature. To be well succeeded the high yield cells cesses in achieved by means of a mapping to Bouchauds behave as a cooperative organism while high rate cells trap model. The mapping highlights the profound diffe- are non-cooperative. This problem can be posed in terms rence in the dynamics of the random walks according to of competition between cooperative and selfish cells. the value of the exponent γ in the activity distribution. Here we investigate the probability that a cooperative Acknowledgements and Financial Support: FAPEMIG, trait can invade a population of defectors (high rate of CAPES and CNPq. consumption, but low yield of conversion of research into [02/11/2015 - P011] energy) by using multilevel selection, where individuals Nonequilibrium phase transitions in a model of (cells) are organised into groups. Selection occurs at tax evasion dynamics, Rafael Mynssem Brum, both individual and group levels. Selection at the group Nuno Crokidakis, Instituto de F´ısica, Universidade level stems from differentiated growth rates that vary Federal Fluminense, Niter´oi/RJ, Brazil The agent- according to the group’s composition. We show in what based modeling of social interactions is one of the most conditions high yield cells can invade a population of interesting problems in statistical physics and it has been high rate metabolic cells. One of the motivation to extensively explored due to its high degree of applicability study this problem is understanding the transition of (opinion dynamics, spreading of diseases, systems with unicellular organism (high rate) to multicellular (high competing dynamics, rumor spreading, etc). Its results yield). are also very important from the point of view of an orga- [02/11/2015 - P010] nized society (stationary behavior of a disease, epidemic Slow relaxation dynamics and aging in ran- cycles, predominance of a certain political position and dom walks on activity driven temporal networks, many others). Angelica´ Sousa da Mata, Universidade Federal de In this work, we study the tax evasion dynamics in a Lavras, Romualdo Pastor-Satorras, Universitat Po- artificial society. In this case, we consider three types ENFE - 02/11/2015 - Abstracts 5 of individuals in relation to tax compliance: honest, tax change. In complex systems bifurcations can be associ- evaders and undecided individuals. We analyze the social ated with phase transitions in many different contexts: interactions that may occur among these individuals and from ecosystems to models of information traffic on the autonomous decisions, and we consider our population Internet, and from epidemic to social models. as a fully-connected network of N nodes, that characte- But if a complex system undergoes a phase transition rize a mean-field-like approach. Through analytical and when the control parameter reaches some critical value, numerical results, we found that even if initially there which are the critical aspects (in the sense of thermodyna- are no tax evaders in the population, these individuals mics / statistical mechanics) of such a transition? Trying may emerge in the society, and its fraction may stabilize to answer this question, it is presented an exploratory in the population, depending on certain socio-economic study, based on Thermodynamic Geometry (TG), of the parameters (efficiency of the influence of certain individu- critical aspects of some common bifurcations (e. g. trans- als over others, autonomous decisions, punishment rules, critical, saddle-node, etc.). etc). We show that this emergency of tax evaders is as- Thermodynamic Geometry (TG), also called Information sociated with a nonequilibrium phase transition. Geometry, is an approach based on Riemannian geome- try for the study of thermodynamical systems in equili- COMPLEXITY & INTERDISCIPLI- brium. The main idea is that the space of equilibrium NARY TOPICS states of a system is described by a metric, which is pro- portional to the Hessian matrix of entropy (or another [02/11/2015 - P012] thermodynamic potential) of the system with respect to On a connection between a class of q-deformed al- the thermodynamic parameters. In this context, the no- gebras and the Hausdorff derivative in a medium tion of “distance” between two states is associated with with fractal metric, Jose´ Weberszpil, Universidade probability of fluctuation between them: the less likely Federal Rural do Rio de Janeiro, UFRRJ-IM/DTL, Av. the fluctuation between the states, more distant they are. Governador Roberto Silveira s/n, Nova Igua¸c´u, Rio de From the metric one can obtain the corresponding geo- Janeiro, RJ, Brazil, Matheus Jatkoske Lazo, Uni- desic equations and curvature scalar R, which in turn is d versidade Federal do Rio Grande-FURG, Instituto de Es- proportional to the correlation volume ξ and therefore tat´ıstica e F´ısica, Rio Grande, RS, Brazil, Jose´ Ab- closely related to phase transitions. dalla Helayel-Neto,¨ Centro Brasileiro de Pesquisas In order to study the critical aspects of bifurcations in F´ısicas-CBPF-Rua Dr Xavier Sigaud 150, 2290-180, Rio such a geometrical approach, the normal form (described de Janeiro, RJ, Brazil Over the recent decades, di- by a nonlinear first-order differential equation) of each bi- verse formalisms have emerged that are adopted to ap- furcation is regarded (after some appropriate transforma- proach complex systems. Amongst those, we may quote tions) as a geodesic equation of some hypothetical model. the q-calculus in Tsallis version of Non-Extensive Sta- With this we obtain the corresponding metric, curvature tistics with its undeniable success whenever applied to scalar and also the exponent α, analogous to the critical a wide class of different systems; Kaniadakis approach, exponent of specific heat in a typical thermodynamical based on the compatibility between relativity and ther- system. modynamics; Fractional Calculus (FC), that deals with [02/11/2015 - P014] the dynamics of anomalous transport and other natural Mean-field approximation for the Sznajd model phenomena, and also some local versions of FC that claim in complex networks, Maycon S. Araujo,´ Andre´ to be able to study fractal and multifractal spaces and to M. Timpanaro, Carmen P. C. do Prado, USP - SP describe dynamics in these spaces by means of fractio- - Brazil, Fabio´ S. Vannucchi, UNESP - SP - Brazil nal differential equations. The question we might ask We will present a work in which we revisited mean fied is whether or not there are common aspects that con- approximations to the Sznajd model for opinion forma- nect these alternative approaches. In this short com- tion in a population connected through a general network. munication, we discuss a possible relationship between A master equation describing the time evolution of opi- q-deformed algebras in two different contexts of Statisti- nions is presented and solved in a hybrid approach. Ba- cal Mechanics, namely, Tsallis framework and Kaniadakis sed in a mean-field approximation, we were able to in- scenario, with local form of fractional derivative operators clude some features of the underling structure of the defined in fractal media, the so-called Hausdorff deriva- network though the estimate of some network parame- tives, mapped into a continuous medium with a fractal ters. Although quite simple, this approximation allows us measure. This connection opens up new perspectives for to capture the most important features regarding the ste- theories that satisfactorily describe the dynamics for the ady states of the model. When spontaneous opinion chan- transport in media with fractal metrics, such as porous ges are included, a discontinuous transition from consen- or granular media. Possible connections with other alter- sus to polarization can be found as the rate of sponta- native definitions of FC are also contemplated. Insights neous change is increased. The main point in this work on complexity connected to concepts like coarse-grained is the presentation of a hybrid mean-field approach, that space-time and physics in general are pointed out. includes interactions between second nearest neighbors, [02/11/2015 - P013] that are necessary to estimate correctly the critical point Critical aspects of bifurcations in a geometrical of the transition. The analytical prediction of the cri- approach, A. Mihara, Universidade Federal de S. tical point is also compared with numerical simulations Paulo Bifurcation is a qualitative change that occurs in a wide variety of networks, in particular Barab´asi- in the behavior of a dynamic system when the value Albert networks, finding reasonable agreement despite of a (control) parameter undergoes a small and smooth the strong approximations involved. The same hybrid 6 Abstracts - ENFE - 02/11/2015 approach that made it possible to deal with second-order tions, several quantities characterizing the mixed-games neighbors could just as well be adapted to treat other pro- are still the same as the ones obtained in the average blems such as epidemic spreading or predator-prey sys- game when the two games are not very different. Also we tems. The work has been published in PRE 91, 022813 find interesting results regarding how the heterogeneity (2015). of the games played can increase the final fraction of coo- perators above the usual mean game limit. We would like [02/11/2015 - P015] Interspike intervals in neuronal networks with to thanks FAPEMIG and CNPq for the financial support self-organized criticality, Osame Kinouchi, Ari- given. adne A. Costa, Lezio´ Bueno, Geraldine Bosco, [02/11/2015 - P017] Universidade de S˜ao Paulo, Mauro Copelli, Universi- Differences between quenched and annealed dade Federal do Pernambuco The distribution of inters- neuronal networks with self-organized critica- pike intervals (ISI) for individual neurons is a standard lity, Ariadne de Andrade Costa, Osame Kinou- measure in neuroscience. Recently the idea that neuro- chi, FFCLRP - USP - SP - Brasil, Joao˜ G. F. Cam- nal networks works near the critical point of an absorbing pos, Mauro Copelli, UFPE - PE - Brasil In a re- state transition has been explored, with a lot of theoreti- cent work, mean field analysis and computer simulations cal and experimental results. Surprisingly, it appears that were employed to analyze critical self-organization in an the distribution of ISI for neurons pertaining to a criti- annealed network of excitable cellular automata (SIRS) cal network has not been measured yet. Here, we model neuronal networks, where randomly chosen synapses are the generation of interspike intervals in neuronal systems depressed after each neuron spike. Calculations agree with self-organized criticality by two methods. First, we with simulations of the annealed version, showing that use a simple stochastic model where interavalanches in- the nominal branching ratio (σ) converges to the critical tervals (IAI) are generated from a power law distribution, value σc = 1 and fluctuations vanish in the thermodyna- and interspike intervals are sums of IAI between two times mic limit, as expected of a self-organized critical system. where a neuron spikes because it pertains to an avalan- However, the question remains whether the same results che (modeled by a neuronal size avalanche distribution occurs to the quenched version of the model (which is bi- 3/2 P (S)= cS− ). Second, we perform a full simulation in ologically more plausible) where neighborhoods are fixed a network of excitable elements (neurons) with dynamical and only the acting synapses are depressed. We have seen synapses which presents well behaved self-organized criti- that simulations of the quenched model yield a stationary cality. In this simulation we define avalanches as sequen- value σ(t )=1.105 which is a significant deviation tial activity above some threshold level of active sites. from σ =→ 1, ∞ due to spatio-temporal correlations produ- This enables us to define interavalanches (IAI) and inters- ced by avalanches. However, the model is shown to be pikes intervals (ISI) and construct histograms for them. critical, as the largest eigenvalue λ of the synaptic matrix We compare these results with ISI distributions that pre- is shown to approach λc = 1, with fluctuations vanishing sent power law tails from real neurons from cortical and in the thermodynamic limit. We also study the influence thalamic areas of freely behaving rats. We find that self- of the recovery and decay synaptic parameters in both organized criticality can explain power laws in the tail of types of models, as well the influence of the number of ISI distributions of neurons and that such power laws in neighbors. As a future work, we intend to study the dis- single neurons could suggest the presence of criticality at tribution of interspike intervals in this kind of neuronal the network level. networks. [02/11/2015 - P016] [02/11/2015 - P018] Cooperation in two-dimensional mixed-games, Real genomic representation for topopatric speci- Marco Antonio Amaral, Jafferson Kamphorst ation, Camilo Rodrigues Neto, Sergio Candido Leal da Silva, Lucas Wardill, Universidade Fe- de Oliveira Junior, EACH, University of Sao Paulo, deral de Minas Gerais - MG - Brasil, University of Brasil The evolutionary theory for speciation has pro- British Columbia - Vancouver Vancouver, BC, Canada duced several models to explain the diversity of life. Four Evolutionary game theory is a common mathematical of these models are called alopatric, peripatric, parapa- framework to study the evolution of cooperation in sel- tric and sympatric, named after the predominant kind fish systems, specially using the Prisoners Dilemma game, of genetic flux disturbance acting among the population. where it is usually assumed that the same game is played Physical barriers and ecological interactions are the two in all interactions. Here, we investigate a model where usual factors. the game that is played by two individuals is uniformly Recently, a new kind of speciation was proposed [Aguiar, drawn from a sample of two different games at each itera- 2009]. The model is based in selective mating determined tion. Using the master equation approach we show that by genetic affinity and spatial proximity. The model does the random mixture of two games is equivalent to play not include any kind of geographical barrier, ecological the average game when (i) the strategies are statistically interaction or natural selection. This new kind of specia- independent of the game distribution and (ii) the transi- tion is named topopatric, as it emphasize the role of the tion rates are linear functions of the payoffs. This result spatial auto organization of the species origin and distri- still holds using Pair-Approximation for a small cluster of bution. The usual approach is to represent a specimen as 8 sites arranged in a square lattice. We also use Monte- a binary string, and defying usual genetic algorithm with Carlo simulations in a two dimensional lattice to investi- crossover and mutation operators, but without the se- gate the scenario when the two above conditions do not lection phase. Differences between individuals reproduc- hold, i.e. we use the Fermi-Dirac distribution for the tive rates are randomly attributed at reproduction time, transition rates. We find that even outside of such condi- rather than being due to any special ability. ENFE - 02/11/2015 - Abstracts 7
We present a modified model with a real-coded genetic N, residual correlations will eventually drive it to ther- algorithm, where the specimens are represented by real modynamic equilibrium (if such equilibrium exists, which numbers, with modified crossover and mutation opera- is not the case for 3d gravitational systems) after a time tors. The reproductive rates are similar to the previous t which scales with N as t N δ, where δ is a system binary model. We studied the number of new species, specific× exponent. On the other× ∼ hand, in the thermody- the abundance of the species, and the distribution of the namic limit, N , the system will remain trapped in species over space in function of the genomic distance to- a stationary state→ forever.∞ In this collisionless limit, the lerance, the searchable radius for mating, the migration relaxation to stationarity is a result of Landau damping, and the random reproductive rate. All simulations start which transfers the energy of collective oscillations to the with an initially uniform population, with the same real individual particles. Once the oscillations of the mean- coded genetic algorithm crossover and mutation opera- field potential die out, the particles will move in a static tors, and were performed in torus and ring like spatial mean-field potential. If a system has sufficient symmetry, setups. The results resemble the ones already reported the motion of particles in a static potential will be inte- in the literature, but are not strictly equivalent. We at- grable, and the ergodicity will be irrevocably broken. In tribute this difference mainly to topological differences in this paper we will explore the role of chaotic dynamics on the space of the genome representation. the time that a system with LR interactions remains trap- Aguiar, M. A. M., et al. ”Global patterns of speciation ped in a QSS. We discover that a small degree of chaos, and diversity.”Nature 460.7253 (2009): 384-387. measured by the Lyapunov exponents, favors a faster re- laxation to equilibrium. Surprisingly, a larger degree of [02/11/2015 - P019] Long-range correlations of the wind velocity chaos hinders the relaxation to equilibrium. in Salvador-BA, Jose´ Vicente Cardoso Santos, [02/11/2015 - P021] Davidson Martins Moreira, Marcelo A. Moret, Entropic simulations of the spin-1/2 Baxter-Wu SENAI CIMATEC The use of wind energy has been model., Lucas Nunes Jorge, Alvaro´ de Al- increasingly adopted worldwide. In Brazil, the Northeast meida Caparica, Universidade Federal de Goi´as - UFG has been a strong investment option in the industry be- Among the various models used to describe spins sys- cause recent wind maps have been made and show high tems the Baxter-Wu model is particularly interesting, attendance rates winds around its coastline and in some since it considers triplets of spins, thus, it does not pre- cases inside. Notwithstanding this great wind potential sents spin-reversal symmetry, as it occurs in the most scenario, Brazil has this alternative source of energy with know models. This model is defined in a triangular two- a low representation in its energy mix. Thus, in order dimensional lattice, and the three-spin interaction is gi- to justify an increase in interest and reliability in this ven by the Hamiltonian, alternative source of energy, is presented in this paper a preliminary analysis of time series representing the in- H = J s s s , (1) BW − i j k tensity records and wind direction in fixed unit of data collection on the drive SENAI / CIMATEC in Salvador, X Bahia. This analysis is done by Destrended Fluctution where the variables of spin are located at the vertices of Analysis method (DFA). Using this method substantia- the triangular lattice and take the values si = 1, J is tes the viability of wind farms installations in the areas the coupling constant that defines the energy scale± and of collection, because it indicates the possibility of long- the sum extends over all the triangular faces. For the range correlation in the distribution of magnitudes analy- spin-1/2 case, the model was exactly solved by Baxter zed what may prove the function of the constancy of the and Wu, and presents the same critical temperature of wind flow, and, with this, enable the facilities of wind tur- the Ising model, but the critical exponents are those of bines more efficiently and effectively. Preliminary results the q = 4 Potts model, so, this model is an excellent ob- indicate that the number of local data is persistent in ject of study to test new Monte Carlo procedures. Monte direction, speed and related thermodynamic quantities, Carlo simulations are an efficient tool to calculate critical which corroborates the feasibility of wind participation temperatures and static critical exponents. In particu- in local energy matrix. lar, the Wang-Landau sampling has become in last years Keywords: Renewable Sources. DFA. Complexity. Wind more and more accurate and robust. In this work we pre- Energy. sent a simulational study of the pure spin-1/2 Baxter-Wu model using a modified Wang-Landau scheme to calcu- FUNDAMENTAL ASPECTS OF STA- late the critical exponents γ, β and ν and the critical TISTICAL MECHANICS temperature Tc in the Baxter-Wu model. In this new procedure, instead of updating the density of states af- [02/11/2015 - P020] ter every spin-flip we adopt the Monte Carlo sweep for Chaos and relaxation to equilibrium in systems updating the density of states, the microcanonical avera- with long-range interactions, Felipe L. Antunes, ges are accumulated only after a few Wang-Landau levels Fernanda P. C. Benetti, Renato Pakter, Yan Le- have already run out, and stop the simulations when a vin, Universidade Federal do Rio Grande do Sul In the checking parameter, ε, which measures the fluctuation of thermodynamic limit, systems with long-range (LR) in- the peak of the specific heat during the simulations, va- 4 teractions do not relax to equilibrium, but become trap- ries below 10− for a complete Wang-Landau level. As a ped in non-equilibrium stationary states. Once a system result, different runs proceed up to different final modi- is trapped in a non-equilibrium state, two outcomes are fication factors. Moreover, the final results are obtained possible: if the system has a finite number of particles as averages over ten independent sets of finite size sca- 8 Abstracts - ENFE - 02/11/2015 ling simulations. Our results are very consistent and we this project given by FAPEMIG and CNPq. compare them with exact data available in literature. [02/11/2015 - P024] [02/11/2015 - P022] The Cluster Expansion in Statistical Mechanics: Mass segregation on Hamiltonian Mean Field mo- Holder inequality, Jose´ Andre´ Lourenc¸o, UFES del, J.R. Steiner, Zolacir T.O.Jr, Universidade In this review, the Glimm-Jaffe-Spencer cluster expan- Estadual de Santa Cruz, T.M. Rocha-Filho, Univer- sion from constructive quantum field theory is adapted sidade de Bras´ılia The dynamical evolution of young to treat quantum statistical mechanical systems of parti- stellar clusters is thought to be entangled with that of cles interacting by finite range potentials. The Hamilto- the mass segregation. This is a common sense in the as- nian H0 + V need be stable in the extended sense that trophysical community. In this work, mass segregation H0 +4V + KN ≧ 0 form some K. In this situation, phenomena (MSP) is investigated, as a dynamical fea- with a mild technical condition on the potentials, the ture, using the Hamiltonian Mean Field (HMF) model. cluster expansions converge and infinite volume limit of The study of MSP in the HMF model is justified by the the correlation functions existes, at low enough density. fact that stellar and galaxies clusters are clearly examples These infinite volume correlation functions cluster expo- of systems with long range interaction as the HMF itself nentially. Following the usual literature, we define a class and exhibits MSP. To achieve this aim, we introduce diffe- of interacting boson and fermion particle theories with a rent masses in the Hamiltonian and perform computatio- matter-like potential, 1/r suitably truncated at large dis- nal simulations for that HMF model. We focus in looking tance. This system would collapse in the absence of the for what happens over the mass distribution in the phase exclusion principle. The potential is unstable, but the space for the system. A comparison with a short range Hamiltonian is stable. This provides an example of a version of HMF, with only first neighbours interaction, system for which this is method proves existence of the that is known as XY-model, also with different masses infinite volume limit, that is not covered by the classic is made. The integration of the equations of motion is work of Ginibre, which requires stable potentials. The conducted using a fourth-order symplectic Omelyan in- main focus of this review is to discuss a key ingredient, tegrator. We analyse what happens through the violent a type of Holder inequality for the expectation values of relaxation period and what stand for the quasi-stationary spatially smeared Euclidian densities, a special interpo- states (QSS) of this dynamics. The results obtained sup- lation theorem. The cluster expansion as developed here port the fact that MSP is observed already in the violent is purely a geometric analysis of the paths that realize relaxation time and is maintained during the QSS’s that the traces in path space. The total path space integral come after that. Some structures are observed in the is split into subsets in which paths avoid certain regions mass distribution function. Another result of this study and must hit other regions. is that the mass distribution is determined by the system [02/11/2015 - P025] dynamics and is independent of the dimensionality of the Holographic considerations on non-gaussian sta- system. MSP occurs in a one dimensional system as a tistics and gravothermal catastrophe, Everton result of the long range forces that acts in the system. M. C. de Abreu, Universidade Federal Rural do Rio de [02/11/2015 - P023] Janeiro, Jorge Ananias Neto, Universidade Federal The BKT phase transition in the diluted XY Mo- de Juiz de Fora, Edesio´ M. Barboza Jr., Universidade del, Tatiana Pena Figueiredo, Julio Cesar´ Si- do Estado do Rio Grande do Norte, Rafael C. Nunes, queira Rocha, Bismarck Vaz da Costa, UFMG Universidade Autonoma de Barcelona The mechanism The Berezinskii - Kosterliz - Thouless phase transition of gravothermal instability, discovered by Antonov is an (BKT) in the diluted XY Model is studied in detail with important phenomena in gravitational thermodynamics. Monte Carlo simulations using the Wang-Landau algo- It has been very helpful for an extensive research con- rithm. The transition temperature, TBKT , was found cerning statistical mechanics of long range interactions by the finite size scaling of the helicity modulus and the systems in several fields in physics. This connection with in-plane magnetic susceptibility, for various magnetics si- thermodynamics and statistical mechanics has motiva- tes density. The spin-spin correlation function was cal- ted us to investigate statistically the gravothermal catas- culated, for all spins, and for those who are within the trophe cluster that percolates, in temperatures after and below At the same time, there are theoretical evidences that the TBKT , for various sizes lattices. We intend, using the understanding of gravity has been greatly benefited these results, discuss the mechanism of the unbinding of from a possible connection with thermodynamics. Pio- the vortices - antivortices pairs in the transition process. neering works of Bekenstein and Hawking have described The position of this pairs and the vortices and antivor- this issue. For example, quantities as area and mass of tices density was calculated for various magnetics sites black-holes are associated with entropy and temperature density. A well accepted theory suggests that this is the respectively. Working on this subject, Jacobson inter- main mechanism responsible for the transition. However, preted Einstein field equations as a thermodynamic iden- we observed that vortices do not move for long distances tity. Padmanabhan gave an interpretation of gravity as through the lattice and in the vacancies neighborhood it an equipartition theorem. can be pinned there for a long period of time, then if the In this paper we have derived the equipartition law of transition occurs even with pinned vortices, these can not energy using Tsallis formalism and the Kaniadakis power be responsible for the transition, and another mechanism law statistics in order to obtain a modified gravitational should be observed. There is another theory which sug- constant. We have applied this result in the gravothermal gests such polymerizations walls domain is responsible for collapse phenomenon. We have discussed the equivalence transition. We would like to thank the partial support of between Tsallis and the Kaniadakis statistics in the con- ENFE - 02/11/2015 - Abstracts 9 text of Verlinde’s entropic formalism. In the same way (MWLS) and, in the case susceptibility, the MWLS we have analyzed the negative heat capacities in the light shows that there is a limit to begin the accumulation of gravothermal catastrophe. The relative deviations of the microcanonicas media. This study also showed the the modified gravitational constants are derived. density of states for the Baxter-Wu Spin-1/2 and Spin- 1, and perform a scale analysis of finite size for the model. [02/11/2015 - P026] Classical Origns of Frequency Probabilities, [1] A. A. Caparica and A. G. Cunha-Netto, Phys. Rev. Guilherme Roncaratti Galanti, Osame Kinou- E 85, 046702. chi, Universidade de S˜ao Paulo A classical open pro- [2] Fugao Wang and D. P. Landau Phys. Rev. Lett. 86, blem in probability theory concerns the so called three 2050. sided dice: suppose a cylinder with diameter d and height h, what should be the ratio h/d so that the frequency to [02/11/2015 - P028] obtain a face is 1/3 (where falling on the cylinder late- Classical dynamics of two electric charges, ral side counts as a face). In a more general situation, Rodrigo R. Silva, Annibal Figuereido, Universi- we can ask for what is the frequency P (S d,h,H,θ,ε) for dade de Bras´ılia From Li´enard-Wiechert fields that obtaining a fall on the lateral side S given| the cylinder describe the classical electromagnetic effect of a moving diameter d and tallness h, the height H of its center of electric point charge, we construct the equations of mo- mass above a table at the moment of launching, the ini- tion for two particles. To build the equations of motion we use the Lorentz force dp = q(E + v B) where E and tial angle θ of the cylinder axis with the horizontal and dt × the elastic coefficient of restitution ε, which depends on B are the Li´enard-Wiechert fields and p is the relativistic the materials of the table and the cylinder. We do not momentum. Rescheduling to simplify this equations and consider here initial conditions with translational or ro- thus get only one parameter, the ratio of masses. As the tational velocities. We made experimental measures for system depending only on the ratio of the masses and the P (S d,h,H,ε) varying h/d with H large (so that influ- boundary conditions we can simulate a electron-positron ences| of initial conditions vanishes) and ε fixed. We also system, a electron-proton system and any two electric model numerically the system as a two dimensional “cy- point charge system varying the ratio of mass. For each linder”of height h and d composed by four masses linked system we can analyze the position and velocity of each by springs. Given initial conditions, the numerical result particle, the center of mass and the decay time. Analy- is deterministic. However, for H > 30cm, the toss out- zing the equations of motion in polar coordinates, and come depends strongly on initial conditions, so that we the results of the numerical solutions can be analytically must average over a cell ∆H.∆θ of initial conditions un- deduce the decay time of the particles depends on the der control of the experimenter. This average furnishes initial radius and the ratio of the mass. From the Lan- a frequency P (S d,h,H,ε) to be compared to the expe- gevin equation that is a stochastic differential equation rimental results.| By using the unknown ε as a free para- we search the noise term that generates the stability of meter, we obtain a very good agreement between the full existing orbits on the electron-proton model. three-dimensional experiment and the two-dimensional [02/11/2015 - P029] simulation. The experimental data can also be fitted by A study of Entropies and Non-linear Cons- a recently proposed “Gibbs curve”. traints in Long Range Interacting Systems, Moises F. Junior, Marco A. Amato, Universidade [02/11/2015 - P027] de Bras´ılia The statistical basis for entropy has been Efficient Wang-Landau Sampling for the Baxter- laid by Boltzmann and Planck giving S = k ln Ω, Wu Model, Maria Lucia´ M. Costa, Universidade N B where S is the total thermodynamic entropy of the sys- Federal do Par´a, Joao˜ Antonio Plascak, Univer- N tem, N the number of entities, Ω the statistical weight sidade Federal da Para´ıba This study analyzes the or number of possible realisations (e.g., microstates) of two-dimensional Baxter-Wu model of Spin-1/2 and the system, of equal probabilities, and k is the Boltz- Spin-1 through the revenue to improve the accuracy B mann constant. For a discrete system one may also write of the Wang-Landau sampling. This scheme of Wang- S = k p ln p . p is the probability of occurrence of Landau was proposed by Caparica and Cunha-Neto[1]. B i i i the i−-th distinguishable outcome or state, from a total of Remembering that Wang-Landau simulation primor- s such states.P In a previous paper those authors provide dial [2] generates the density of states g(E), i.e. the a natural extension of the Boltzmann counting method number of all possible states (or configurations) for any in order to obtain generalized entropies, which leads to a energy level E of the system, allowing determine the statistical interpretation based on the occupational sta- canonical average of any thermodynamic variable, as βE tistics of a stochastic process. In this paper we impose E qc, destroys the consensus found in epidemic spreading could be observed. However, ENFE - 02/11/2015 - Abstracts 11 many complex networks are, theoretically, infinite Figueiredo, and J. R. Steiner; Phys. Rev. E 89, 032116 dimensional systems, while GPs are expected for finite (2014). dimensions. It was recently observed Griffiths effects in the contact process (CP) on weighted heterogeneous [02/11/2015 - P035] networks in the form of slow relaxations. On the other FRACTAL RESERVOIRS, Tawan Tayron de An- hand, the fundamental susceptible-infected-susceptible drade Carvalho, Valdemiro da Paz Brito, Uni- (SIS) model, which have a null critical point for infinite versidade Federal do Piau´ı, Marcelo Andrade de Fil- networks, exhibit multiple transition involving localized gueiras Gomes, Universidade Federal de Pernambuco configuration when investigated in a single finite-size Since the 1980s, when the study of fractal structures network. In this work, we performed simulations of the quickly began to move forward, there are many unexplo- SIS model on a large ensemble of random networks red aspects about how a fractal object interacts with cer- N γ tain types of environments. In some of these problems ex- with power-law degree distributions P (k) k− , with and without a hard upper cutoff, the latter∼ rendering tensively studied in the area of out-of-equilibrium growth networks with outliers, which were pointed as the origin models, there is an interaction of a fractal surface with of multiple localized transitions. We observed that for external fields, especially a diffusive field [1, 2, 3]. In a fixed and small infection rate many samples were in another area, the study of the interaction of molecules the inactive phase, but some exhibited quasistationary with enzymes, proteins, or biological tissues is modeled states with very long and highly fluctuating lifespans, by fractals, in catalytic processes or adsorption, an area analogously to the rare region effect. In the absence of known as chemistry between 2 and 3 dimensions [4]. More cutoff, a logarithmic decay to the absorbing state was recently, Balankin and coworkers studied the kinetics of observed. For a hard cutoff we found power-law decays the water escape from aluminum crumpled surfaces [5]. with non-universal exponents in agreement with a GP. This process occurs until the mass of water absorbed by Even being observable on very large networks (N 108), the surface reaches the threshold of the water retained a finite size scaling shows that the GP-like∼ regime in the structure by surface tension forces. In the pre- disappears in the infinite size limit. We acknowledge the sent work, we are interested to study experimentally the financial support of CNPq and FAPEMIG. process of absorption and retention of water by crumpled wire balls, a much more effective system for the retention [02/11/2015 - P034] of water than the aluminum wrinkled surfaces. These An study of the scaling of the dynamics of homo- balls were made with wire with 1.5 mm in diameter and geneous states of one-dimensional gravitational had radii varying from 0.40 to 5.15 cm. Such structures systems, Lydiane F. Souza, M. A. Amato, T. M. are obtained in approximately spherical shapes by the Rocha Filho, UnB Quasi-Stationary States of long- use of nearly isotropic external compressing forces and range interacting systems have been studied at length have a high porosity and behave as fractals. The ex- over the last fifteen years. A pair interaction potential periment consists of (1) generation and characterization is said to be long ranged if it decays at long distances of the balls, (2) immersion of the samples in water by α as r− with α d where d is the spatial dimension. ≤ following a fixed protocol, (3) evaluation of the retained Kinetic equations for long-range interacting systems usu- water as a function of the size of the balls, the volume of ally can be obtained from the BBGKY hierarchy [1] by pores, and the surface roughness, and (4) comparison of taking into account contributions from the two-body cor- the experimental results with a simple mean-field model. relation functions, which are of order 1/N [2] that result [1] M. Plischke, Z. R´acz, Phys. Rev. Lett. 53, 415 in a time scale of collisional relaxation proportional to (1984). N. The Balescu-Lenard equation for a one-dimensional [2] A. Coniglio, H. E. Stanley, Phys. Rev. Lett. 52, 1068 homogeneous system vanishes identically due to the Di- (1984). rac delta function. Therefore higher order terms must be [3] M. A. F. Gomes, G. L. Vasconcelos, G. C. Nascimento, kept when truncating the hierarchy, leading to a different J. Phys. A: Math. Gen. 20, L1167 (1987). scaling of the time evolution of a homogeneous state. It [4] D. Avnir, D. Farin, J. Chem. Phys. 79, 3566 (1983). would be natural to expect that in the present case the [5] A. S. Balankin et alli, Phys. Rev. E 83, 036310 (2011). predominant collisional corrections to the kinetic equa- tion come from higher order terms proportional to 1/N 2, [02/11/2015 - P036] this implies a relaxation scaling proportional to N 2. Ina Analysis of the Shear Viscosity of Binary Gaseous previous report [3] it is shown that the scaling from the- Mixtures near from the Chemical Equilibrium, oretical considerations for the HMF and Ring Model is Adriano Willian da Silva, Kayk Bueno Martins, proportional to the square of the number of particles and At´ıria Sbrissia, Instituto Federal do Paran´a- Cˆampus have also, in the former case, confirmed by computatio- Curitiba A binary gaseous mixture with reversible re- nal calculations. In this report we propose an extension action of type A+A=B+B is studied with Boltzmann of the theoretical calculations given in Ref. [3] for a 1D equation, assuming hard spheres cross sections for elastic gravitational system in order to provide a kinetic equa- collisions and two models for reactive interactions: line- tion for such systems. of-centers model and modified line-of-centers model. The [1] R. L. Liboff, Kinetic Theory - Classical, Quantum, and Chapman-Enskog method is used to obtain the solution Relativistic Descriptions, 3rd ed, Springer-Verlag (New of the Boltzmann equation in a chemical regime for wich York, 2003). the reactive interactions are of the same order as the elas- [2] R. Balescu, Statistical Dynamics - Matter out of Equi- tic one, i.e. in the system is closed to the final stage of librium, Imperial College Press (London, 1997). a chemical reaction where the affinity is considered to be [3] T. M. Rocha Filho, M. A. Amato, A. E. Santana, A. a small quantity and the system tends to the chemical 12 Abstracts - ENFE - 02/11/2015 equilibrium. This kind of reaction is known as fast reac- Nowadays, new and elaborated experimental, technolo- tions. The internal degrees of freedom of the particles of gical , and industrial situations require new and advan- the gas are not taken into account. The value of the re- ced phisico-chemical theoretical formalisms. We consider action heat distort the Maxwellian distribution function, here one such case, which appears to provide a good il- for large values the effect becomes more importants. The lustration: the so called Therma laser Stereolithography. resulting integral equation is solved with the expansion of This is a recent technological process that allows solid the distribution function in Sonine polynomials. The aim physical parts to be made directly and rapidly from com- of this paper is to evaluate the influence of the chemical puter data. In a description of order one it is presented reactions on the shear viscosity coefficient of the mixture. an analysis of the conditions necessary for a satisfactory It was verified the reaction heat changes the shear visco- characterization of the technological process of thermal sity and these change differ for exotermic and endotermic prototyping. We also consider the nonequilibrium ther- reactions. The change is bigger for endotermics reactions modynamic aspects of the related techno-industrial pro- and reactive interactions of line-of-centers model. cess of thermal laser stereolithography. [02/11/2015 - P037] [02/11/2015 - P039] Anomalous temperature relaxation for polymeric Nonlinear Ehrenfest’s Urn Model, chains, Rogelma M. S. Ferreira, Universidade Gabriela A. Casas, Fernando D. Nobre, Evaldo Federal do Recˆoncavo da Bahia - UFRB, David L. M. F. Curado, Centro Brasileiro de Pesquisas F´ısicas Azevedo, Fernando A. Oliveira, Universidade de - CBPF The Ehrenfest’s urn model (sometimes also Bras´ılia - UnB We analyse the temperature relaxation called Ehrenfest’s flea model) has played an important phenomena of small polymeric chains in contact with a role in clarifying the foundations of statistical mecha- thermal reservoir. We simulate the chains in a fluid [1], nics, providing an interpretation of irreversibility in a and we show that they reproduce a behaviour predicted statistical manner. The model is defined by N balls by recent theoretical investigations [2]. The temperature distributed in two urns (or boxes) 1 and 2, such that decay reveals the existence of an anomalous cooling in at each discrete instant of time s, a ball is chosen at which the temperature may oscillate [2]. This effect is random and moved from the box in which it is found a consequence of collective behaviour of the monomers to the other box. At the beginning of the 20th century, in the chain, which builds up correlation [3,4]. This such a simple model was useful in explaining the heat anomalous behaviour however, does not violated the exchange between two bodies at unequal temperatures, second law of thermodynamics [5]. We analyse as well where the temperatures are mimicked by the number the scaling [4,6] dynamical properties of the chains. of balls in each box, and the heat exchange becomes a random process. In the present work the Ehrenfest’s [1] A. M. Maroja, F. A. Oliveira, M. Ciesla, and L. Longa, urn model is modified by introducing nonlinear terms in Phys. Rev. E 63, 061801 (2001). the associated transition probabilities. It is shown that [2] L. C. Lapas, R.M. S. Ferreira, J. M Rub´ı,and F. A. these modifications lead, in the continuous limit, to a Oliveira, J. Chem. Phys. 142, 104106 (2015). Fokker-Planck equation characterized by two competing [3] F. A. Oliveira, Phys. Rev. B 57, 10576 (1998). diffusion terms, namely, the usual linear one, as well [4] C. L. Dias, M. Dube, F. A. Oliveira, and M. Grant, as a nonlinear diffusion term, typical of anomalous Phys. Rev. E 72, 011918 (2005). diffusion. By considering a generalized H-theorem, the [5] L. C. Lapas, R. Morgado, M. H. Vaintein, J. M. Rub´ı, associated entropy is calculated, resulting in a sum of and F. A. Oliveira, Phys. Rev. Lett. 101, 230602 (2008). Boltzmann-Gibbs and Tsallis entropic forms. It is shown [6] R.M. S. Ferreira, M. V. S. Santos, C. C. Donato, J. that the stationary state of the associated Fokker-Planck S. Andrade, Jr. and F. A. Oliveira, Phys. Rev. E, 86, equation satisfies precisely the same equation obtained 0211211 (2012). by extremization of the entropy. Moreover, the effects of [02/11/2015 - P038] the nonlinear contributions on the entropy production Thermal Transport in a Higher-Order Generali- phenomenon are also analyzed. zed Hydrodynamics, Carlos A. B. Silva, Instituto [02/11/2015 - P040] Tecnol´ogico de Aeron´autica SP Brasil, Cloves´ G. Ro- Statistical Thermodynamics of the Fr¨ohlich-Bose- drigues, Pontif´ıcia Universidade Cat´olica de Goi´as De- Einstein Condensation of Non-Equilibrium Mag- partamento de F´ısica Goiˆania GO Brasil, J. Galvao˜ nons, Fabio S. Vannucchi, Campus do Litoral Pau- Ramos, R. Luzzi, Institute of Physics Gleb Wataghin lista, Universidade Estadual Paulista - UNESP, Ro- State University of Campinas 13083-859 Campinas SP berto Luzzi, Instituto de F´ısica “Gleb Wataghin”, Uni- Brasil Thermal transport in classical fluids is analy- versidade Estadual de Campinas - UNICAMP A non- zed in terms of a Higher-Order Generalized Hydrodyna- equilibrium statistical-thermodynamic approach to the mics (or Mesoscopic Hydro-Thermodynamics) , that is, study of a Fr¨ohlich-Bose-Einstein condensation of mag- depending on the evolution of the energy density and its nons under radio-frequency radiation pumping is presen- fluxes of all orders. Its derived in terms of a kinetic the- ted. Such a system displays a complex behavior consis- ory based on the Non-Equilibrium Ensemble Formalism. ting in steady-state conditions leading to the emergence The general system of coupled evolution equations is deri- of a synergetic dissipative structure resembling the Bose- ved. Maxwell times which are of large relevance to deter- Einstein condensation of systems in equilibrium, due to mine the character of the motion are derived. They also a peculiar and fundamental contribution of a non-linear have a quite important role for the choice for the contrac- character related with the magnon-lattice interaction. tion of description (limitation in the number of fluxes to A kind of “two fluid model”arises: the “normal”non- be retained) in the study of the hydrodynamic motion. equilibrium structure and Fr¨ohlich condensate or “non- ENFE - 02/11/2015 - Abstracts 13 equilibrium”one, which is shown to be an attractor to the of theoretical models proposed in the literature have a system. After a brief description of the system in terms rectification which decreases with increasing system size of its Hamiltonian, the presentation of the relevant vari- and thus vanishes in the thermodynamic limit. So, since ables and the associated kinetic equations, and the study these ingredients may not suffice to maintain a finite of stability of the steady-state solutions, we analyze some thermal rectification, we introduce a new ingredient, aspects of the irreversible thermodynamics of this dissipa- namely energy-conserving noise that randomly flips tive complex system. At first, the informational entropy the sign of the velocity of the system’s particles with and its production is calculated, and an order parameter a certain rate λ. With this new ingredient, we show is introduced in terms of the scaled rate of pumping and that a finite and non-zero thermal rectification in the the Fr¨ohlich parameter. Then, Glansdorff-Prigogine cri- thermodynamic limit can be obtained. Our analysis is teria for evolution and (in)stability are verified, a genera- done numerically, with the simulation of a harmonic lized H-theorem is established, a Boltzmann-like relation chain subject to a quartic local potential (pinning) and for the non-equilibrium statistical entropy is derived, as coupled at its ends to thermal reservoirs by Langevin well as expressions for the fluctuations in non-equilibrium equations. conditions and the associated Maxwell relations. [02/11/2015 - P043] [02/11/2015 - P041] Geometrical relations during coarsening for Critical properties of the susceptible- the Potts model, Marcos Paulo de O. Loureiro, exposed-infected model on a square lattice, Universidade Federal de Vi¸cosa - Campus Rio Parana´ıba, Alexander H. O. Wada, Taniaˆ Tome,´ Mario´ J. Jeferson J. Arenzon, Universidade Federal do Rio de Oliveira, Instituto de F´ısica da Universidade de Grande do Sul, Leticia F. Cugliandolo, Universit´e S˜ao Paulo The epidemiological model susceptible- Pierre et Marie Curie When taken out of equilibrium infected-exposed (SEI) is studied on a square lattice. by an instantaneous temperature quench, from above The SEI model is defined by its transition probabili- to below the critical temperature, several systems form ties, in which only susceptible sites with at least one a time evolving complex pattern (coarsening) in which infected site on its near neighborhood are allowed to several equilibrium phases compete. The energy excess change its state. A susceptible site can change its state is concentrated at the interfaces (hulls) separating with probability equal the fraction of nearest infected these several states while the curvature-driven dynamics neighbors, if so becomes infected with probability p or attempts do decrease the total length of these interfaces. exposed with 1 p. Infected and exposed sites remain Despite the very different nature of these systems, − forever in its states. Starting from the initial condition many of them satisfy the dynamic scaling hypothesis with only one infected in a lattice full of susceptibles, that states that the behavior becomes universal when the dynamics of this model grows a single cluster of a proper rescaling is performed using the characteristic infected sites until the absorbing state, in which there is length R(t) that increases in time as t1/2 when the order no pairs susceptible-infected and the cluster of infected parameter is not conserved. The study of the topological sites is completely surrounded by exposed sites. By and geometrical properties in liquid crystals, soap froths, interpreting infected and exposed sites as occupied cellular tissues, magnetic materials, superconductors and and vacant sites respectively, we show, by means of polycrystalline microstructures has attracted attention numerical simulations, that the clusters generated have for several decades. The morphology of the coarse- the same critical properties of the percolation clusters, ning patterns in these experimental systems can be in other words , both the SEI model and the percolation reproduced by the q-states Potts model. Several pheno- problem have the same critical threshold and exponents. menological laws have been proposed from the analysis Furthermore, we analyze the time series of the SEI of area and perimeter (e.g.,the Lewis and Fetham’s model up to a lattice of linear size L =215 in the critical laws) and confronted with data from two-dimensional point, calculating the dynamical critical exponents with biological tissues and metal grains. Here we follow high precision and classifying the SEI as belonging to the formation and evolution of patterns generated by the Dynamical Percolation universality class. Monte Carlo simulations of the two-dimensional Potts [02/11/2015 - P042] model for several values of q after a deep quench in Thermal Rectification in Anharmonic order to check the validity of those empirical laws in Chains under Energy-Conserving Noise, less isotropic systems and the dependence on the order Pedro Henrique Guimaraes˜ , Mario´ Jose´ de of the underlying phase transition and the amount of Oliveira, IFUSP - SP - Brasil, Gabriel T. Landi, correlation present in the initial state. UFABC - SP - Brasil Thermal rectification is the phenomenon in which the heat flux in a given system [1] Loureiro,MPO, Arenzon,JJ, Cugliandolo,LF, Si- depends on the direction the flux is applied. This cilia,A, ‘Curvature-driven coarsening in the two- phenomenon has been widely investigated in recent years dimensional Potts model’, Phys. Rev. E 81, 021129 due to its great academic and technological relevance. In (2010) order to present thermal rectification, at least two main [2] Loureiro,MPO, Arenzon,JJ, Cugliandolo,LF, ‘Geome- conditions are necessary to such systems: an inherent trical properties of the Potts model during the coarsening spacial asymmetry, which breaks the invariance under regime’, Phys. Rev. E 85, 021135 (2012) bath reversal and a temperature dependent thermal conductivity, which induces different phonon spectra Thanks FAPEMIG for financial support. when the baths are reversed. However, in disagreement with the results of experimental works, most results NONLINEAR DYNAMICS & CHAOS 14 Abstracts - ENFE - 02/11/2015
shrimp shape rotated around a focal point, called perio- [02/11/2015 - P044] dicity hubs, with intricate connections between different Exploiting the weakness of preferential attach- shrimps. ment networks, Tiago M. Vieira, Gandhi M. The purpose of this work is to further develop our recent Viswanathan, Luciano R. da Silva, UFRN - RN - study and we focus on periodicity hubs. Such hubs shows Brasil We address the general problem of how to at- focal points where your localization are associated with tack and destroy a network by node removal given limi- the generation recombination rate. In the present model ted or no prior information about the edges. Networks the coefficients of g-r processes, like as impact ioniza- have been used to describe many kinds of systems. In tion and field enhanced trapping, depend on the electric general, nodes represent systems components and edges field. Our results shows the impact when a small variati- the interactions between them. How the edges are arran- ons on some parameter of the coefficients g-r can causes ged in a network has great importance because quantities on position of the focal points and consequently over all of interest depend on edge placement, e.g. connectivity self-organized periodic structure around them. Our re- distribution, clustering coefficient, resilence to node and sults show that the position of the focal points and the edge removal, spreading processes, and small-world ef- structures of the periodic hubs have dependence with the fects. The rules controlling edge placement define the change of the coefficients g-r. The stable and unstable network structure and they can be exploited by agents hubs and many mergers between shrimps with different that wish to attack weaknesses of the networks. In our periodicity are presented. study, we consider a family of strategies in which nodes [1] da Silva, S. L., Viana, E. R., de Oliveira, A. G., Ri- are randomly chosen, but not removed. Instead a ran- beiro, G. M. and da Silva, R. L. - Int. J. Bifurcation dom acquaintance (i.e., a first neighbour) of the chosen Chaos 25, 1530004 (2015). node is removed from the network. Our approach is a ge- [2] Viana, E. R., Rubinger, R. M., Albuquerque, H. A., neralization of the strategy introduced by Cohen et. al. Dias, F. O., de Oliveira, A. G. and Ribeiro., G. M. - [Phys. Rev. Lett., 91 (2003)], in which the acquaintance Nonlinear Dyn. 67, 385-392 (2012). of a randomly chosen node is promptly removed from the network as soon as it was chosen. Instead of the imme- diate removal, a given node needs to be pointed by other [02/11/2015 - P046] Dynamical Properties of Soft Elliptical Billiard, randomly chosen nodes more than once before being re- Tiago Kroetz, Universidade Tecnol´ogica Federal do moved. As a result, we observe that our approach leads Paran´a, Pato Branco - PR - Brasil, Hercules Alves the network to be destroyed more quickly, i.e., it’s neces- de Oliveira Junior, Universidade Tecnol´ogica Fede- sary to remove a lower number of nodes, in comparison ral do Paran´a, Ponta Grossa - PR - Brasil Two- to the original strategy. dimensional billiards can be considered as special cases [02/11/2015 - P045] of two-dimensional potentials. These potentials must be Self-organization of periodicity hubs and spi- constant at the inner part of the billiard and present an rals in a high resolution parameter space from abrupt variation of their values at the coordinates on the the two-level model on semi-insulating GaAs, border of the billiard. Thus, the force exerted on a parti- Samir Lacerda da Silva, Instituto Federal do Esp´ırito cle subjected to this kind of potential is null into the billi- Santo - Campus Vit´oria, Rodrigo Lacerda da Silva, ard area and is infinite at the border of the billiard. Also, Instituto Federal Fluminense - Campus Bom Jesus do the direction of the force (and thus the potential gradient) Itabapoana, Emilson Ribeiro Viana, Universidade must be normal to the frontier of the billiard. In this work Tecnol´ogica Federal do Paran´a Semi-insulating Gal- we obtained a soft version of a two-dimensional billiard. lium Arsenide (SI-GaAs) samples experimentally show, Differently from the hard billiard, the particle confined under high electric fields and even at room tempera- in the soft billiard suffers the influence of a force during ture, negative differential conductivity in N-shaped form a time interval greater than zero. Due to this reason, the (NNDC). In recent work [1], we proposed a physical mo- particle trajectory is smooth at the reflections and differs del, the two-valley model, which describes electrical con- from straight lines between consecutive reflections. The duction in SI-GaAs. The model is based on the minimal obtaining of the soft billiard was made considering a par- set of generation-recombination (g-r) equations for two ticle subjected to a two-dimensional potential with a pa- valleys inside of the conduction band, and an equation rameter capable to change the values of gradient function for the drift velocity as a function of the applied elec- without however alter the shape of equipotential curves. tric field, that covers the physical properties of the non- With this, we can investigate the continuous transition linear electrical conduction of the SI-GaAs system. The of the dynamics from soft two-dimensional potential to model generated theoretically the NNDC region for the the corresponding hard two-dimensional billiard. We opt first time and the nonlinear dynamics were investigated to perform this investigation considering the elliptical ge- in this region by building of high-resolution parameter- ometry of the equipotential curves, where the values of space of the periodicity using a Periodicity-Detection rou- eccentricities are the same for each equipotential and can tine [2]. In the parameter-space we observed too many be controlled by a parameter in the potential expression. self-organized periodic structures embedded in the cha- Using this procedure we can reveal the changes of the otic regions, like as a “shrimp”shaped in a spiral form, numerical results by varying the hardness of the border that forms a “snail shell”. This structure established until recover the well known phase space of hard ellip- a direct communication between the windows in order tical billiard. We investigate the two-dimensional space within chaotic regions, producing new routes of bifurca- of parameters identifying the transitions order-chaos in tion. The snail structure show three regions where the there. ENFE - 02/11/2015 - Abstracts 15
[02/11/2015 - P047] the financial support of FAPESP (grants 2014/07043-0 Three unequal masses on a ring and soft triangu- and 2011/19296-1). lar billiards, Hercules A. Oliveira, Universidade Tecnol´ogica Federal do Paran´a, Ponta Grossa, Marcus [02/11/2015 - P049] W. Beims, Universidade Federal do Paran´a Collisions Chaos in 3D Brans-Dicke Model, with hard (infinite) walls in billiards systems are usually Thiago Gilberto do Prado,Hercules´ Alves described by instantaneous reversal of the particles linear de Oliveira Junior, Marcos Cesar Verges,` Uni- momentum. From this, simple analytical relations of ve- versidade Tecnol´ogica Federal do Paran´a The evidence locities and angles before and after the collisions with of an accelerated expansion of the universe, observed by the walls are obtained. However, in order to analyze the WMAP, opened the discussion of the general relativity transition to soft walls, which are more realistic, it is es- eventual limits. There are many options for alternative sential to have well defined equations of motion since, theories of gravity, and among them wecan cite scalar- in general, no simple analytical solutions are obtained. tensor theories like supergravity, Kaluza-Klein theories, The present work suggests that appropriated soft walls dual string theories, M-Theory, etc . One particular potentials are those for which the corresponding forces kind of scalar-tensor theory to describe an accelerate become “delta functions”in the limit of hard walls. This expansion of the universe, called Brans-Dicke theory, was allows for better numerical investigation of the soft-hard proposed in the early sixties. This theory uses the princi- transition. A general scaled Hamiltonian is derived for ple of Mach and the hypothesis of Dirac, considering an three unequal masses interacting particles on a friction- eventual variation in time of the Newton’s gravitational less ring, which nicely describes the transition and shows constant, thus ensuring the universality of free fall that the dynamics occurs inside a soft triangular billiard. (equivalence principle). Most of the works which have The dynamics of three soft interacting particles on a ring been published in this theory up to now consider four is shown to correspond to the motion of one particle in- flat dimensions, and some of them have tried to associate side a soft triangular billiard. The dynamics inside the the scalar field of the Brans-Dicke theory as quintessence soft billiard depends only on the masses ratio between field as a type of K-essence field. Others have tried to particles and softness ratio of the particles interaction. find a solution for the observed accelerated expansion The transition from soft to hard interactions can be ap- using a dimensional reduction of the 5D Brans-Dicke propriately explored using potentials for which the cor- theory without matter. Concerning 3 dimensions, a responding equations of motion are well defined in the broad study has been done in gravitational theories since hard wall limit. Numerical examples are shown for the the publication of BTZ Black Hole, motivated by the soft Toda-like interaction and the error function. fact that 3D theories avoid some complications found in higher dimensions. However, there are not as many [02/11/2015 - P048] results about 3D scalar-tensor theories, and it would Plasma Structures in Texas Helimak, be interesting to find some results in this subject, more F. A. C. Pereira, Z. O. Guimaraes-Filho,˜ I. specifically in 3D Brans-Dicke theory. For instance, we L. Caldas, Instituto de F´ısica da Universidade de S˜ao can see some problems like the association of the scalar Paulo, D. L. Toufen, Instituto Federal de Educa¸c˜ao, field of the Brans-Dicke theory to K-essence fields which Ciˆencia e Tecnologia de S˜ao Paulo, Campus Guarulhos, models the dark energy. The goal in this paper is study K. W. Gentle, Institute for Fusion Studies of the the dynamical behavior of a 3D Brans-Dicke model. For University of Texas at Austin Intermittent structures this we use tools such as calculation of the exponent with extreme events (bursts) have been detected in the Lyaponov and analysis of phase space for the study of turbulence of the scrape of layer (SOL) of Tokamaks this dynamic system. From the results we identified and these structures have a major role in plasma confi- regions whose theory has chaotic behavior. nement. The Texas Helimak is a toroidal plasma device with one-dimensional equilibrium, magnetic curvature [02/11/2015 - P050] and shear, thus resembling closely a SOL of a Tokamak. Algebraic period 3 orbits: rotation on The Texas Helimak vacuum vessel has a rectangular complex plane and statistical quantities., cross section with 0.6 m of internal radius, 1.6 m of Antonioˆ Joao˜ Fidelis´ , IFC - Campus Luzerna, Lu- external radius, and 2 m of height and the plasma ciano Camargo Martins, UDESC - CCT For the is generated by electron cyclotron resonance heating. period-3 window of the logistic map xn+1 = rxn(1 xn) The Helimak has 16 bias plates, where a large set of it is presented algebraic orbits, for both stable− and Langmuir probes is mounted and from where is possible unstable ones. From the solution of xn+3 = xn it is to impose an external electrical bias. This electrical bias obtained the polynomial that rules the periodic orbits can change the turbulence properties and even suppress inside this window. The roots of this polynomial are the appearance of the intermittent structures. In this the orbits, and they are functions of the fixed parameter work, we study the statistical characteristics of bursts value r. As r is increased, the value of the roots are in the ion saturation current, and we show that a shot modified: some increase and the others decrease. For the noise model can reproduce the signal statistics. We also same branch of the three possible orbits, the value of the characterize both the spacial and temporal profiles of roots present opposite behavior with respect to stable the bursts and their propagation, as a function of the and unstable orbits. The roots of the polynomial, i.e., external imposed bias and the radial position inside the orbits, are presented in two different ways: a sum the Helimak. Both profiles and the bursts propagation of complex numbers xi = a + bc + bc, and via Euler’s are analyzed by conditional analysis and the cross formula xi = a + 2 b cos(θ) – the overbar indicates correlation between the probes signals. We acknowledge complex conjugation.| The| algebraic orbits are obtained 16 Abstracts - ENFE - 02/11/2015 for three different fixed control parameter values of r: 2t and height 2b, as usual, we set b = 1 [1]. Here, we at tangent bifurcation (birth), at super-stability and study numerically the vicinity of a particular line in the at ending pitchfork bifurcation (death). The algebraic paramater space a t, namely t = t (a) = √a2 1. If × c 0 − expressions of the constants a, b, c, b and θ are given t > tc, there is chaos almost everywhere[1]. If t