---
address:
- |
    Niels Bohr Institute, Blegdamsvej 17, DK - 2100 Copenhagen, Denmark\
    and\
    Institute of Physics and Astronomy, University of Aarhus, DK - 8000 Aarhus\
    E-mail: neergard@nbi.dk
- |
    Niels Bohr Institute, Blegdamsvej 17, DK - 2100 Copenhagen, Denmark\
    E-mail: bondorf@nbi.dk
- |
    Frankfurt University, D-60054, Germany\
    and\
    The Kurchatov Institute, Russian Research Center, 123182 Moscow, Russia\
    E-mail: mishustin@nbi.dk
author:
- 'G. Neergaard'
- 'J. P. Bondorf'
- 'I. N. Mishustin'
title: Thermodynamics of explosions
---

Introduction
============

The assumption of thermodynamic equilibrium at an intermediate stage of a heavy-ion collision is often incorporated in models of the colliding nuclear matter. These models range from statistical models of nuclear multifragmentation to the fluid dynamical models of the quark gluon plasma. In contrast, microscopic models of molecular dynamics type (e.g. RQMD, FMD and NMD), which are based upon constituent interactions, do not contain this assumption. Such models are appropriate for testing to what extent thermodynamic equilibrium is actually achieved. And if it is not, the application of thermostatic concepts such as temperature and entropy becomes questionable. In this study we employ a very simple model, and focus on the thermodynamic or “overall” description of the system.

The model
=========

Our model consists of a number $A$ of identical balls of radius $r_{hc}$ having mass $m$. They perform classical non-relativistic hard-sphere scatterings, conserving energy, momentum and angular momentum. Initially the $A$ balls are placed randomly within a sphere of radius $R=R_0 A^{1/3}$, and the initial velocities are chosen as a superposition of thermal (Maxwell-Boltzmann) and collective motion. We use a spherically symmetric Hubble-like flow field for the initial collective motion: () = -v\_[0f]{}/R \[hubbleflow\] where $v_{0f}$ is a model parameter, $v_{0f}>0$ for ingoing flow and $v_{0f}<0$ for outgoing flow. We fix the total energy $E=E_{fl}+E_{th}$, and vary the fraction $\eta$ of the flow energy, $\eta=E_{fl}/E$, where $E_{fl}=\frac{m v_{0f}^2}{2 R^2}\sum_{i=1}^A\vec{r}_i^2$ and $E_{th}$ are the flow energy and the thermal energy, respectively. Because of the way in which the system is built up, these energies will fluctuate from event to event with a relative uncertainty of the order of $A^{-\frac{1}{2}}$.

In our simulations we have chosen nuclear-scale parameters: $m=940$ MeV, $r_{hc}=0.5$ fm, $R_0=1.2$ fm, $0 \leq v_{0f} \leq 0.5$ (in units of the velocity of light, $c=1$), but since the behavior of the model only depends on the two combined parameters $mv_{0f}^2$ and $r_{hc}/R$, the choice of nuclear scale is not crucial. We choose $A=50$, so with these parameters the initial radius of the system is $4.2$ fm.

We focus on four different types of event: $\bullet$ ’th20’: The particles are started in $100\%$ thermal motion inside a spherical container of radius $4.2$ fm, at $t=20$ fm/c the container walls are removed. $\bullet$ ’in’: $100\%$ ingoing flow. After interacting, the particles will move out again. This implosion-explosion process is intended to simulate some features of a heavy-ion collision. $\bullet$ ’50/50’: $50\%$ thermal motion + $50\%$ outgoing flow, simulating an explosion from a non-thermalized state. $\bullet$ ’100out7’: 100% outgoing flow inside a spherical container of radius $7$ fm.

The results are averaged over an ensemble of $20$ events of each kind.

Thermodynamic considerations
============================

It is clear that we cannot use ordinary thermodynamics (or its well-known extensions to small systems [@hill] or to small deviations from equilibrium [@zubarev]) for the description of the overall behavior of our model. First, it is not clear that equilibrium prevails, even locally. Indeed we wish to investigate to what extent equilibrium is reached in the course of an implosion-explosion process. Second, our system has no fixed volume, it expands freely into the vacuum. It is the combination of these two facts, no temperature and no volume, that makes our approach different from much previous work on the subject [@prevstud].

Equilibrium thermodynamics is linked to the motion of the individual constituents making up the macroscopic system via the entropy [@kennard; @ll]. A natural starting point for the investigation of the overall, [*i.e.*]{} the “thermodynamic”, behavior of our system is therefore to apply an expression similar to the entropy, but in a way that makes sense in this highly non-equilibrium system.

