## Brazilian Journal of Physics

*Print version* ISSN 0103-9733

### Braz. J. Phys. vol.37 no.2c São Paulo June 2007

#### http://dx.doi.org/10.1590/S0103-97332007000500006

**FLUCTUATIONS AND CORRELATIONS**

**Source chaoticity in relativistic heavy ion collisions at SPS and RHIC **

**Kenji Morita ^{I,} ^{*}; Shin Muroya^{II}; and Hiroki Nakamura^{I}**

^{I}Department of Physics, Waseda University, Tokyo 169-8555, Japan

^{II}Matsumoto University, Matsumoto 390-1295, Japan

**ABSTRACT**

We investigate degree of coherence of pion sources produced in relativistic heavy ion collisions using multi-particle interferometry. In order to obtain ''true'' chaoticity, l^{true} from two-pion correlation functions measured in experiments, we make a correction for long-lived resonance decay contributions. Using this l^{true} and the weight factor which are obtained from parameter fitted to two- and three-pion correlation function, we calculate a chaotic fraction e and the number of coherent sources a for different colliding energies. The result gives constraints on the source and shows an increase of the minimum number of a with multiplicity, although multiplicity independent chaoticity is not excluded.

**Keywords:** Relativistic heavy ion collisions; Pion interferometry

**I. INTRODUCTION**

Relativistic heavy ion collisions provide us an unique opportunity to explore nature of hot and dense nuclear matter on the earth. In the highest energy collisions at the BNL-RHIC, it is expected that the matter created soon after the collision of two heavy nuclei can be the strongly interacting quark-gluon plasma, which gradually cools down and then becomes hadronic matter via phase transitions. To understand the nature of the QCD matter, it is important to know what kind of information experimental observables contain. Pion interferometry has been one of the most important observables because it can give us information on sizes of the source which pions come from, through the HBT effect. The HBT effect is a quantum-mechanical effect due to symmetrization of two-boson wave function and occurs if the source is not completely coherent. The strength of the two-boson momentum intensity correlation takes its maximum value in the case of the perfectly chaotic source. Because the chaoticity can be related to the degree of thermal equilibrium and to composition of the source, analyses of the chaoticity can provide information on a state of hadronic matter which may reflect how hadrons are produced.

Experimentally measured two-pion chaoticity l = *C*_{2}(*p,p*)-1 does not usually reflect real coherence of the source because of long-lived resonance decay contributions. Three-pion correlation function is more useful for this purpose. The three-pion correlator measured in some experiments is

where *Q _{ij}* = ,

*Q*

_{3}= , and

*C*denotes the n-particle correlation function. The weight factor w =

_{n}*r*

_{3}(0)/2 characterizes the degree to how pions sources are chaotic. For a completely chaotic source, w = 1. If the source chaoticity is really a physical quantity, both two- and three-pion correlation function should give quantitatively consistent value of the chaoticity. However, experimental results do not seem so, because of the apparent reduction of the l due to long-lived resonances. But it has been shown that we can impose stronger constraints by using both of two- and three-pion correlation data [1].

In this work, we find that such decay contributions can be eliminated by making use of a statistical model [2] and applied it to various experimental data measured in SPS and RHIC experiments [3]. Then, we obtain the weight factor w from two- and three-pion correlation function data in these experiments. Using the resonance-corrected l which we call l^{true} and w , we calculate a chaotic fraction e and the number of coherent sources based on a partially multicoherent source model [4]. From this result, we discuss how the structure of the sources changes from low energy collisions at SPS and higher ones at RHIC [5].

**II. CORRECTION OF l FACTOR FOR LONG-LIVED RESONANCE DECAY**

In the presence of long-lived resonance decay contributions to two-pion correlation function, the two-pion chaoticity for a chaotic source is reduced as

where *N*_{ p } is the total number of pions and is the number of pions coming from long-lived resonance decay [6]. Hence, if one can estimate the ratio of pions from the decay to the total number, one can obtain apparent reduction factor caused by the decay contribution. Indeed, it is possible with the help of the statistical model which has been known to give a nice description of observed particle number ratio. Under the assumption of the constant chemical freeze-out temperature and chemical potential and boost-invariance along collision axis, the ratio of particle number can be written as the ratio of number density, *i.e.,N _{i }/ N_{j} = n_{i }/ n_{j}*, where

*n*is given by

_{i}with *f*(*E,T*, µ ) being the equilibrium distribution functions [7]. Using these formulae, we calculate the l^{eff} for S+Pb and Pb+Pb collisions at the SPS and Au+Au collisions at the RHIC and obtain l^{true} through a relation l^{true} = l^{exp}/ l^{eff} where l^{exp} is momentum-averaged experimental data [2, 3]. The temperature and baryonic chemical potential are determined by the c^{2} fitting to experimental data of particle ratio. See Ref. [3] for details. Here we summarize the result in Table I.

**III. EXTRACTION OF THE WEIGHT FACTOR**

