Variational Inference MPC Using Tsallis Divergence

Total Page:16

File Type:pdf, Size:1020Kb

Variational Inference MPC Using Tsallis Divergence Variational Inference MPC using Tsallis Divergence Ziyi Wang∗y, Oswin So∗, Jason Gibson, Bogdan Vlahov, Manan S. Gandhi Guan-Horng Liu and Evangelos A. Theodorou Autonomous Control and Decision Systems Lab Georgia Institute of Technology, Atlanta, Georgia ∗These authors contributed equally yCorrespondence to: [email protected] Abstract—In this paper, we provide a generalized framework the variational distribution and target distribution. for Variational Inference-Stochastic Optimal Control by using the In existing VI-SOC works, the KL divergence is used as non-extensive Tsallis divergence. By incorporating the deformed the distributional distance metric due to its simplicity. On the exponential function into the optimality likelihood function, a novel Tsallis Variational Inference-Model Predictive Control other hand, recent advances in VI research involve extending algorithm is derived, which includes prior works such as Varia- the framework to other statistical divergences, such as the tional Inference-Model Predictive Control, Model Predictive Path α-divergence [7] and χ-divergence [8]. In Wang et al. [9], Integral Control, Cross Entropy Method, and Stein Variational Regli and Silva [10], the authors proposed variants of the Inference Model Predictive Control as special cases. The pro- α-divergence to improve the performance and robustness of posed algorithm allows for effective control of the cost/reward transform and is characterized by superior performance in terms the inference algorithm. Wan et al. [11] further extended the of mean and variance reduction of the associated cost. The VI framework to f-divergence, which is a broad statistical aforementioned features are supported by a theoretical and divergence family that recovers the KL, α and χ-divergence numerical analysis on the level of risk sensitivity of the proposed as special cases. algorithm as well as simulation experiments on 5 different robotic Tsallis divergence is another generalized divergence rooted systems with 3 different policy parameterizations. in non-extensive statistical mechanics [12]. Tsallis divergence and Tsallis entropy are generalizations of the KL divergence I. INTRODUCTION and Shannon entropy respectively. The Tsallis entropy is non- Variational Inference (VI) is a powerful tool for approxi- additive (hence non-extensive) in the sense that for a system mating the posterior distribution of the unobserved random composed of (probabilistically) independent subsystems, the variables [1]. VI recasts the approximation problem as an total entropy differs from the sum of the entropies of the optimization problem. Instead of directly approximating the subsystems [13]. Over the last two decades, an increasing target distribution p(zjx) of the latent variable z, VI minimizes number of complex natural, artificial and social systems the Kullback-Leibler (KL) divergence between a tractable have verified the predictions and consequences derived from variational distribution q(z) and the target distribution. Due the Tsallis entropy and divergence [14, 15, 16, 17]. Lee to its faster convergence and comparable performance to et al. [17] demonstrated that Tsallis entropy regularization Markov Chain Monte Carlo sampling methods, VI has received leads to improved performance and faster convergence in RL increasing attention in machine learning and robotics [2, 3, 4]. applications. VI has been applied to Stochastic Optimal Control (SOC) In this paper, we provide a generalized formulation of the problems recently. In Okada and Taniguchi [5], the authors VI-SOC framework using the Tsallis divergence and introduce formulated the SOC problem as a VI problem by setting a novel Model Predictive Control (MPC) algorithm. The main arXiv:2104.00241v1 [cs.LG] 1 Apr 2021 the desired policy distribution as the target distribution. The contribution of our work is threefold: VI-SOC framework works directly in the space of policy • We propose the Tsallis VI-MPC algorithm, which allows distributions instead of specific policy parameterizations in most for additional control of the shape of the cost transform SOC and Reinforcement Learning (RL) frameworks. This gives compared to previous VI-MPC algorithms using KL rise to a unified derivation for a variety of parametric policy dis- divergence. tributions, such as unimodal Gaussian and Gaussian mixture in • We provide a holistic view of connections between Tsallis Okada and Taniguchi [5]. Lambert et al. [6] derived the VI-SOC VI-SOC and the state-of-the-art Model Predictive Path algorithm for non-parametric policy distribution using Stein Integral (MPPI) control, Cross Entropy Method (CEM) variational gradient descent. Apart from the policy distribution, and Stochastic Search (SS) methods. the VI-SOC framework is characterized by two components, the • We show, both analytically and numerically, that the optimality likelihood function and distributional distance metric. proposed Tsallis VI-SOC framework achieves lower cost The optimality likelihood function measures the likelihood of variance than MPPI and CEM. We further demonstrate trajectory samples to be optimal and defines the cost/reward the superior performance of Tsallis VI-SOC in mean cost transform to allow for different algorithmic developments. The and variance minimization on simulated systems from distributional distance metric is the measure of distance between control theory and robotics for 3 different choices of policy distributions. TABLE I: Notation comparison between Variational Inference and Variational Inference-Stochastic Optimal Control. The rest of this paper is organized as follows: in SectionII, we review KL VI-SOC and derive the Tsallis VI-SOC frame- VI VI-SOC work. In Section III, we discuss the connections between Tsallis Notation Meaning Notation Meaning VI-SOC and related works. A reparameterization and analysis x data o optimality of the framework is included in SectionIV, and we propose the z latent variable τ trajectory novel Tsallis VI-MPC algorithm in SectionV. We showcase the p(z) prior distribution p(τ) prior distribution performance of the proposed algorithm against related works q(z) variational distribution q(τ) controlled distribution in SectionVI and conclude the paper in Section VII. p(xjz) generative model p(ojτ) optimality likelihood p(zjx) posterior distribution p(τjo) optimal distribution II. TSALLIS VARIATIONAL INFERENCE-STOCHASTIC OPTIMAL CONTROL In this section, we review the VI-SOC formulation using KL where the constant p(o) is dropped. The first term in the divergence [5] and derive the novel Tsallis divergence VI-SOC objective measures the likelihood of a trajectory being optimal framework. while the second term serves as regularization. Splitting the A. SOC Problem Formulation expectation in the first term, we get ∗ n In discrete time SOC, we work with trajectories τ := fX; Ug q (U) = arg min − Eq(U)[Ep(XjU)[log p(ojτ)]] nx×T +1 q(U) of state, X := fx0; x1; : : : ; xT g 2 R , and control, (3) o : nu×T U = fu0; u1; : : : ; uT −1g 2 R over a finite time horizon + DKL (q(U) k p(U)) ; T > 0. The goal in SOC problems is to minimize the QT −1 expected cost defined by an arbitrary cost function J : where p(U) and q(U) represent t=0 p(ut) and nx×T × nu×T −1 ! + with initial state distribution p(x ) QT −1 R R R 0 t=0 q(ut) respectively due to independence and T −1 and state transition probability p(xt+1jxt; ut) corresponding Q p(XjU) = p(x0) t=0 p(xt+1jxt; ut). n to stochastic dynamics xt+1 = F (xt; ut; t) where t 2 R Optimality Likelihood: The optimality likelihood in Eq. (3) is the system stochasticity. can be parameterized by a non-increasing function of the trajec- B. KL VI-SOC tory cost J(X; U) as p(ojτ) := f(J(X; U)). The monotonicity requirement ensures that trajectories incurring higher costs We can formulate the SOC problem as an inference problem are always less likely to be optimal. Common choices of f and apply methods from VI. To apply VI to SOC, we introduce include f(x) = exp(−x) and f(x) = 1fx≤γg. In this paper, a dummy optimality variable o 2 f0; 1g with o = 1 indicating we choose f(x) = exp(−x) such that log p(ojτ) = −J(X; U). that the trajectory τ = fX; Ug is optimal. TableI compares To avoid excessive notation, we abuse the notation to define the differences in notation between conventional VI methods J(U) := Ep(XjU)[J(X; U)] and the N-sample empirical and SOC formulated as a VI problem. approximation of the expectation J = PN [J (n)]. Hence, The objective of VI-SOC is to sample from the posterior n=1 Eq. (3) takes the form distribution ∗ p(o = 1jτ)p(τ) q (U) = arg min Eq(U)[J(X; U)] + DKL (q(U) k p(U)) ; p(τjo) = q(U) p(o) (4) T −1 (1) p(o = 1jτ) Y which can be interpreted as an application of the famous = p(x ) p(x jx ; u )p(u ): p(o = 1) 0 t+1 t t t maximum entropy principle. t=0 C. Tsallis VI-SOC Here p(xt+1jxt; ut) is the state transition probability, p(x0) is the initial state distribution and p(ut) is some prior control In this subsection, we use the Tsallis divergence as the reg- distribution (e.g. zero mean Gaussian or uniform distribution). ularization function and derive the Tsallis-VI-SOC algorithm. For simplicity, we use o to indicate o = 1 and o0 for o = 0 First, we define the deformed logarithm and exponential as from here on. We can now formulate the VI-SOC objective as xr−1 − 1 log (x) = ; (5) minimizing the distance between a controlled distribution q(τ) r r − 1 and the target distribution p(τjo) 1 exp (x) = (1 + (r − 1)x) r−1 ; (6) ∗ r + q (τ) = arg min DKL (q(τ) k p(τjo)) q(τ) where (·)+ := max(0; ·) and r > 0. Using logr and expr we " # p(XjU)q(U) can define the Tsallis entropy and the corresponding Tsallis = arg min log Eq(τ) p(ojτ) divergence [12] as q(τ) (2) p(o) p(XjU)p(U) Sr(q(z)) := − q[log q(z)] " T −1 # E r X q(ut) Z = arg min − log p(ojτ) + log ; 1 r (7) Eq(τ) p(u ) = − q(z) dz − 1 ; q(τ) t=0 t r − 1 and by the empirical distribution q~∗(U) with weights w(n): q(z) Dr (q(z) k p(z)) :=Eq logr N p(z) ∗ X (n) q~ (U) = w 1 (n) ; (12) r−1 ! fU=U g 1 Z q(z) n=1 = q(z) dz − 1 : r − 1 p(z) q∗(U (n)) w(n) = : (13) (8) PN ∗ (n0) n0=1 q (U ) Note that as r ! 1, logr ! log, expr ! exp, Dr !DKL We now solve Eq.
Recommended publications
  • Luis Gregorio Moyano
    Luis Gregorio Moyano Rep´ublicadel L´ıbano 150, Mendoza, 5500 Mendoza, Argentina [email protected] Current position Facultad de Ciencias Exactas y Naturales Mendoza, Argentina Universidad Nacional de Cuyo Adjunt Professor May 2016 - Today CONICET Adjunt Researcher May 2016 - Today Education and training Universidad Carlos III de Madrid Madrid, Spain Postdoctoral Position at 2007 - 2008 Mathematics Department & GISC Centro Brasileiro de Pesquisas F´ısicas Rio de Janeiro, Brazil PhD in Physics 2001 - 2006 Supervisor: Prof. Constantino Tsallis. PhD Thesis: \Nonextensive statistical mechanics in complex systems: dynamical foundations and applications" Instituto Balseiro Bariloche, Argentina MSc in Physics 1997- 2000 Supervisors: Prof. Dami´anZanette and Prof. Guillermo Abramson. Msc. Thesis: \Learning in coupled dynamical systems" Research Interests • Machine Learning applications on networks, Natural Language Processing, Big Data applications. • Complex Networks, Social Networks, Complex Systems. • Statistical mechanics, applications to economical, social and biological systems. • Bayesian inference, Markovian systems, congestion and information transfer dynamics in networks. • Evolutionary game theory, emergence of cooperation in complex networks. • Nonlinear dynamics, coupled chaotic systems, synchronization. Luis G. Moyano - Curriculum Vitæ 1 Positions, fellowships and awards • IBM Research Brazil, Rio de Janeiro, Brazil - Staff Research Member December, 2013 - April, 2016 Staff Research Member at the Social Data Analytics
    [Show full text]
  • Associating an Entropy with Power-Law Frequency of Events
    entropy Article Associating an Entropy with Power-Law Frequency of Events Evaldo M. F. Curado 1, Fernando D. Nobre 1,* and Angel Plastino 2 1 Centro Brasileiro de Pesquisas Físicas and National Institute of Science and Technology for Complex Systems, Rua Xavier Sigaud 150, Rio de Janeiro 22290-180, Brazil; [email protected] 2 La Plata National University and Argentina’s National Research Council (IFLP-CCT-CONICET)-C. C. 727, La Plata 1900, Argentina; plastino@fisica.unlp.edu.ar * Correspondence: [email protected]; Tel.: +55-21-2141-7513 Received: 3 October 2018; Accepted: 23 November 2018; Published: 6 December 2018 Abstract: Events occurring with a frequency described by power laws, within a certain range of validity, are very common in natural systems. In many of them, it is possible to associate an energy spectrum and one can show that these types of phenomena are intimately related to Tsallis entropy Sq. The relevant parameters become: (i) The entropic index q, which is directly related to the power of the corresponding distribution; (ii) The ground-state energy #0, in terms of which all energies are rescaled. One verifies that the corresponding processes take place at a temperature Tq with kTq ∝ #0 (i.e., isothermal processes, for a given q), in analogy with those in the class of self-organized criticality, which are known to occur at fixed temperatures. Typical examples are analyzed, like earthquakes, avalanches, and forest fires, and in some of them, the entropic index q and value of Tq are estimated. The knowledge of the associated entropic form opens the possibility for a deeper understanding of such phenomena, particularly by using information theory and optimization procedures.
    [Show full text]
  • The Entropy Universe
    entropy Review The Entropy Universe Maria Ribeiro 1,2,* , Teresa Henriques 3,4 , Luísa Castro 3 , André Souto 5,6,7 , Luís Antunes 1,2 , Cristina Costa-Santos 3,4 and Andreia Teixeira 3,4,8 1 Institute for Systems and Computer Engineering, Technology and Science (INESC-TEC), 4200-465 Porto, Portugal; [email protected] 2 Computer Science Department, Faculty of Sciences, University of Porto, 4169-007 Porto, Portugal 3 Centre for Health Technology and Services Research (CINTESIS), Faculty of Medicine University of Porto, 4200-450 Porto, Portugal; [email protected] (T.H.); [email protected] (L.C.); [email protected] (C.C.-S.); andreiasofi[email protected] (A.T.) 4 Department of Community Medicine, Information and Health Decision Sciences-MEDCIDS, Faculty of Medicine, University of Porto, 4200-450 Porto, Portugal 5 LASIGE, Faculdade de Ciências da Universidade de Lisboa, 1749-016 Lisboa, Portugal; [email protected] 6 Departamento de Informática, Faculdade de Ciências da Universidade de Lisboa, 1749-016 Lisboa, Portugal 7 Instituto de Telecomunicações, 1049-001 Lisboa, Portugal 8 Instituto Politécnico de Viana do Castelo, 4900-347 Viana do Castelo, Portugal * Correspondence: [email protected] Abstract: About 160 years ago, the concept of entropy was introduced in thermodynamics by Rudolf Clausius. Since then, it has been continually extended, interpreted, and applied by researchers in many scientific fields, such as general physics, information theory, chaos theory, data mining, and mathematical linguistics. This paper presents The Entropy
    [Show full text]
  • Nonextensive Statistical Mechanics: a Brief Review of Its Present Status
    Anais da Academia Brasileira de Ciências (2002) 74(3): 393–414 (Annals of the Brazilian Academy of Sciences) ISSN 0001-3765 www.scielo.br/aabc Nonextensive statistical mechanics: a brief review of its present status CONSTANTINO TSALLIS* Centro Brasileiro de Pesquisas Físicas, 22290-180 Rio de Janeiro, Brazil Centro de Fisica da Materia Condensada, Universidade de Lisboa P-1649-003 Lisboa, Portugal Manuscript received on May 25, 2002; accepted for publication on May 27, 2002. ABSTRACT We briefly review the present status of nonextensive statistical mechanics. We focus on (i) the cen- tral equations of the formalism, (ii) the most recent applications in physics and other sciences, (iii) the a priori determination (from microscopic dynamics) of the entropic index q for two important classes of physical systems, namely low-dimensional maps (both dissipative and conservative) and long-range interacting many-body hamiltonian classical systems. Key words: nonextensive statistical mechanics, entropy, complex systems. 1 CENTRAL EQUATIONS OF NONEXTENSIVE STATISTICAL MECHANICS Nonextensive statistical mechanics and thermodynamics were introduced in 1988 [1], and further developed in 1991 [2] and 1998 [3], with the aim of extending the domain of applicability of statistical mechanical procedures to systems where Boltzmann-Gibbs (BG) thermal statistics and standard thermodynamics present serious mathematical difficulties or just plainly fail. Indeed, a rapidly increasing number of systems are pointed out in the literature for which the usual functions appearing in BG statistics appear to be violated. Some of these cases are satisfactorily handled within the formalism we are here addressing (see [4] for reviews and [5] for a regularly updated bibliography which includes crucial contributions and clarifications that many scientists have given along the years).
    [Show full text]
  • Q-Gaussian Approximants Mimic Non-Extensive Statistical-Mechanical Expectation for Many-Body Probabilistic Model with Long-Range Correlations
    Cent. Eur. J. Phys. • 7(3) • 2009 • 387-394 DOI: 10.2478/s11534-009-0054-4 Central European Journal of Physics q-Gaussian approximants mimic non-extensive statistical-mechanical expectation for many-body probabilistic model with long-range correlations Research Article William J. Thistleton1∗, John A. Marsh2† , Kenric P. Nelson3‡ , Constantino Tsallis45§ 1 Department of Mathematics, SUNY Institute of Technology, Utica NY 13504, USA 2 Department of Computer and Information Sciences, SUNY Institute of Technology, Utica NY 13504, USA 3 Raytheon Integrated Defense Systems, Principal Systems Engineer 4 Centro Brasileiro de Pesquisas Fisicas, Rua Xavier Sigaud 150, 22290-180 Rio de Janeiro - RJ, Brazil 5 Santa Fe Institute, 1399 Hyde Park Road, Santa Fe, NM 87501, USA Received 3 November 2008; accepted 25 March 2009 Abstract: We study a strictly scale-invariant probabilistic N-body model with symmetric, uniform, identically distributed random variables. Correlations are induced through a transformation of a multivariate Gaussian distribution with covariance matrix decaying out from the unit diagonal, as ρ/rα for r =1, 2, …, N-1, where r indicates displacement from the diagonal and where 0 6 ρ 6 1 and α > 0. We show numerically that the sum of the N dependent random variables is well modeled by a compact support q-Gaussian distribution. In the particular case of α = 0 we obtain q = (1-5/3 ρ) / (1- ρ), a result validated analytically in a recent paper by Hilhorst and Schehr. Our present results with these q-Gaussian approximants precisely mimic the behavior expected in the frame of non-extensive statistical mechanics.
    [Show full text]
  • Some Comments on Boltzmann-Gibbs Statistical Mechanics
    CBPF-NF-001/94 Some Comments on Boltzmann-Gibbs Statistical Mechanics Constantino TSALLIS Centro Brasileiro de Pesquisas Fsicas Rua Dr. Xavier Sigaud, 150 22290-180 { Rio de Janeiro { RJ, Brazil Abstract A nonexhaustive review is presented of the limits of the impressive and vastly known success of Boltzmann-Gibbs statistics and normal thermo dynamics. These limits naturally op en the do or for the research of generalized formalisms that could enlarge the domain of validity of standard statistical mechanics and thermo dynam- ics. A p ossible such generalization (recently prop osed by the author) is commented along this p ersp ective. CBPF-NF-001/94 1 1 Limitations of Boltzmann-Gibbs Statistics 1.1 Intro duction The qualitative and quantitative success of Boltzmann-Gibbs (BG) Statistical Mechan- ics (and, naturally, of its particular cases, the Fermi-Dirac and Bose-Einstein quantum statistics, with their common high temp erature asymptotic limit, the classical Maxwell- Boltzmann statistics) is so ubiquitous, persistent and delicate that not few physicists and chemists have a kind of strong (not necessarily rationalized) feeling that this brilliant for- malism is, in practical terms, universal, eternal and in nitely precise. A more balanced analysis reveals some of its sp eci c limitations and inadequacies, and consequently the fragility or nonuniversality of some of the basic hyp othesis of its foundations. Among these foundation stones, a privileged p osition is detained by the entropy S ,asintro duced R and used by Boltzmann and Gibbs (S = k dxf (x)lnf (x) for classical distribution B laws f (x) de ned in phase space), further generalized byvon Neumann (S = k Tr ln , B b eing the density op erator de ned in Hilb ert or Fock spaces), with its diagonal form W (S = k p ln p , p b eing the probability of the i-th among W microstates, and its B i i i i=1 famous equal-probability particular form S = k ln W ) as nely discussed, in the context B of Information Theory,by Shannon.
    [Show full text]
  • Curriculum Vitae
    CURRICULUM VITAE Nome e sobrenome / Name: CONSTANTINO TSALLIS Nascimento / Birth: 05/11/1943, Atenas - Grecia Nacionalidade / Citizenship: Brasileira / Grega Familia / Family: 3 filhos Endere¸coprofissional / Office address: Centro Brasileiro de Pesquisas Fisicas / MCTI Rua Xavier Sigaud 150 - Urca 22290-180 Rio de Janeiro - Brasil Tel.: (55)(21) 2141 7190 / 2141 7213 Fax: (55)(21) 2141 7515 / 2141 7400 E-mail: [email protected] Internet: http://tsallis.cat.cbpf.br Endere¸coresidencial / Home address: Av. Epitacio Pessoa 4344 - Apto. 502 - Bl.B - Lagoa 22471-001 Rio de Janeiro - Brasil Tel.: (55)(21) 2286 9503 / 98152 5051 Carteira de Identidade / Identity card: RG 5.173.806-0 (expedida em 04/01/85 - Instituto Felix Pacheco-RJ) C.P.F./ Social Security: 119.361.181-49 Cargo atual / Present position (CBPF): Pesquisador Emerito / Emeritus Professor 1 DISTINCOES~ / HONORS AND DISTINCTIONS • Chefias Departamento / CBPF [Head of Department] Fisica Teorica (1979 - 1980) Fisica da Materia Condensada (1981 - 1983) Fisica da Materia Condensada e Fisica Estatistica (1991 - 1993) Fisica da Materia Condensada e Fisica Estatistica (1999 - 2000) Fisica Teorica (2009 - 2013) • Guggenheim Fellow (USA, l982) • Pesquisador do CNPq N´ıvel 1 A (1984 - ...) [CNPq highest granting level] • "Distinguished Scientist of the Greek Diaspora", Governo da Gr´ecia(1985) • Membro do "International Advisory Committee " STATPHYS 16 / IUPAP (USA, 1986) STATPHYS 17 / IUPAP (Brasil, 1989) [Conference Chairman] STATPHYS 18 / IUPAP (Alemanha, 1992) STATPHYS 19 / IUPAP (China, 1995)
    [Show full text]
  • Tsallis Entropy, Likelihood, and the Robust Seismic Inversion
    entropy Article Tsallis Entropy, Likelihood, and the Robust Seismic Inversion Igo Pedro de Lima 1,†, Sérgio Luiz E. F. da Silva 2,*,† , Gilberto Corso 1,3,† and João M. de Araújo 1,2,† 1 Programa de Pós-Graduação em Ciência e Engenharia de Petróleo - Universidade Federal do Rio Grande do Norte, Natal RN 59078-970, Brazil; [email protected] (I.P.d.L.); [email protected] (G.C.); joaomedeiros@fisica.ufrn.br (J.M.d.A.) 2 Departamento de Física Teórica e Experimental, Universidade Federal do Rio Grande do Norte, Natal RN 59078-970, Brazil 3 Departamento de Biofísica e Farmacologia, Universidade Federal do Rio Grande do Norte, Natal RN 59078-970, Brazil * Correspondence: [email protected] † These authors contributed equally to this work. Received: 2 April 2020; Accepted: 17 April 2020; Published: 19 April 2020 Abstract: The nonextensive statistical mechanics proposed by Tsallis have been successfully used to model and analyze many complex phenomena. Here, we study the role of the generalized Tsallis statistics on the inverse problem theory. Most inverse problems are formulated as an optimisation problem that aims to estimate the physical parameters of a system from indirect and partial observations. In the conventional approach, the misfit function that is to be minimized is based on the least-squares distance between the observed data and the modelled data (residuals or errors), in which the residuals are assumed to follow a Gaussian distribution. However, in many real situations, the error is typically non-Gaussian, and therefore this technique tends to fail. This problem has motivated us to study misfit functions based on non-Gaussian statistics.
    [Show full text]
  • Thermodynamics and Statistical Mechanics for Complex Systems – Foundations and Applications∗
    THERMODYNAMICS AND STATISTICAL MECHANICS FOR COMPLEX SYSTEMS { FOUNDATIONS AND APPLICATIONS∗ Constantino Tsallis Centro Brasileiro de Pesquisas F´ısicas and National Institute of Science and Technology for Complex Systems Rua Xavier Sigaud 150, 22290-180 Rio de Janeiro, Brazil and Santa Fe Institute, 1399 Hyde Park Road, Santa Fe, New Mexico 87501, USA The complexity (or simplicity) of a system can be characterized through very many paths. A quite interesting and operational one, which we follow here, consists in making use of entropic concepts. The thermodynamical and statistical-mechanical foundations of this approach, as well as its (an- alytical, computational, observational, and experimental) applications for natural, artificial and social systems will be briefly reviewed. 1. Introduction In the 1860s Clausius introduced [1] the concept of entropy (noted S) within thermodynamics. This most relevant step basically concluded the foundations of this science, and enabled its expression in terms of the so called principles of thermodynamics (nowadays referred to as zeroth, first, second, third principles). One decade later Boltzmann postulated [2,3] an entropic functional (noted SBG where BG stands for Boltzmann-Gibbs, in order to also incorporate the refinements introduced by Gibbs [4]), which connects the macroscopic entropy introduced by Clausius with the micro- scopic world. The entropy SBG is essentially consistent with the entire structure of classical thermodynamics, and constitutes the basis of the the- ory (very important in contemporary physics) referred to as statistical me- chanics. During nearly 140 years, BG statistical mechanics keeps providing uncountable excellent predictions for many kinds of thermal and similar phenomena. This has somehow sedimented { wrongly! { in the mind of ∗ Invited review presented at the XXVIIth Marian Smoluchowski Symposium on Sta- tistical Physics, held in Zakopane, Poland, in 22-26 September 2014.
    [Show full text]
  • Tsallis 70 Complex Systems Foundations and Applications
    Complex Systems Foundations and Applications Scientific Program Also: A celebration of the 70th birthday of Constantino Tsallis CBPF - Rio de Janeiro - Brazil From October 29 to November 1 – 2013 The digital version of this book of abstracts can be downloaded at: http://www.cbpf.br/~compsyst/index.php Cover photo: C. Tsallis in his CBPF office, 2013, Rio de Janeiro, Brazil. Credits: J. Ricardo. Foreword Dear Colleagues, We are very pleased to meet you here in CBPF, Rio de Janeiro, for the Complex Systems – Foundations and Applications conference. This will be an international event covering many topics in the area of complex systems, like nonlinear phenomena, econophysics, foundations and ap- plications of non-extensive statistical mechanics, biological complexity, far-from-equilibrium phenomena, nonequilibrium in social and natural sciences, and interdisciplinary applications, among others. In this event, we will also celebrate the 70th birthday of Constantino Tsallis. We wish you all a great conference. Andr´eM. C. Souza (UFS) Evaldo M. F. Curado (CBPF) Fernando D. Nobre (CBPF) Roberto F. S. Andrade (UFBA) ORGANIZING COMMITTEE Scientific Program DAY I DAY II DAY III DAY IV TUESDAY WEDNESDAY THURSDAY FRIDAY 29/10/2013 30/10/2013 31/10/2013 01/11/2013 08h00 – 09h00 Registration 09h00 – 09h30 Opening ORAL COMMUNICATIONS 09h30 – 10h30 ORAL COMMUNICATIONS 10h30 – 11h00 Coffee Break 11h00 – 12h30 ORAL COMMUNICATIONS 12h30 – 14h30 LUNCH 14h30 – 16h00 ORAL COMMUNICATIONS 16h00 – 16h30 Coffee Break Coffee Break and Coffee Break 16h30 – 17h00 ORAL COM. POSTER SESSIONS‡ ORAL COM. 17h00 – 17h30 ORAL COMMUNICATIONS C. TSALLIS ORAL COM. 17h30 – 18h00 ORAL COMMUNICATIONS SPECIAL Closure 18h00 – 18h30 Cocktail TALK 20h00 Dinner All talks: 25 minutes (presentation) + 5 minutes (questions).
    [Show full text]
  • Nonadditive Entropy and Nonextensive Statistical Mechanics - an Overview After 20 Years
    Brazilian Journal of Physics, vol. 39, no. 2A, August, 2009 337 Nonadditive entropy and nonextensive statistical mechanics - An overview after 20 years Constantino Tsallis Centro Brasileiro de Pesquisas F´ısicas and National Institute of Science and Technology for Complex Systems Rua Xavier Sigaud 150, 22290-180 Rio de Janeiro, Brazil and Santa Fe Institute 1399 Hyde Park Road, Santa Fe, NM 87501, USA (Received on 17 April, 2009) Statistical mechanics constitutes one of the pillars of contemporary physics. Recognized as such — together with mechanics (classical, quantum, relativistic), electromagnetism and thermodynamics —, it is one of the mandatory theories studied at virtually all the intermediate- and advanced-level courses of physics around the world. As it normally happens with such basic scientific paradigms, it is placed at a crossroads of various other branches of knowledge. In the case of statistical mechanics, the standard theory — hereafter referred to as the Boltzmann-Gibbs (BG) statistical mechanics — exhibits highly relevant connections at a variety of microscopic, mesoscopic and macroscopic physical levels, as well as with the theory of probabilities (in particular, with the Central Limit Theorem, CLT). In many circumstances, the ubiquitous efects of the CLT, with its Gaussian attractors (in the space of the distributions of probabilities), are present. Within this complex ongoing frame, a possible generalization of the BG theory was advanced in 1988 (C.T., J. Stat. Phys. 52, 479). The extension of the standard concepts is intended to be useful in those “pathological”, and nevertheless very frequent, cases where the basic assumptions (molecular chaos hypothesis, ergodicity) for applicability of the BG theory would be violated.
    [Show full text]
  • Black Hole Entropy: a Closer Look
    entropy Review Black Hole Entropy: A Closer Look Constantino Tsallis 1,2,3 1 Centro Brasileiro de Pesquisas Físicas and National Institute of Science and Technology of Complex Systems, Rua Xavier Sigaud 150, Rio de Janeiro 22290-180, RJ, Brazil; [email protected] 2 Santa Fe Institute, 1399 Hyde Park Road, Santa Fe, NM 87501, USA 3 Complexity Science Hub Vienna, Josefstädter Strasse 39, 1080 Vienna, Austria Received: 25 November 2019; Accepted: 19 December 2019; Published: 22 December 2019 Abstract: In many papers in the literature, author(s) express their perplexity concerning the fact that the (3 + 1) black-hole ‘thermodynamical’ entropy appears to be proportional to its area and not to its volume, and would therefore seemingly be nonextensive, or, to be more precise, subextensive. To discuss this question on more clear terms, a non-Boltzmannian entropic functional noted Sd was applied [Tsallis and Cirto, Eur. Phys. J. C 73, 2487 (2013)] to this complex system which exhibits the so-called area-law. However, some nontrivial physical points still remain open, which we revisit now. This discussion is also based on the fact that the well known Bekenstein-Hawking entropy can be expressed as being proportional to the event horizon area divided by the square of the Planck length. Keywords: black holes; nonadditive entropies; thermodynamics; complex systems 1. Introduction We frequently verify perplexity by various authors that the entropy of a black hole appears to be proportional to its area whereas it was expected to be proportional to its volume. Such perplexity is tacitly or explicitly based on a sort of belief (i) that a black hole is a three-dimensional system, and (ii) that its thermodynamic entropy is to be understood as its Boltzmann-Gibbs (BG) one.
    [Show full text]