To study one-body observables, we reduce the $6A$ dimensional phase-space of the $A$ particles to $6$ dimensions in the standard way [@kennard]. Then we introduce a finite grid in the reduced phase-space, dividing each of the $6$ axes into $D$ segments. Instead of working with a fixed grid in phase-space, which would give us the usual entropy[^1], we let the entire grid expand or contract along with the swarm of points in phase-space in a uniform way: The outer grid edge follows the outermost point, the boxes are of equal size, and the number of boxes is kept fixed, thus the physical size ([*e.g.*]{} in units of $\hbar^3$) of each box in phase-space varies with time. This is to deal with the [*no volume*]{} problem, we mentioned above. We then introduce the [*pseudo-entropy*]{} as $$\Sigma = -\frac{1}{\xi} \sum_i p_i \log(p_i)
\label{pS}$$ where $$p_i=\frac{{number\ of\ points\ in\ box}\ i}{total\ number\ of\ points\ in\ phase\ space}
\label{plogp}$$ and $\xi$ is a normalization constant. We choose $\xi$ as the theoretical maximum value of $-\sum_i p_i \log(p_i)$ so that $\Sigma \in [0,1].$ In the current set-up, for the case of $N$ points in a reduced phase-space divided into $D^6$ boxes, $\xi$ is the smaller number of $\log(N)$ and $6\log(D)$. We refer to the quantity $\Sigma$ as the pseudo-entropy instead of entropy, since some important features of the entropy, [*e.g.*]{} that it increases, are not retained in this formulation. Nevertheless, we shall see that $\Sigma$ has some nice properties, including that of characterizing the degree of equilibrium. For a system of fixed volume in equilibrium, $\Sigma$ is the usual one-body entropy (apart from normalization).

Results
=======

In the calculations presented here, we have chosen $D=7$, so the $6$ dimensional phase-space is divided into $7^6$ boxes. Because we are dealing with a small number of events (typically $20$) in the ensembles, the precise values of $\Sigma$ depend on the choice of grid. We have, however, verified that $\Sigma$ behaves qualitatively similar to what is shown here over a range of grid-sizes, varying $D$ between $2$ and $9$.

To give an idea of the dynamics and the timescales, we show in Fig. \[fig:scatrate\] the scattering rate.

=[18pc]{}

Fig. \[fig:pseudos\] shows how the pseudo-entropy behaves in each of the four cases.

=[18pc]{}

In the ’th20’ case, the pseudo-entropy $\Sigma \simeq 1$ as long as the particles are in equilibrium at fixed volume inside the container. Then at $t=20$ fm/c, when the container is removed and the system starts to expand, the pseudo-entropy decreases, reflecting the fact that the system goes out of equilibrium[^2].

In the case ’100out7’, where the particles are started in an extreme non-equilibrium situation, the pseudo-entropy is low ($\Sigma=0.8$ is a low value in this context), but increases towards $\Sigma=1$ as the scatterings equilibrize the system. By comparison with Fig. \[fig:scatrate\] one can see that the first jump in $\Sigma$ at $t\simeq 5$ fm/c is due to particles scattering against the container wall (when particles hit the wall their velocity is reversed, so in this process many new states in phase-space are being populated), and the second “jump” around $t\simeq 20$ fm/c is due to the many particle-particle scatterings around this time.

The interesting case is the implosion-explosion (’in’) scenario, since here we do not know in advance if the system reaches a state of equilibrium or not. From the fact that the pseudo-entropy in Fig. \[fig:pseudos\] reaches a value $\Sigma \simeq 1$, the same value as the known equilibrium case ’th20’, we infer that the system [*is*]{} in a state of equilibrium around $t \simeq 6$ fm/c (which is also the time of maximum compression). We have checked that the speed distribution of the particles becomes nearly Maxwellian from $t \simeq 6$ fm/c with a temperature in the compressed state of $47$ MeV, which is also the theoretical value of the temperature assuming that all of the initial flow energy is converted to thermal energy.

Another interesting feature of the pseudo-entropy is that it seems to decay in a characteristic fashion when the system expands from a state of equilibrium, see Fig. \[fig:logpseudos\].

=[18pc]{}

Conclusions
===========

In this note we address the problem of thermodynamic equilibration in the context of heavy-ion collisions. We have defined a variable inspired by the entropy which, at least for the cases we have considered, seems to characterize the degree of equilibrium in an [*a priori*]{} highly non-equilibrium process such as an explosion. Now, more theoretical work needs to be done in order to understand [*why*]{} $\Sigma$ behaves in this seemingly interesting way.

Acknowledgements {#acknowledgements .unnumbered}
================

This work was in part supported by the Danish Natural Science Research Council. GN thanks the Niels Bohr Institute for kind hospitality and The Leon Rosenfeld Foundation for financial support during the preparation of this work. JPB thanks the Nuclear Theory Group at GSI, where part of this work was done. INM thanks The Humbold Foundation for financial support.

[99]{} See e.g. T.L. Hill: [*Thermodynamics of Small Systems I-II*]{}, W.A. Benjamin, NY, 1963-64. See e.g. N.D. Zubarev: [*Nonequilibrium Statistical Thermodynamics*]{}, Consultants Bureau, NY, 1974. To mention a few: J. Randrup, Nucl. Phys. [**A314**]{} (1979) 429–453; J. D. Bjorken, Phys. Rev. D [**27**]{} (1983) 140–151; G. Baym, Phys. Lett. [**138B**]{} (1984) 18–22; H. Heiselberg and X.-N. Wang, Phys. Rev. C [**53**]{} (1996) 1892–1902 and references herein. E.H. Kennard: [*Kinetic Theory of Gases*]{}, McGraw-Hill, 1938. Landau and Lifshitz: Statistical Physics I, Pergamon Press, 1968.

[^1]: in the limit $D\rightarrow \infty$ in an ensemble of infinitely many events

[^2]: The particles stay almost [*thermalized*]{}, though, in the sense that they retain their Maxwell-Boltzmann velocity distribution. But they are certainly not in [*equilibrium*]{}, since this means that the phase-space distribution is independent of time.
