Information theory of open fragmenting systems

An information theory description of finite systems explicitly evolving in time is presented. We impose a MaxEnt variational principle on the Shannon entropy at a given time while the constraints are set at a former time. The resulting density matrix contains explicit time odd components in the form of collective flows. As a specific application we consider the dynamics of the expansion in connection with heavy ion experiments. Lattice gas and classical molecular dynamics simulations are shown.


INTRODUCTION
The microscopic foundations of thermodynamics are well established using the Gibbs hypothesis of statistical ensembles maximizing the Shannon entropy [1].At the thermodynamic limit, the various Gibbs ensembles converge to a unique thermodynamic equilibrium.On the other hand in finite systems the various Gibbs ensembles are not equivalent [2] and the physical meaning and relevance of these different equilibria has to be investigated.
A common interpretation of a statistical ensemble for a finite system is given by the Boltzmann ergodic assumption.In this interpretation the statistical ensemble represents the collection of successive snapshots of a physical system evolving in time, and the state variables are identified with the conserved observables.Nevertheless this image is blured when finite observation times are considered.Moreover the specification of the proper boundary conditions (see below) require infinite information.Finally, in heavy ion collisions, the systems are free to expand and thus nonstationary equilibrium is present.
In this communication we will address the problem of time dependent process in which the density matrix as well as boundary conditions depend on time leading to the appearance of new time-odd constraints or collective flows.

STATISTICAL EQUILIBRIA
When the system is characterized by L observables known in average < Â >= Tr D Â , statistical equilibrium corresponds to the maximization of the constrained entropy where D = ∑ (n) Ψ (n) p (n) Ψ (n) is the density matrix, and λ = {λ } are Lagrange multipliers.Gibbs equilibrium is then given by where Z λ is the associated partition sum.It should be noticed that microcanonical thermodynamics [3] can also be obtained from the variation of the Shannon entropy in the special case of a fixed energy subspace.In this case the maximum of the Shannon entropy can be identified with the Boltzmann entropy max (S) = logW (E), where W is the total number of states corresponding to an energy E. In the following we shall confine ourselves to the Gibbs formulation (1), which is more general than the microcanonical ansatz.Indeed the microcanonical density matrix corresponds to an even occupation of the whole energy shell while non ergodic components can be already included within the Gibbs formalism through the introduction of extra constraints.

Boundary condition problem in finite systems
The statistical physics formalism recalled above is valid for any system size .However, as soon as one Â contains differential operators such as a kinetic energy, eq. ( 1) is not defined, unless boundary conditions are specified.Only at the thermodynamic limit boundary conditions are irrelevant, as only in this limit surface effects are negligeable.The definition of any density with a finite number of particles requires the definition of a finite volume.To this aim, a fictitious container is generally introduced [4].The volume and shape of this unphysical box has no influence on the thermodynamics of self-bound systems, but in the presence of continuum states the situation is different.Let us consider the standard case of the annulation of the wavefunction on the surface S of a containing box V. Introducing the projector, PS , over S and its exterior, the boundary conditions reads PS Ψ (n) = 0 for all microstates (n).Using P2 S = PS , we can see that this condition imposes an extra constraint to the statistical ensemble < PS >= Tr D PS = 0.The density matrix then reads which shows that the thermodynamics of the system depends on the whole surface S. To specify the density matrix, the projector PS has to be exactly known at each point of the surface.This infinity of points corresponds to an infinite amount of information to be known to define the density matrix (2).This requirement is in contradiction with the statistical mechanics principle of minimum information.

Incomplete knowledge on the boundaries
One way to get around the difficulties encountered to take into account our incomplete knowledge on the boundaries, is to introduce a hierarchy of observables describing the size and shape of the matter distribution.
For example, if only the average system size < R2 > is known, the minimum information principle implies which is akin an isobar canonical ensemble, since the additional Lagrange multiplier λ R 2 imposing the size information has the dimension of a pressure when divided by a typical scale R 0 and by the temperature, λ R 2 = β PR 0 .
A typical application of this concept is the so called freeze-out hypothesis in nuclear collisions : at a given time t 0 the main evolution (i.e. the main entropy creation) is assumed to stop and partitions are supposed to be essentially frozen.Typically thermal and chemical equilibrium is assumed, meaning that the information at t 0 on the energetics and particle numbers is limited to the observables < Ĥ > and < Nf > for the different species f [5,4].Freeze-out occurs when the system has expanded to a finite size.Then at least one measure of the system's compactness should be included .The limited knowledge of the system extension leads to a minimum biased density matrix given by eq.( 3) [6].

