Tsallis Statistics of the Magnetic Field in the Heliosheath Lf Burlaga,1 Af Vin

Total Page:16

File Type:pdf, Size:1020Kb

Tsallis Statistics of the Magnetic Field in the Heliosheath Lf Burlaga,1 Af Vin The Astrophysical Journal, 644:L83–L86, 2006 June 10 ᭧ 2006. The American Astronomical Society. All rights reserved. Printed in U.S.A. TSALLIS STATISTICS OF THE MAGNETIC FIELD IN THE HELIOSHEATH L. F. Burlaga,1 A. F. Vin˜ as,1 N. F. Ness,2 and M. H. Acun˜a3 Received 2006 March 27; accepted 2006 May 3; published 2006 May 25 ABSTRACT The spacecraft Voyager 1 crossed the termination shock on 2004 December 16 at a distance of 94 AU from the Sun, and it has been moving through the heliosheath toward the interstellar medium since then. The distributions of magnetic field strength B observed in the heliosheath are Gaussian over a wide range of scales, yet the measured profile appears to be filamentary with occasional large jumps in the magnetic field strengthB(t) . All of the probability distributions of changes in B,dBn { B(t ϩ t) Ϫ B(t) on scales t from 1 to 128 days, can be fit with the symmetric Tsallis distribution of nonextensive statistical mechanics. At scales ≥32 days, the distributions are Gaussian, but on scales from 1 through 16 days the probability distribution functions have non-Gaussian tails, suggesting that the inner heliosheath is not in statistical equilibrium on scales from 1 to 16 days. Subject headings: magnetic fields — solar wind — turbulence — waves 1. INTRODUCTION is expressed by the form of u in the distribution function. In The solar wind moves supersonically beyond a few solar Tsallis statistical mechanics, the probability that a system is in a state with energy u is proportional toexp (Ϫ u) . In classical radii from the Sun to a “termination shock” at ≈94 AU, carrying q b statistical mechanics, the probability that a system is in a state solar magnetic fields with it. The heliosheath is a subsonic Ϫ region extending from the termination shock to the interstellar with energy u is proportional toexp ( bu) , the Boltzmann-Gibbs distribution. For example, in a gas of molecules, the probability medium (Axford 1972; Parker 1963; Hundhausen 1972). Voy- p ager 1 entered the heliosheath by crossing the termination of finding molecules with energyu E is proportional to exp (ϪE/kT), where T is the temperature and k is Boltzmann’s shock on 2004 December 16 (Stone et al. 2005; Gurnett & p Kurth 2005; Decker et al. 2005; Burlaga et al. 2005), and it constant (b 1/kT ). has been moving through the inner heliosheath since that time. This Letter demonstrates that the PDFs of the fluctuations Magnetic field observations have been made by the MAG ex- of B observed in the heliosheath from DOY 1 through 256, periment on Voyager 1 (Behannon et al. 1977) from the launch 2005 on scales from 1 to 128 days can be described by the in 1977 to the present. Tsallis distribution. We show that the PDF of the heliosheath A generalization of Boltzmann-Gibbs statistical mechanics magnetic field strength is Gaussian (consistent with statistical 116 days, but it is not in Gaussian at was introduced by Tsallis (1988, 2004a) and has been applied equilibrium) at scales scales between 1 and 4 days. The observation of the Tsallis to many complex physical systems that are not in thermal equi- distribution of fluctuations of B in the heliosheath and the mea- librium (Tsallis 2004b; Tsallis & Brigatti 2004; Gell-Mann & surement of the variation of the parameters of this distribution Tsallis 2004). The statistical mechanics of Tsallis introduces a with scale describe the multiscale state of magnetic strength new entropy function, which is appropriate for nonequilibrium, fluctuations in the heliosheath and place strong constraints on nonextensive systems in which phase space has a hierarchical any model of these fluctuations. multiscale structure (rather than being uniformly occupied as it is for an equilibrium Gaussian distribution). By extremizing the entropy function subject to a normalization condition and 2. OBSERVATIONS OF DAILY AVERAGES OF MAGNETIC FIELD a constraint on energy, Tsallis derived a probability distribution STRENGTH IN HELIOSHEATH AND THEIR DISTRIBUTIONS function that describes many types of nonlinear physical sys- Daily averages of the magnetic field strength B as a function tems, e.g., “chaotic systems,” turbulence, cosmic-ray spectra, of time (DOY, day of year) from DOY 1 to 256, 2005 are plasma physics, and solar wind physics. shown in the top panel of Figure 1. The error bars on B are The probability distribution function (PDF) derived by Tsal- ≈0.015 nT, the estimated 1 j error due to both systematic and lis is proportional to the “q-exponential function” random effects. Two components of B are determined by rolls about the z-axis (directed toward the Sun) at the times indicated Ϫ { ϩ Ϫ 1/(1Ϫq) expq ( bu) [1 (q 1)bu] , (1) by the asterisks in the top panel of Figure 1. The third com- ponent (along the z-axis) is estimated at the time of a roll by p where b and q are constants. This function reduces to an ex- assumingABz S 0 in an interval of 52 days centered on a roll. ponential in the variable u when q approaches 1, and it ap- Corrections for varying contributions to B from the spacecraft proaches a power law in u when u is large. When q approaches and instrument are made between the rolls, using the infor- p 22Ϫ Ϫ 1 andu y ,expq ( bu)exp( approachesby ) , which is mation from both the primary and secondary instruments (Be- proportional to a Gaussian distribution. hannon et al. 1977). The physics for a particular application of the theory of Tsallis The fluctuations in daily averages of B from DOY 1 to 256 have a distribution that is approximately Gaussian with a mean 1 Laboratory for Solar and Space Physics, Code 612.2, NASA Goddard ABS p 0.11 nT, as shown in the bottom panel of Figure 1. The Space Flight Center, Greenbelt, MD 20771; [email protected]. p ϩ 2 solid curve is a fit to the Gaussian function B B0 Institute for Astrophysics and Computational Sciences, Catholic University 1/2Ϫ Ϫ 2 of America, Washington, DC 20064; [email protected]. {A/[w(p/2) ]} exp { 2[(B ABS)/w]}; the dotted curves are 3 Planetary Magnetospheres Laboratory, Code 695, NASA Goddard Space the 95% confidence bands based on a t-test computed by the Flight Center, Greenbelt, MD 20771; [email protected]. plotting program Origin. L83 L84 BURLAGA ET AL. Vol. 644 Fig. 1.—Top: Daily averages of the magnetic field strength B vs. time measured by Voyager 1 from DOY 1 to 256, 2005 in the heliosheath. Bottom: Distribution of B (filled squares) as well as a fit to a Gaussian distribution function (solid curve) and the 95% confidence bands (dotted curves). Fig. 2.—The symbols are PDFs of relative changes in the magnetic field strength measured by Voyager 1 in the heliosheath from DOY 1 to 256, 2005 The quality of the fit is given by the coefficient of determi- on scales from 1 to 128 days. The solid curves are fits of the data to the Tsallis nationR 2 p 0.89 , and the parameters derived from the fit are distribution function. p p ע p ע p Ϫ B0 1.8 6.0 nT,ABS 0.099 0.004 nT, w 2 j ע p ע 0.10 0.02 nT, andA 5.6 1.6 . bins that were used in constructing the histograms; each his- Gaussian distributions of B were also found at smaller scales togram is plotted a factor of 100 above the one below it for and in different intervals in the heliosheath (Burlaga et al. clarity. The error bars, the short horizontal segments above and 2005). However, the magnetic field strength fluctuations in the below the observed points in Figure 2, were computed as the top panel of Figure 1 have a filamentary structure rather than square root of the countsCi in bin i; thus, the error bars are 1/2 ע simple Gaussian fluctuations. For example, there are relatively plotted atCii(C ) with a shift of 100 for adjacent distri- large jumps at steps and peaks inB(t) , which we shall describe butions introduced for clarity. quantitatively and examine in some detail. These jumps and Using the Levenberg (1944) and Marquardt (1963) method, other changes related to various heliosheath features can be we fit these observed PDFs with the symmetric Tsallis described quantitatively by the probability distributions, as dis- distribution cussed in § 3. They are also associated with the multifractal structure observed in the interval under consideration (L. F. Ϫ Ϫ A exp [Ϫ (dBn)]22{ A [1 ϩ (q Ϫ 1) dBn ]1/(q 1) , (2) Burlaga et al. 2006, in preparation). qqb q qb q where q is the “entropic index” or “nonextensivity” factor, 3. MULTISCALE PROBABILITY DISTRIBUTION FUNCTIONS Ϫϱ ! q ≤ 3, which depends on scale (Tsallis 1988; Burlaga & From the daily averages of observations of B in the helio- Vin˜as 2004, 2005).ABqq and refer to a Tsallis distribution with sheath (Fig. 1), we computed a set of PDFs (histograms) de- a particular q; they too depend on scale. Forq p 1 , the Tsallis { ϩ Ϫ scribing increments in B,ondBn(tin; t ) B(t it n) B(t i) distribution reduces to the Boltzmann-Gibbs (Gaussian) distri- p n p scalestn 2 days, wheren 0 , 1, 2, 3, 4, 5, 6, and 7. Thus, bution. Burlaga & Vin˜as (2004, 2005, 2006) showed that this we consider scales from 1 to 128 days. The plasma is being Tsallis distribution (eq. [2]) provides good fits to PDFs of the convected past the spacecraft, but the speed is not measured magnetic field strength B observed in the supersonic solar wind directly; hence, we cannot determine the length scale associated on scales from 1 to 128 days from 1 to ≈87 AU.
Recommended publications
  • On Escort Distributions, Q-Gaussians and Fisher Information Jean-François Bercher
    On escort distributions, q-gaussians and Fisher information Jean-François Bercher To cite this version: Jean-François Bercher. On escort distributions, q-gaussians and Fisher information. 30th International Workshop on Bayesian Inference and Maximum Entropy Methods in Science and Engineering, Jul 2010, Chamonix, France. pp.208-215, 10.1063/1.3573618. hal-00766571 HAL Id: hal-00766571 https://hal-upec-upem.archives-ouvertes.fr/hal-00766571 Submitted on 18 Dec 2012 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. On escort distributions, q-gaussians and Fisher information J.-F. Bercher Laboratoire d’Informatique Gaspard Monge, Université Paris-Est, ESIEE, 5 bd Descartes, 77454 Marne-la-Vallée Cedex 2, France Abstract. Escort distributions are a simple one parameter deformation of an original distribution p. In Tsallis extended thermostatistics, the escort-averages, defined with respect to an escort dis- tribution, have revealed useful in order to obtain analytical results and variational equations, with in particular the equilibrium distributions obtained as maxima of Rényi-Tsallis entropy subject to constraints in the form of a q-average. A central example is the q-gaussian, which is a generalization of the standard gaussian distribution.
    [Show full text]
  • Arxiv:1707.03526V1 [Cond-Mat.Stat-Mech] 12 Jul 2017 Eq
    Generalized Ensemble Theory with Non-extensive Statistics Ke-Ming Shen,∗ Ben-Wei Zhang,y and En-Ke Wang Key Laboratory of Quark & Lepton Physics (MOE) and Institute of Particle Physics, Central China Normal University, Wuhan 430079, China (Dated: October 17, 2018) The non-extensive canonical ensemble theory is reconsidered with the method of Lagrange multipliers by maximizing Tsallis entropy, with the constraint that the normalized term of P q Tsallis' q−average of physical quantities, the sum pj , is independent of the probability pi for Tsallis parameter q. The self-referential problem in the deduced probability and thermal quantities in non-extensive statistics is thus avoided, and thermodynamical relationships are obtained in a consistent and natural way. We also extend the study to the non-extensive grand canonical ensemble theory and obtain the q-deformed Bose-Einstein distribution as well as the q-deformed Fermi-Dirac distribution. The theory is further applied to the general- ized Planck law to demonstrate the distinct behaviors of the various generalized q-distribution functions discussed in literature. I. INTRODUCTION In the last thirty years the non-extensive statistical mechanics, based on Tsallis entropy [1,2] and the corresponding deformed exponential function, has been developed and attracted a lot of attentions with a large amount of applications in rather diversified fields [3]. Tsallis non- extensive statistical mechanics is a generalization of the common Boltzmann-Gibbs (BG) statistical PW mechanics by postulating a generalized entropy of the classical one, S = −k i=1 pi ln pi: W X q Sq = −k pi lnq pi ; (1) i=1 where k is a positive constant and denotes Boltzmann constant in BG statistical mechanics.
    [Show full text]
  • Beyond Behavioural Economics
    Beyond Behavioural Economics A Process Realist Perspective Mathematics of Behavioural Economics and Finance James Juniper November, 2015 Ultimate scientific explanation is attractive and important, but financial engineering cannot wait for full explanation. That is, it is legitimate to strive towards a second best: a “descriptive phenomenology” that is organized tightly enough to bring a degree of order and understanding. Benoit B. Mandelbrot (1997) Overview • Keynes on DM under uncertainty & role of conventions • Ontological considerations • Objective not just subjective • Applications: • Bounded sub-additivity (Tversky & Wakker) • Equivalent approach—fuzzy measure theory (approximation algebras) • Concept lattices and “framing” • Tsallis arithmetic: q-generalized binomial approximation […] civilization is a thin and precarious crust erected by the personality and will of a very few and only maintained by rules and conventions skillfully put across and guilefully preserved” (Keynes, X: 447) The “love of money as a possession” is a “somewhat disgusting morbidity, one of those semi-criminal, semi-pathological propensities which one hands over with a shudder to the specialists in mental disease (IX: 329) For we shall enquire more curiously than is safe today into the true character of this ‘purposiveness’ with which in varying degrees Nature has endowed almost all of us. For purposiveness means that we are more concerned with the remote future results of our actions than with their own quality or their immediate effects on the environment. The ‘purposive’ man is always trying to secure a spurious and delusive immortality for his acts by pushing his interest in them forward into time. He does not love his cat, but his cat’s kittens; nor, in truth, the kittens, but only the kittens’ kittens, and so on forward to the end of catdom.
    [Show full text]
  • Fractal Structure and Non-Extensive Statistics
    entropy Article Fractal Structure and Non-Extensive Statistics Airton Deppman 1,* ID , Tobias Frederico 2, Eugenio Megías 3,4 and Debora P. Menezes 5 1 Instituto de Física, Universidade de São Paulo, Rua do Matão Travessa R Nr.187, Cidade Universitária, CEP 05508-090 São Paulo, Brazil 2 Instituto Tecnológico da Aeronáutica, 12228-900 São José dos Campos, Brazil; [email protected] 3 Departamento de Física Teórica, Universidad del País Vasco UPV/EHU, Apartado 644, 48080 Bilbao, Spain; [email protected] 4 Departamento de Física Atómica, Molecular y Nuclear and Instituto Carlos I de Física Teórica y Computacional, Universidad de Granada, Avenida de Fuente Nueva s/n, 18071 Granada, Spain 5 Departamento de Física, CFM, Universidade Federal de Santa Catarina, CP 476, CEP 88040-900 Florianópolis, Brazil; [email protected] * Correspondence: [email protected] Received: 27 June 2018; Accepted: 19 August 2018; Published: 24 August 2018 Abstract: The role played by non-extensive thermodynamics in physical systems has been under intense debate for the last decades. With many applications in several areas, the Tsallis statistics have been discussed in detail in many works and triggered an interesting discussion on the most deep meaning of entropy and its role in complex systems. Some possible mechanisms that could give rise to non-extensive statistics have been formulated over the last several years, in particular a fractal structure in thermodynamic functions was recently proposed as a possible origin for non-extensive statistics in physical systems. In the present work, we investigate the properties of such fractal thermodynamical system and propose a diagrammatic method for calculations of relevant quantities related to such a system.
    [Show full text]
  • Thermodynamics from First Principles: Correlations and Nonextensivity
    Thermodynamics from first principles: correlations and nonextensivity S. N. Saadatmand,1, ∗ Tim Gould,2 E. G. Cavalcanti,3 and J. A. Vaccaro1 1Centre for Quantum Dynamics, Griffith University, Nathan, QLD 4111, Australia. 2Qld Micro- and Nanotechnology Centre, Griffith University, Nathan, QLD 4111, Australia. 3Centre for Quantum Dynamics, Griffith University, Gold Coast, QLD 4222, Australia. (Dated: March 17, 2020) The standard formulation of thermostatistics, being based on the Boltzmann-Gibbs distribution and logarithmic Shannon entropy, describes idealized uncorrelated systems with extensive energies and short-range interactions. In this letter, we use the fundamental principles of ergodicity (via Liouville's theorem), the self-similarity of correlations, and the existence of the thermodynamic limit to derive generalized forms of the equilibrium distribution for long-range-interacting systems. Significantly, our formalism provides a justification for the well-studied nonextensive thermostatistics characterized by the Tsallis distribution, which it includes as a special case. We also give the complementary maximum entropy derivation of the same distributions by constrained maximization of the Boltzmann-Gibbs-Shannon entropy. The consistency between the ergodic and maximum entropy approaches clarifies the use of the latter in the study of correlations and nonextensive thermodynamics. Introduction. The ability to describe the statistical tion [30] or the structure of the microstates [27]. Another state of a macroscopic system is central to many ar- widely-used approach is to generalize the MaxEnt prin- eas of physics [1{4]. In thermostatistics, the statistical ciple to apply to a different entropy functional in place of BGS state of a system of N particles in equilibrium is de- S (fwzg).
    [Show full text]
  • Tsallis Entropy and the Vlasov-Poisson Equations
    Brazilian Journal of Physics, vol. 29, no. 1, March, 1999 79 Tsallis Entropy and the Vlasov-Poisson Equations 1;2;3 3;4 A. R. Plastino and A. Plastino 1 Faculty of Astronomy and Geophysics, National University La Plata, C. C. 727, 1900 La Plata, Argentina 2 CentroBrasileiro de Pesquisas F sicas CBPF Rua Xavier Sigaud 150, Rio de Janeiro, Brazil 3 Argentine Research Agency CONICET. 4 Physics Department, National University La Plata, C.C. 727, 1900 La Plata, Argentina Received 07 Decemb er, 1998 We revisit Tsallis Maximum Entropy Solutions to the Vlasov-Poisson Equation describing gravita- tional N -b o dy systems. We review their main characteristics and discuss their relationshi p with other applicatio ns of Tsallis statistics to systems with long range interactions. In the following con- siderations we shall b e dealing with a D -dimensional space so as to b e in a p osition to investigate p ossible dimensional dep endences of Tsallis' parameter q . The particular and imp ortant case of the Schuster solution is studied in detail, and the p ertinent Tsallis parameter q is given as a function of the space dimension. In the sp ecial case of three dimensional space we recover the value q =7=9, that has already app eared in many applications of Tsallis' formalism involving long range forces. I Intro duction R q 1 [f x] dx S = 1 q q 1 where q is a real parameter characterizing the entropy Foravarietyofphysical reasons, muchwork has b een functional S , and f x is a probability distribution de- devoted recently to the exploration of alternativeor q N ned for x 2 R .
    [Show full text]
  • Physica a X-Ray Binary Systems and Nonextensivity
    Physica A 389 (2010) 854–858 Contents lists available at ScienceDirect Physica A journal homepage: www.elsevier.com/locate/physa X-ray binary systems and nonextensivity Marcelo A. Moret a,b,∗, Valter de Senna a, Gilney F. Zebende a,b, Pablo Vaveliuk a a Programa de Modelagem Computational - SENAI - CIMATEC 41650-010 Salvador, Bahia, Brazil b Departamento de Física, Universidade Estadual de Feira de Santana, CEP 44031-460, Feira de Santana, Bahia, Brazil article info a b s t r a c t Article history: We study the x-ray intensity variations obtained from the time series of 155 light curves Received 9 April 2009 of x-ray binary systems collected by the instrument All Sky Monitor on board the satel- Received in revised form 14 August 2009 lite Rossi X-Ray Timing Explorer. These intensity distributions are adequately fitted by q- Available online 20 October 2009 Gaussian distributions which maximize the Tsallis entropy and in turn satisfy a nonlinear Fokker–Planck equation, indicating their nonextensive and nonequilibrium behavior. From PACS: the values of the entropic index q obtained, we give a physical interpretation of the dynam- 05.10.Gg ics in x-ray binary systems based on the kinetic foundation of generalized thermostatistics. 05.90.+m 89.75.Da The present findings indicate that the binary systems display a nonextensive and turbulent 97.80.-d behavior. ' 2009 Elsevier B.V. All rights reserved. Keywords: Astrophysical sources Tsallis statistics In recent years, the generalized thermostatistical formalism (GTS) proposed by Tsallis [1] has received increasing attention due its success in the description of certain phenomena exhibiting atypical thermodynamical features.
    [Show full text]
  • Nonextensive Statistical Mechanics and High Energy Physics
    EPJ Web of Conferences 71, 00132 (2014) DOI: 10.1051/epjconf/20147100132 C Owned by the authors, published by EDP Sciences, 2014 Nonextensive statistical mechanics and high energy physics Constantino Tsallis1,2,a and Zochil Gonzalez Arenas1,b 1Centro Brasileiro de Pesquisas Fisicas, and National Institute of Science and Technology for Complex Systems - Rua Xavier Sigaud 150, 22290-180 Rio de Janeiro-RJ, Brazil 2Santa Fe Institute - 1399 Hyde Park Road, Santa Fe, NM 87501, USA Abstract. The use of the celebrated Boltzmann-Gibbs entropy and statistical mechanics is justified for ergodic-like systems. In contrast, complex systems typically require more powerful theories. We will provide a brief introduction to nonadditive entropies (char- acterized by indices like q,which,intheq → 1 limit, recovers the standard Boltzmann- Gibbs entropy) and associated nonextensive statistical mechanics. We then present some recent applications to systems such as high-energy collisions, black holes and others. In addition to that, we clarify and illustrate the neat distinction that exists between Lévy distributions and q-exponential ones, a point which occasionally causes some confusion in the literature, very particularly in the LHC literature. 1 Introduction Boltzmann-Gibbs (BG) statistical mechanics constitutes one of the pillars of contemporary physics, together with Newtonian, quantum and relativistic mechanics, and Maxwell’s electromagnetism. Cen- tral points of any physical theory are when and why it works, and when and why it fails. This is so for any human intellectual construct, hence for the BG theory as well. We intend to briefly discuss here when we must (or must not) use the BG entropy W W 1 S = k p ln p = 1 , (1) BG i p i i=1 i i=1 where k is a constant (either taken equal to Boltzmann constant kB, or to unity).
    [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]
  • Non-Gaussian Stochastic Models and Their Applications in Econophysics
    UNIVERSITÀ DEGLI STUDI DI PAVIA DOTTORATO DI RICERCA IN FISICA – XX CICLO Non-Gaussian Stochastic Models and Their Applications in Econophysics Enrica Vera Cisana Tesi per il conseguimento del titolo Università degli Dipartimento di Fisica Istituto Nazionale di Studi di Pavia Nucleare e Teorica Fisica Nucleare DOTTORATO DI RICERCA IN FISICA – XX CICLO Non-Gaussian Stochastic Models and Their Applications in Econophysics dissertation submitted by Enrica Vera Cisana to obtain the degree of DOTTORE DI RICERCA IN FISICA Supervisor: Prof. Guido Montagna Referee: Prof. Rosario N. Mantegna Cover Left: Traders at work in the New York Stock Exchange (NYSE). Right, top: Simulated paths of the time evolution of volatility process in the Heston model. More details can be found in Fig. 4.1 of this work. Right, bottom: Pictorial representation of the frontier of two-dimensional Brownian motion. Non-Gaussian Stochastic Models and Their Applications in Econophysics Enrica Vera Cisana PhD thesis – University of Pavia Printed in Pavia, Italy, November 2007 ISBN 978-88-95767-06-2 To Maria To Silvana E intanto il tempo passa e tu non passi mai... Se potessi far tornare indietro il mondo farei tornare poi senz'altro te per un attimo di eterno e di profondo in cui tutto sembra e niente c' `e Negramaro, Estate & Immenso Contents Introduction 1 1 Why non-Gaussian models for market dynamics? 5 1.1 Financial markets: complex systems and stochastic dynamics . 6 1.2 The Black and Scholes-Merton paradigm . 7 1.3 Stylized facts of real market dynamics . 9 1.3.1 Non-Gaussian nature of log-return distribution .
    [Show full text]
  • Useful Mathematical Formulae
    Appendix A Useful Mathematical Formulae x1−q − 1 ln x ≡ (x > 0, q ∈ R)(A.1) q 1 − q 1−q lnq x = x ln2−q x (x > 0; ∀q)(A.2) lnq (1/x) + ln2−q x = 0(x > 0; ∀q)(A.3) q q lnq x + ln(1/q)(1/x ) = 0(x > 0; ∀q)(A.4) ⎧ < < − / − , ⎪ 0ifq 1 and x 1 (1 q) ⎪ ⎪ ⎪ 1 ⎨⎪ + − 1−q < ≥− / − , 1 [1 (1 q) x] if q 1 and x 1 (1 q) x 1−q e ≡ [1 + (1 − q) x]+ ≡ q ⎪ ⎪ ex if q = 1(∀x) , ⎪ ⎪ ⎩ 1 [1 + (1 − q) x] 1−q if q > 1 and x < 1/(q − 1) . (A.5) x −x = ∀ eq e2−q 1(q)(A.6) x q −qx = ∀ (eq ) e(1/q) 1(q)(A.7) x+y+(1−q)xy = x y ∀ eq eq eq ( q)(A.8) 329 330 Appendix A Useful Mathematical Formulae x ⊕q y ≡ x + y + (1 − q) xy (A.9) For x ≥ 0 and y ≥ 0: ⎧ −q −q ⎪ 0ifq < 1andx 1 + y1 < 1 , ⎪ ⎪ ⎪ 1 ⎪ 1−q 1−q 1−q 1−q 1−q 1 ⎨ [x + y − 1] if q < 1andx + y ≥ 1 , 1−q 1−q 1−q x ⊗ y ≡ [x + y − 1]+ ≡ q ⎪ ⎪ xy if q = 1 ∀(x, y) , ⎪ ⎪ ⎩ − − 1 − − [x 1 q + y1 q − 1] 1−q if q > 1andx 1 q + y1 q > 1 . (A.10) 1 − x ⊗q y = [1 + (1 − q)(lnq x + lnq y)] 1 q (A.11) ⊕ x q y = x y ∀ eq eq eq ( q) (A.12) x+y = x ⊗ y ∀ eq eq q eq ( q) (A.13) d ln x 1 q = (x > 0; ∀q) (A.14) dx xq dex q = (ex )q (∀q) (A.15) dx q x q = qx ∀ (eq ) e2−(1/q) ( q) (A.16) x a = ax ∀ (eq ) e1−(1−q)/a ( q) (A.17) / − 2 − x 1/(q−1) 1 − b (q 1) a b b a− x x e = x q−1 e (b > 0; q > 1) (A.18) q q − 1 q Appendix A Useful Mathematical Formulae 331 1 1 3 1 2 1 ex = ex [1 − (1 − q)x2 + (1 − q)2x3(1 + x) − (1 − q)3x4(1 + x + x2) q 2 3 8 4 3 12 1 65 5 5 + (1 − q)4x5(1 + x + x2 + x3) 5 72 24 384 1 11 17 1 1 − (1 − q)5x6(1 + x + x2 + x3 + x4) + ...](q → 1; ∀x) 6 10 48 24 640 (A.19) 1 1 1 ln x = ln x [1 + (1 − q)lnx + (1 − q)2 ln2 x + (1 − q)3 ln3 x q 2 6 24 1 1 + (1 − q)4 ln4 x + (1 − q)5 ln5 x + ...](q → 1; x > 0) 120 720 (A.20) x ⊗q y = xy 1 − (1 − q)(ln x)(ln y) 1 + (1 − q)2 [(ln2 x)(ln y) + (ln x)(ln2 y) + (ln2 x)(ln2 y)] 2 1 − (1 − q)3 [2(ln3 x)(ln y) + 9(ln2 x)(ln2 y) + 2(ln x)(ln3 y) 12 +6(ln3 x)(ln2 y) + 6(ln2 x)(ln3 y) + 2(ln3 x)(ln3 y)] 1 + (1 − q)4 [(ln4 x)(ln y) + 14(ln3 x)(ln2 y) 24 +14(ln2 x)(ln3 y) + (ln x)(ln4 y) + 4 2 + 3 3 + 2 4 7(ln x)(ln y) 24(ln x)(ln y) 7(ln x)(ln y) +6(ln4 x)(ln3 y) + 6(ln3 x)(ln4 y) + (ln4 x)(ln4 y)] + ..
    [Show full text]
  • Asymmetric Tsallis Distributions for Modelling Financial Market Dynamics
    Asymmetric Tsallis distributions for modelling financial market dynamics Sandhya Devi† 509 6th Ave S, Edmonds, WA 98020, United States of America Email: [email protected] Abstract: Financial markets are highly non-linear and non-equilibrium systems. Earlier works have suggested that the behavior of market returns can be well described within the framework of non-extensive Tsallis statistics or superstatistics. For small time scales (delays), a good fit to the distributions of stock returns is obtained with q-Gaussian distributions, which can be derived either from Tsallis statistics or superstatistics. These distributions are symmetric. However, as the time lag increases, the distributions become increasingly non-symmetric. In this work, we address this problem by considering the data distribution as a linear combination of two independent normalized distributions – one for negative returns and one for positive returns. Each of these two independent distributions are half q-Gaussians with different non-extensivity parameter q and temperature parameter beta. Using this model, we investigate the behavior of stock market returns over time scales from 1 to 80 days. The data covers both the .com bubble and the 2008 crash periods. These investigations show that for all the time lags, the fits to the data distributions are better using asymmetric distributions than symmetric q-Gaussian distributions. The behaviors of the q parameter are quite different for positive and negative returns. For positive returns, q approaches a constant value of 1 after a certain lag, indicating the distributions have reached equilibrium. On the other hand, for negative returns, the q values do not reach a stationary value over the time scales studied.
    [Show full text]