January 2015 BI-TP 2015/09 The Legacy of Rolf Hagedorn: Statistical Bootstrap and Ultimate Temperature 5 ∗ 1 0 2 n a Krzysztof Redlich1 and Helmut Satz2 J 9 2 1 InstituteofTheoretical Physics,UniversityofWrocław,Poland ] h 2 Fakulta¨tfur Physik,Universita¨tBielefeld, Germany p - p e h [ Abstract: 1 v 3 Inthelatterhalfofthelastcentury,itbecameevidentthatthereexistsaneverincreasingnumber 2 5 of different states of the so-called elementary particles. The usual reductionist approach to 7 this problem was to search for a simpler infrastructure, culminating in the formulation of the 0 . quark model and quantum chromodynamics. In a complementary, completely novel approach, 1 0 Hagedorn suggested that the mass distribution of the produced particles follows a self-similar 5 composition pattern, predicting an unbounded number of states of increasing mass. He then 1 concluded that such a growth would lead to a limiting temperature for strongly interacting : v matter. We discuss the conceptual basis for this approach, its relation to critical behavior, and i X itssubsequentapplicationsin differentareas ofhighenergy physics. r a To appear in MeltingHadrons,BoilingQuarks, R. Hagedorn and J. Rafelski (Ed.), Springer ∗ Verlag 2015,open access. Aprophetisnotwithouthonour, butinhisowncountry. TheNew Testament,Mark6,4. 1 Rolf Hagedorn The development of physics is the achievement of physicists, of humans, persisting against oftenconsiderableodds. Eveninphysics,fashionratherthanfactisfrequentlywhatdetermines thejudgementand recognition. When Rolf Hagedorn carried out his main work, now quite generally recognized as truly pio- neering, much of the theoretical community not only ignored it, but even considered it to be nonsense. “Hagedorn istein Narr”, heis afool,was asummaryofmanyleadingGerman theo- ristsofhistime. Wheninthe1990’sthequestionwasbroughtupwhetherhecouldbeproposed for the Max Planck Medal, the highest honor of the German physics community, even then, whenhisachievementswerealreadyknownworld-wide,theanswerwasstill“proposed,yes...” At the time Hagedorn carried out his seminal research, much of theoretical physics was ideo- logically fixed on “causality, unitarity, Poincare´ invariance”: from these three concepts, from axiomatic quantum field theory, all that is relevant to physics must arise. Those who thought thatscienceshouldprogressinsteadbycomparisontoexperimentwerederogatedas“fittersand plotters”. Galileo was almost forgotten... Nevertheless, one of the great Austrian theorists of the time, Walter Thirring, himself probably closer to the fundamentalists, noted: “If you want to dosomethingreallynew,you first haveto haveanewidea”. Hagedorndid. He had a number of odds to overcome. He had studied physics in Go¨ttingen under Richard Becker, where he developed a life-long love for thermodynamics. When he took a position at CERN,shortlyaftercompletinghisdoctorate,itwastoperformcalulationsfortheplanningand constructionof theprotonsynchrotron. When that was finished, heshifted to thestudyof mul- tihadron production in proton-proton collisions and to modelling the results of these reactions. It took a while before various members of the community, including some of the CERN The- ory Division, were willing to accept the significance of his work. This was not made easier by Hagedorn’s stronglyfocussed region of interest, but eventuallyit became generally recognized that here was someone who, in this perhaps a little similar to John Bell, was developing truly novelideas whichat first sightseemed quitespecific, butwhich eventuallyturned outtohavea lastingimpactalsoon physicswell outsideitsregionsoforigin. We findthat RolfHagedorn’swork centers ontwo themes: the statistical bootstrap model, a self-similar scheme for the composition and decay of • hadronsand theirresonances;forHagedorn, thesewerethe“fireballs”. theapplicationoftheresultingresonancespectruminanidealgascontainingallpossible • hadrons and hadron resonances, and to the construction of hadron production models based onsuch athermalinput. We will address these topics in the first two sections, and then turn to their role both in the thermodynamics of strongly interacting matter and in the description of hadron production in 2 elementary as well as nuclear collisions. Our aim here is to provide a general overview of Hagedorn’sscientificachievements;otheraspectswillbecoveredinotherchaptersofthisbook. SomeofwhatwewillsaytranscendsHagedorn’slife. Butthen,toparaphraseShakespeare, we havecometopraiseHagedorn,nottoburyhim;wewanttoshowthathisideasarestillimportant and very muchalive. 2 The Statistical Bootstrap Around 1950, the world still seemed in order for those looking for the ultimate constituents of matter in the universe. Dalton’s atoms had been found to be not really atomos, indivisible; Rutherford’smodeloftheatomhadmadethemlittleplanetarysystems,withthenucleusasthe sun and the electrons as encircling planets. The nuclei in turn consisted of positively charged protonsandneutralneutronsastheessentialmasscarriers. Withanequalnumberofprotonsand electrons, the resulting atoms were electrically neutral, and the states obtained by considering thedifferent possiblenucleus compositionsreproduced the periodictable ofelements. So for a short time, the Greek dream of obtaining the entire complex world by combining three simple elementaryparticles in different ways seemed finally feasible: protons, neutrons and electrons were thebuildingblocks ofouruniverse. But therewere thosewhorediscoveredan oldproblem,first formulatedby theRomanphiloso- pher Lucretius: if your elementary particles, in our case the protons and neutrons, have a size and a mass, as both evidently did, it was natural to ask what they are made of. An obvious way to find out is to hit them against each other and look at the pieces. And it turned out that there werelotsoffragments, themoretheharderthecollision. But theywere notreally pieces, sincethedebrisfoundafteraproton-protoncollisionstillalsocontainedthetwoinitialprotons. Moreover, the additional fragments, mesons and baryons, were in almost all ways as elemen- tary as protons and neutrons. The study of such collisions was taken up by more and more laboratories and at ever higher collision energies. As a consequence, the number of different “elementary”particles grew byleaps and bounds,fromtens to twentiestohundreds. Thelatest compilationoftheParticleDataGroupcontainsoverathousand. What to do? One approach was in principle obvious: just as Dalton’s atoms could be con- structed from simpler, more elementary constituents, so one had to find a way of reproducing allthehadrons,theparticlesformedinstronginteractioncollisions,intermsoffewerandmore elementarybuildingblocks. Thisconceptuallystraight-forwardproblemwas,however,farfrom easy: asimplyadditivecompositionwasnotpossible. Nevertheless,inthelate1960’sthequark model appeared, in which three quarks and their antiquarks were found to produce in a non- Abelian composition all the observed states and more, predicted and found. Not much later, quantum chromodynamics (QCD) appeared as the quantum field theory governing strong in- teractions; moreover, it kept the quarks inherently confined, without individual existence: they occurred only as quark-antiquark pair (a meson) or as a three-quark state (a baryon). And wealth of subsequent experiments confirmed QCD as the basic theory of strong interactions. Theconventionalreductionistapproach hadtriumphedoncemore. Let us, however, return to the time when physics was confronted by all those elementary par- ticles, challenging its practioners to find a way out. At this point, in the mid 1960’s, Rolf 3 Hagedorncameupwithatrulynovelidea[1]-[5]. Hewasnotsomuchworriedaboutthespe- cificpropertiesoftheparticles. Hejustimaginedthataheavyparticlewassomehowcomposed out of lighterones, and theseagain in turn ofstilllighterones, and so on, until onereached the pion as the lightest hadron. And by combining heavy ones, you would get still heavier ones, again: and so on. The crucial input was that the composition law should be the same at each stage. Todaywecallthatself-similarity,andithadbeenaroundinvariousformsformanyyears. AparticularlyelegantformulationwaswrittenahundredyearsbeforeHagedornbytheEnglish mathematicianAugustusdeMorgan,thefirst presidentoftheLondonMathematicalSociety: Greatfleashavelittlefleasupontheirbackstobite’em, andlittlefleashavelesserstill,andsoadinfinitum. Andthegreatfleasthemselves,inturn,havegreaterfleastogoon, whiletheseagainhavegreaterstill,andgreaterstill,andsoon. Hagedorn proposed that “a fireball consists of fireballs, which in turn consist of fireballs, and so on...” The concept later reappeared in various forms in geometry; in 1915, it led to the celebratedtriangledevisedbythePolishmathematicianWacławSierpinski: “atriangleconsists of triangles, which in turn consist of triangles, and so on...”, in the words of Hagedorn. Still later,shortlyafterHagedorn’sproposal,theFrenchmathematicianBenoitMandelbroitinitiated thestudyofsuchfractalbehaviouras anewfield ofmathematics. TheSierpinskiTriangle Hagedorn had recalled a similar problem in number theory: how many ways are there of de- composinganintegerintointegers? Thiswassomethingalreadyadressedin1753byLeonhard Euler, and more than a century later by the mathematician E. Schro¨der in Germany. Finally G. H.HardyandS.RamanujaninEnglandprovidedanasymptoticsolution. Letushere,however, considerasimplified,easilysolvableversionoftheproblem[6],inwhichwecountallpossible different ordered arrangements p(n) ofan integern. So wehave 1=1 p(1)=1=2n 1 − 2=2 ,1+1 p(2)=2=2n 1 − 3=3, 2+1,1+2, 1+1+1 p(3)=4=2n 1 − 4=4, 3+1,1+3, 2+2, 2+1+1,1+2+1, 1+1+2, 1+1+1+1 p(4)=8=2n 1 − and so on. In otherwords, thereare 1 p(n)=2n 1 = enln2 (1) − 2 4 ways of partitioning an integer n into ordered partitions: p(n) grows exponentially in n. In this particular case, the solution could be found simply by induction. But there is another way of getting it, more in line with Hagedorn’s thinking: “large integers consist of smaller integers, which in turn consist of still smaller integers, and so on...” This can be formulated as an equation, r (n)=d (n 1)+ (cid:229) n 1 (cid:213) k r (n) d (S n n). (2) i i i − k! − k=2 i=1 Itisquiteevidentherethattheformofthepartitionnumberr (n)isdeterminedbyaconvolution of many similar partitions of smaller n. The solution of the equation is in fact just the number ofpartitionsofn thatwe hadobtainedabove, r (n)=z p(n) (3) up to a normalization constant of order unity (for the present case, it turns out that z 1.25). ≃ For Hagedorn, eq. (2) expressed the idea that the structure of r (n) was determined by the structureofr (n)–wenowcallthisself-similar. HeinsteadthoughtofthelegendaryBaronvon Mu¨nchhausen, who had extracted himself from a swamp by pulling on his own bootstraps. So forhim,eq. (2)becameabootstrapequation. The problem Hagedorn had in mind was, of course, considerably more complex. His heavy resonancewasnotsimplyasumoflighteronesatrest,butitwasasystemoflighterresonances in motion,with therequirementthat thetotal energy ofthissystemadded up to themassofthe heavy one. And similarly, the masses of the lighter ones were the result of still ligher ones in motion. Thebootstrapequationforsuch asituationbecomes r (m,V )=d (m m ) +(cid:229) 1 V0 N−1 (cid:213) N [dm r (m) d3p ] d 4(S p p), (4) 0 − 0 N!(cid:20)(2p )3(cid:21) Z i i i i i− N i=1 where the first term corresponds to thecase of just one lightestpossibleparticle, a “pion”. The factorV ,theso-calledcompositionvolume,specifiedthesizeoftheoverallsystem,anintrinsic 0 fireball size. Since the mass of any resonance in the composition chain is thus determined by the sum over phase spaces containing lighter ones, whose mass is specified in the same way, Hagedorn called thisformofbootstrap“statistical”. After a number of numerical attempts by others, W. Nahm [7] solved the statistical bootstrap equationanalytically,obtaining r (m,V )=const.m 3exp m/T . (5) 0 − H { } So even though the partitioning now was not just additive in masses, but included the kinetic energy of themovingconstituents,the increase was again exponentialin mass. The coefficient oftheincrease, T 1,is determinedbytheequation H− V T3 0 H(m /T )2K (m /T )=2ln2 1, (6) 2p 2 0 H 2 0 H − in terms of two parameters V and m . Hagedorn assumed that the composition volume V , 0 0 0 specifying the intrinsic range of strong interactions, was determined by the inverse pion mass 5 asscale,V0 (4p /3)m−p 3. ThisleadstoascalefactorTH 150MeV.Itshouldbeemphasized, ≃ ≃ however,thatthisisjustonepossibleway to proceed. In thelimitm 0,eq. (6)gives 0 → TH =[p 2(2ln2−1)]1/3V0−1/3 ≃1/rh, (7) whereV =(4p /3)r3 and r denotes the range of strong interactions. With r 1 fm, we thus 0 h h h ≃ have T 200 MeV. From this it is evident that the exponential increase persists also in the H ≃ chirallimitmp 0andisinfactonlyweaklydependentonm0,providedthestronginteraction → scaleV iskeptfixed. 0 Theweightsr (m)determinethecompositionaswellasthedecayof“resonances”, offireballs. The basis of the entire formalism, the self-similarity postulate – here in the form of the statis- tical bootstrap condition – results in an unending sequence of ever-heavier fireballs and in an exponentiallygrowingnumberofdifferentstatesofagivenmassm. Before we turn to the implications of such a pattern in thermodynamics, we note that not long afterHagedorn’sseminalpaper,itwasfoundthataratherdifferentapproach,thedualresonance model [8–10] led to very much the same exponential increase in the number of states. In this model, any scattering amplitude, from an initial two to a final n hadrons, was assumed to be determinedbytheresonancepolesinthedifferentkinematicchannels. Thisresultedstructurally again in a partition problem of the same type, and again the solution was that the number of possible resonance states of mass m must grow exponentially in m, with an inverse scale factor of the same size as obtained above, some 200 MeV. Needless to say, this unexpected supportfromtheforefrontoftheoreticalhadrondynamicsconsiderablyenhancedtheinterestin Hagedorn’swork. 3 The Limiting Temperature of Hadronic Matter Consider a relativistic ideal gas of identical neutral scalar particles of mass m contained in a 0 boxofvolumeV,assumingBoltzmannstatistics. Thegrandcanonicalpartitionfunctionofthis systemisgivenby N Z(T,V)=(cid:229) 1 V d3p exp p2+m2 /T , (8) N!(cid:20)(2p )3Z {−q 0 }(cid:21) N leadingto VTm2 m lnZ(T,V)= 0 K ( 0). (9) 2p 2 2 T Fortemperatures T m , theenergy densityofthesystembecomes 0 ≫ 1 ¶ lnZ(T,V) 3 e (T)= T4, (10) −V ¶ (1/T) ≃ p 2 theparticledensity 6 ¶ ln Z(T,V) 1 n(T)= T3, (11) ¶ V ≃ p 2 and so theaverageenergy perparticleis givenby w 3T. (12) ≃ The important feature to learn from these relations is that, in the case of an ideal gas of one speciesofelementaryparticles,anincreaseoftheenergyofthesystemhasthreeconsequences: it leadsto ahighertemperature, • moreconstituents,and • moreenergeticconstituents. • Ifwenowconsideraninteractinggasofsuchbasichadronsandpostulatethattheessentialform of the interaction is resonance formation, then we can approximate the interacting medium as a non-interacting gas of all possible resonance species [11,12]. The partition function of this resonancegas is lnZ(T,V)=(cid:229) VTm2i r (m )K (mi) (13) 2p 2 i 2 T i wherethesumbeginswiththestablegroundstatem andthenincludesthepossibleresonances 0 m ,i = 1,2,... with weights r (m) relative to m . Clearly the crucial question here is how i i 0 to specify r (m ), how many states there are of mass m . It is only at this point that hadron i i dynamics enters, and it is here that Hagedorn introduced the result obtained in his statistical bootstrapmodel. As we had seen above in eq. (5), the density of states then increases exponentially in m, with a coefficient T 1 determined by eq. (6) in terms of two parameters V and m . If we replace H− 0 0 the sum in the resonance gas partition function (13) by an integral and insert the exponentially growingmassspectrum(5), eq. (13)becomes VT m lnZ(T,V) dm m2r (m )K ( i) ≃ 2p 2Z i 2 T 3/2 T 1 1 V dm m 3/2exp m . (14) ∼ (cid:20)2p (cid:21) Z − {− (cid:20)T −T (cid:21)} H Evidently, the result is divergent for all T > T : in other words, T is the highest posssible H H temperature of hadronic matter. Moreover, if we compare such a system with the ideal gas of onlybasicparticles(a “pion”gas),wefind piongas resonancegas np e 3/4 nres e ∼ ∼ w p e 1/4 w res const. ∼ ∼ Here n denotes the average number density of constituents, w the average energy of a con- stituent. In contrast toto thepiongas,an increase ofenergy nowleadsto 7 afixed temperaturelimit,T T , H • → themomentaoftheconstituentsdo notcontinuetoincrease, and • moreand morespecies ofeverheavierparticles appear. • We thus obtain a new, non-kinetic way to use energy, increasing the number of species and their masses, not the momentum per particle. Temperature is a measure of the momentum of theconstituents,andifthatcannotcontinuetoincrease, thereisahighestpossible,a“limiting” temperatureforhadronicsystems. Hagedorn originally interpreted T as the ultimate temperature of strongly interacting matter. H It is clear today that T signals the transition from hadronic matter to a quark-gluon plasma. H Hadron physicsalonecanonlyspecifyitsinherentlimit;togobeyondthislimit,weneedmore information: weneed QCD. As seen ineq. (5), thesolutionofthestatisticalbootstrapequationhas thegeneral form r (m,V ) m aexp m/T , (15) 0 − H ∼ { } withsomeconstanta;theexactsolutionofeq.(4)byNahmgavea=3. Itispossible,however, to consider variations of the bootstrap model which lead to different a, but always retain the exponentialincreaseinm. WhiletheexponentialformmakesT theupperlimitofpermissible H temperatures, the power law coefficient a determines the behavior of the system at T = T . H For a = 3, the partition function (14) itself exists at that point, while the energy density as first derivativein temperature diverges there. This is what made Hagedorn conclude that T is H indeed the highest possible temperature of matter: it would require an infinite energy to reach it. Onlyafew yearslateritwas,however,pointedoutbyN.CabibboandG.Parisi[13]thatlarger a shifted the divergence at T =T to ever higher derivatives. In particular, for 4>a>3, the H energy density would remain finite at that point, shifting the divergence to the specific heat as nexthigherderivative. Suchcriticalbehaviorwasinfactquiteconventionalinthermodynamics: it signalled a phase transition leading to the onset of a new state of matter. By that time, the quark model and quantum chromodynamics as fundamental theory of strong interactions had appeared and suggestedthe existenceof a quark-gluonplasmaas the relevantstateof matterat extreme temperature or density. It was therefore natural to interpret the Hagedorn temperature T as the critical transition temperature from hadronic matter to such a plasma. This interpre- H tation is moreover corroborated by a calculation of the critical exponents [15] governing the singularbehavioroftheresonancegas thermodynamicsbasedon aspectrumoftheform(15). It shouldbenoted,however,thatin somesenseT didremain thehighestpossibletemperature H of matter as we know it. Our matter exists in the physical vacuum and is constructed out of fundamental building blocks which in turn have an independent existence in this vacuum. Our matter ultimately consists of and can be broken up into nucleons; we can isolate and study a single nucleon. The quark-gluon plasma, on the other hand, has its own ground state, distinct fromthephysicalvacuum,anditsconstituentscanexistonlyinadensemediumofotherquarks – wecan neverisolateand studyasinglequark. That does not mean, however, that quarks are eternally confined to a given part of space. Let us start with atomic matter and compress that to form nuclear matter, as it exists in heavy nu- clei. At thisstage,wehavenucleonsexistingin thephysicalvacuum. Each nucleonconsistsof 8 three quarks, and theyare confined to remain closeto each other; thereis no way to break up a given nucleon into its quark constituents. But if we now continue to compress, then eventually the nucleons will then penetrate each other, until we reach a dense medium of quarks. Now each quark finds in its immediate neighborhood many other quarks besides those which were with it in the nucleon stage. It is therefore no longerpossibleto partitionquarks into nucleons; the medium consists of unbound quarks, whose interaction becomes ever weaker with increas- ing density, approaching the limit of asymptotic freedom predicted by QCD. Any quark can now move freely throughout the medium: we have quark liberation through swarm formation. Wherever a quark goes, there are many other quarks nearby. The transition from atomic to quark matteris schematicallyillustratedinFig. 1. (a) (b) (c) Figure 1: Schematic view of matter for increasing density, from atomic (a) to nuclear (b) and then toquark matter(c). Wehavehereconsideredquark matterformationthroughthecompressionofcoldnuclearmat- ter. A similareffect is obtained if we heat a meson gas; with increasing temperature, collisions andpairproductionleadtoaneverdensermediumofmesons. AndaccordingtoHagedorn,also of ever heavier mesons of an increasing degeneracy. For Hagedorn, the fireballs where point- like, so that the overlap we had just noted simply does not occur. In the real world, however, they do have hadronicsize, so that theywill in fact interpenetrateand overlapbefore the diver- gence of the Hagedorn resonance gas occurs [16]. Hence now again there will be a transition from resonance gas to a quark-gluon plasma, now formed by the liberation of the quarks and gluonsmakinguptheresonances. Atthispoint,itseemsworthwhiletonoteanevenearlierapproachleadingtoalimitingtemper- ature for hadronic matter. More than a decade before Hagedorn, I. Ya. Pomeranchuk [17] had pointedoutthat acrucial feature ofhadronsis theirsize, and hencethedensityofanyhadronic medium is limited by volume restriction: each hadron must have its own volume to exist, and once the density reaches the dense packing limit, it’s the end for hadronic matter. This simply ledtoatemperaturelimit,andforanidealgasofpionsof1fmradius,theresultingtemperature was again around 200 MeV. Nevertheless, these early results remained largely unnoticed until thework ofHagedorn. Such geometric considerations do, however, lead even further. If hadrons are allowed to in- terpenetrate, to overlap, then percolation theory predicts two different states of matter [18,19]: hadronicmatter,consistingofisolatedhadronsorfinitehadronicclusters,andamediumformed 9 asaninfinitesizedclusterofoverlappinghadrons. Thetransitionfromonetotheothernowbe- comes agenuinecritical phenomenon,occurring at acriticalvalueofthehadron density. We thus conclude that the pioneering work of Rolf Hagedorn opened up the field of critical behavior in strong interaction physics, a field in which still today much is determined by his ideas. On a more theoretical level, the continuation of such studies was provided by finite temperature lattice QCD, and on the more experimental side, by resonance gas analyses of the hadron abundances in high energy collisions. In both cases, it was found that the observed behaviorwas essentiallythatpredicted by Hagedorn’sideas. 4 Resonance gas and QCD thermodynamics With the formulation of Quantum Chromodynamics (QCD) as a theoretical framework for the strong interaction force among elementary particles it became clear that the appearance of the ultimate Hagedorn temperature T signals indeed the transition from hadronic phase to a new H phaseofstronglyinteractingmatter,thequark-gluonplasma(QGP)[20]. AsQCDisanasymp- totically free theory, the interaction between quarks and gluons vanishes logarithmically with increasingtemperature,thusatveryhightemperaturestheQGPeffectivelybehaveslikeanideal gasofquarks and gluons. Today we have detailed information, obtained from numerical calculations in the framework of finite temperature lattice Quantum Chromodynamics [21–25], about the thermodynamics of hot and dense matter. We know the transition temperature to the QGP and the tempera- ture dependence of basic bulk thermodynamic observables such as the energy density and the pressure [25,26]. We also begin to have results on fluctuations and correlations of conserved charges [27–29]. The recent increase in numerical accuracy of lattice QCD calculations and their extrapolation tothecontinuumlimitmakeitpossibletoconfrontthefundamentalresultsofQCDwithHage- dorn’sconcepts [1,5], which providea theoreticalscenario forthethermodynamicsofstrongly interactinghadronicmatter[29–32]. In particular, the equation of state calculated on the lattice at vanishingand finite chemical po- tential,and restricted totheconfined hadronicphase, can bedirectlycompared tothat obtained from thepartitionfunction(14)ofthehadronresonancegas,usingtheform(15)introducedby Hagedorn for a continuum mass spectrum. Alternatively, as first approximation, one can also consider a discrete mass spectrum which accounts for all experimentally known hadrons and resonances. In this case the continuum partition function of the Hagedorn model is expressed by Eq. (13) with r (m ) replaced by the spin degeneracy factor of the ith hadron, and with the i summationtaken overallknownresonancespecies listedby theParticleDataGroup [33]. With the above assumption on the dynamics and the mass spectrum, the Hagedorn resonance gas partition function [1,5] can be calculated exactly and expressed as a sum of one–particle partitionfunctionsZ1 ofall hadronsand resonances, i lnZ(T,V)=(cid:229) Z1(T,V). (16) i i 10