MULTIPLE TIME STATISTICAL ENSEMBLES
As soon as one of the constraining observables Â is not a constant of the motion, the statistical ensemble (1) is not stationary.A single time description may still look appropriate in the freeze-out configuration discussed in the last section.Indeed in many physical cases one can clearly identify a specific time at which the information concentrated in a given observable is frozen.However this freeze-out time is in general fluctuating and different for different observables.
Let us now suppose that the different informations on the system, < Â >, are known at different times t : < Â > t = Tr D (t ) Â .A generalization of the Gibbs idea is that at a time t the least biased state of the system is the maximum of the Shannon entropy, considering all informations as constraints.
The maximization of the entropy at time t with the various constraints < Â > t known at former times t corresponds to the free maximization of where the λ are the Lagrange parameters associated with all the constraints.This maximization will lead to a density matrix which can be considered as a generalization to time dependent processes of the Gibbs ensembles (1).
Let us consider the case of a deterministic evolution where {., .}are Poisson bracket in classical physics and commutators divided by i h in quantum physics.The minimum biased density matrix, which corresponds to the solution of δ S c = 0, is given in matrix notation by [7] D where the partition sum is defined as and Â represent the time evolution of the constraining observables Â in the Heisenberg representation Â takes into account the time difference between the various observations.Eq. ( 6) can be interpreted as the introduction of additional constraints B(p) and additional Lagrange parameters ν (p) associated with the time evolution of the system Eq.( 7) is an exact solution of the complete many body evolution problem eq.( 5) with a minimum information hypothesis on the final time t having made few observations < Â > at previous times t , which shows the wide domain of applicability of information theory.A generalization of this theory to non deterministic evolutions can be found in ref. [7].We can see from eq.( 7) that in general an infinite amount of information, i.e. an infinite number of Lagrange multipliers are needed if we want to follow the system evolution for a long time.However, different interesting physical situations exist, for which the series can be analytically summed up.In this case, a limited information (the knowledge of a small number of average observables) will be sufficient to describe the whole density matrix at any time, under the unique hypothesis that the information was finite at a given time.

THE DYNAMICS OF THE EXPANSION
Let us now apply the above formalism to transient unconfined systems.We shall assume that at a given freeze out time t 0 the system can be modelized as a non interacting ensemble of n = 1, . .., N particles or fragments, and a definite value for the mean square ) characterizes the ensemble of states.Then we have to introduce the constraining observable Â2 = R2 associated with a Lagrange multiplier λ 0 .If time is not taken into account, the maximum entropy solution is given by Dβλ Eq.( 8) is akin to a system of non-interacting particles trapped in an harmonic oscillator potential with a string constant k = 2λ 0 /β .From the partition sum, the EOS are easily derived: is not an external confining potential but only a finite size constraint, the minimum biased distribution ( 8) is not stationary.To take into account the time evolution, we must introduce additional constraining observables B(1) R } = 0, all the other B(p) R with p > 2 are zero.The density matrix is given by with The density matrix (9) can be interpreted as a radially expanding ideal gas.Indeed the distribution can be written as where α = ν 0 (t)/β e f f (t) represents a Hubblian factor and the confining Lagrange multiplier is transformed into The term mα (t) ˆ r n correcting the momentum in eq.( 11) can be interpreted as a momentum produced by a radial velocity α (t) ˆ r n .This proportionality of the velocity with rn shows that the motion is self similar.As a consequence, when this collective motion is subtracted from the particle momentum, the density matrix (11) corresponds at any time to a standard equilibrium (8) in the local rest frame.
In this case the infinite information which is a priori needed to follow the time evolution of the density matrix according to eq.( 7), reduces to the three observables ˆ r 2 , ˆ p 2 , ˆ r • ˆ p + ˆ p • ˆ r.Indeed these operators form a closed Lie algebra, and the exact evolution of (11) preserves it algebraic structure.
Radial flow emerges as a necessary ingredient of any statistical description of an unconfined finite system, which is missed in the standard Gibbs ansatz in a confining box [4].One should then consider a more general equilibrium of a finite-size expanding system with β , α and λ 0 as free parameters.Then, if the observed minimum biased distribution at time t is coming from a confined system at time t 0 , the three parameters β , α and λ 0 should be linked to the time t 0 the initial temperature β −1 and the initial λ 0 by equations ( 10) and ( 12) [8].