In order to obtain the weight factor w , we have to extrapolate experimentally measured *r*_{3}(*Q*_{3})[Eq. (1)] to *Q*_{3} = 0. Using a simple source function in which instantaneous emission and spherically symmetric source are assumed, we construct the two- and the three-pion correlation functions with the formulae

where *f _{ij}* = 1/ is the source function and

*R*, l

_{inv}and n are parameters which should be determined by a simultaneous c

^{2}fit to the two- and the three-pion correlation functions[10]. From a set of these parameters, we can calculate w using Eqs. (4), (5), and (6). The results for w are shown in Table II.

**IV. ANALYSIS WITH THE PARTIALLY MULTICOHERENT MODEL**

From considerations in previous sections, we have obtained the two quantities, l^{true} and w as experimental results. Next, we investigate the coherence of the sources using the partially multicoherent model [4]. In this model, there are two characteristic parameters which are related to the l^{true} and w as

By solving the above equations with respect to e and a (this can be analytically), we can obtain allowed regions for e and a corresponding to the available range of l^{true} and w .

The result is shown in Fig. 1. In each of figures, the lightest shade area labelled ''A'' and the second one labelled ''B'' denote the allowed regions coming from l^{true} and w , respectively. The darkest areas labelled ''C'' are overlap of area ''A'' and ''B'', then correspond to the allowed parameter region for e and a . The best fit points are indicated by the filled box.

From Fig. 1, it seems to be difficult to find systematic change of the allowed regions. This result mainly comes from the fact that l^{true}'s are close to unity in the Pb+Pb data. However, it has been suggested that Coulomb correction to the two-pion correlation functions is over-corrected one; the values of l^{exp} can decrease if we take account of the partial Coulomb correction [8, 9]. Hence, if appropriate corrections were made, obtained l^{true} becomes smaller. In order to obtain a rough sketch of a tendency, we also draw the allowed regions for which is 0.8 times smaller than l^{true}. For S+Pb data, we multiply the original l^{true} by factor 0.7 because Gamow (point-like source) correction is made for this data.

From Fig. 2, we can see that all allowed regions (the darkest shaded areas) are narrow but there seem to exist systematics. The best fit point seems to move upper left (small chaotic fraction and large number of coherent sources) side, except for the NA44 Pb+Pb result. Most important feature of this result is that the upper limit of e and lower limit of a are determined by the lower limit of w . We plot the maximum and the minimum values of a in Fig. 3. The clear increase of minimum number of the coherent sources can be seen as a function of multiplicity while maximum number of those shows no such clear tendency.

In summary, we have given an analysis of the degree how chaotic the pion sources are in relativistic heavy ion collisions at the SPS and the RHIC. The analysis can be done by using both two-pion correlations and three-pion correlations. We find that the correction for long-lived resonance decay contributions to the two-pion correlation function can be subtracted with the help of the statistical model. From a point of view in which multicoherent sources and a background chaotic source are produced, we show that the model gives constraints on the structure of the source. Although the maximum number of the coherent sources does not show a clear multiplicity dependence, the minimum number of coherent source increases as the multiplicity increases.

**Acknowledgements**

The authors would like to thank Prof. I. Ohba and Prof. H. Nakazato for their encourgement. K. M's work is supported by a Grant for the 21st Century COE Program at Waseda University from Ministry of Education, Culture, Sports, Science and Technology of Japan and BK21 (Brain Korea 21) program of the Korean Ministry of Education.

[1] H. Nakamura and R. Seki, Phys. Rev. C **66**, 027901 (2002). [ Links ]

[2] K. Morita, S. Muroya, and H. Nakamura, Prog. Theor. Phys. **114**, 583 (2005). [ Links ]

[3] K. Morita, S. Muroya, and H. Nakamura, Prog. Theor. Phys. **116**, 329 (2006). See also references therein. [ Links ]

[4] H. Nakamura and R. Seki, Phys. Rev. C **61**, 054905 (2000). [ Links ]

[5] Because of space limitation, here we only show result for the partially multicoherent source model though we also analyzed the data with partial coherent model and multicoherent model in Refs.[2] and [3].

[6] T. Csörgo, B. Lörstad, and J. Zimányi, Z. Phys. C **71**, 491 (1996). [ Links ]

[7] J. Cleymans and K. Redlich, Phys. Rev. C **60**, 054908 (1999). [ Links ]

[8] D. Adamová et al. (CERES Collaboration), Nucl. Phys. A **714**, 124 (2003). [ Links ]

[9] S. S. Adler et al. (PHENIX Collaboration), Phys. Rev. Lett. **93**, 152302 (2004). [ Links ]

[10] Eq. (5) is a simple reduction of the correlator part of the fully chaotic case shown in Ref. [4]. See Ref. [2] for discussion.

Received on 31 October, 2006; Revised version received on 8 January, 2007; Third version on 26 February, 2007

* Present address: Institute of Physics and Applied Physics, Yonsei University, Seoul 120-749, Korea. Electronic address: morita@phya.yonsei.ac.kr