NUMERICAL SIMULATIONS
As we have already mentioned in section , in the hypothesis of negligeable interaction between the system's constituents the expansion is self-similar, implying that the situation is equivalent to a standard Gibbs equilibrium in the local rest frame.In the expanding ensemble the total average kinetic energy per particle is simply the sum of the thermal energy e th = 3/(2β ) and the radial flow e f l = mα 2 r 2 /2.FIGURE 1. Upper left: effective pressure (full line) and average volume normalized to the ground state volume V 0 = A (dashed line) as a function of the collective radial flow.Upper right and lower part: fragment size distributions in the expanding Lattice Gas model.Distributions without flow (full black lines) are compared with distributions with 11% (dashed black), 43% (full grey) and 150% (dashed grey) contribution of radial flow at the same temperature (upper right) and at the same total energy (lower part).
This scenario is often invoked in the literature [4] to justify the treatment of flow as a collective radial velocity superimposed on thermal motion; however eq.( 11)contains also an additional term ∝ r 2 which corresponds to an outgoing pressure .The phase diagram and fragment observables are therefore expected to be modified by the presence of flow even in the self similar approximation.To quantify this statement, we have performed calculations in the Lattice Gas Model [8], and the results are shown in Figure 2.
The effective pressure λ e f f as well as the associated average volume are shown in the upper left part of figure 2 as a function of the collective radial flow for a given pressure and temperature .The Lagrange parameter λ e f f being a decreasing function of α, the critical point is moved towards higher pressures in the presence of flow [13].However one can see that the effect is very small up to mα ≈ .6 (which corresponds to Δp 2 /p 2 = e f l / e th ≈ 67% flow contribution).In this regime the cluster size distributions displayed in the upper right part of figure 2 are only slightly affected.On the other side if collective flow overcomes a threshold value Δp 2 /p 2 ≈ 100% the average volume shows an exponential increase and the outgoing flow pressure leads to a complete fragmentation of the system (dashed grey line in the lower part of fig.2).We can also observe that an oriented motion is systematically less effective than a random one to break up the system.This is shown in the lower part of figure 2 which compares for a given λ distributions with and without radial flow at the same total deposited energy: for any value of radial flow equilibrium in the standard microcanonical ensemble corresponds to more fragmented configurations.The above results are relevant for the experimental situation if and only if the interfragment interactions can be neglected when the system is still compact enough to bear pertinent information on the phase diagram.Indeed only in this case the series (7) can be analytically summed up and the expansion dynamics can be reduced to a self similar flow [7].The validity of the ideal gas approximation eq.( 9) for the expansion dynamics is tested in Figure 1 [10] in the framework of classical molecular dynamics [11].A Lennard Jones system is initially confined in a small volume and successively freely expanding in the vacuum.We can see that after a first phase of the order of ≈ 10 Lennard Jones time units, where interparticle interactions cannot be neglected, the time evolution predicted by eq.( 9) is remarkably fulfilled for all total energies.This result is due to the fact that the system's size and dynamics is dominated by the free particles, while deviations from a self similar flow can be seen if the analysis is restricted to bound particles [10].We expect eq.( 9) to describe the system evolution even better if the degrees of freedom n are changed from particles to clusters, as suggested by the Fisher model of condensation [12].

CONCLUSIONS
In this paper we have introduced an extension of Gibbs ensembles to time dependent constraints .This formalism gives a statistical description of a system observed at a time at which the entropy has not reached its saturating value yet, as it may be the case in intermediate energy heavy ion reactions [9].Another physical application concerns systems for which the relevant observables pertain to different times, as in high energy nuclear collisions [5].
Our most important result is that any statistical description of a finite unbound system must necessarily contain a local collective velocity term.Indeed the knowledge of the average spatial extension of the system at a given time, naturally produces a flow constraint at any successive time.Conversely a collective flow measurement at a given time can be translated into an information on the system density at a former time.

FIGURE 2 .
FIGURE 2. Time evolution of < R 2 >, < P 2 > and < R • P > for an initially constrained Lennard Jones system of 147 particles freely expanding in the vacuum, at different total energies.Lower right: expansion dynamics (symbols) compared to the prediction of eq.(9) (lines).