1. Introduction
NA61/SHINE is a fixed target experiment at the CERN Super Proton Synchrotron (SPS). One of its main aims is to study the phase diagram of QCD. In order to accomplish that, different collision systems at multiple energies are investigated. The NA61/SHINE detector is equipped with four large time projection chambers (TPC) [
1], these are covering the full forward hemisphere providing excellent tracking down to a transverse momentum of 0 GeV/
c. The experiment also features a modular calorimeter, located on the beam axis after the TPCs. This detector is called the projectile spectator detector, and it measures the forward energy which determines the collision centrality of the events. A setup of the NA61/SHINE detector system is shown in
Figure 1.
In order to study the QCD phase diagram and search for the critical end point (CEP), vastly different collision systems (p + p, p + Pb, Be + Be, Ar + Sc, Xe + La, Pb + Pb) are investigated at various beam momenta (13A, 20A, 30A, 40A, 75A and 150A GeV/c). There are many observables to accomplish this goal. In the analysis described in this paper we measure Bose–Einstein (or Hanbury–Brown and Twiss (HBT)) correlations of identical pions in Be + Be collisions at 150A GeV/c. These, based on the principles of quantum-statistical correlations, reveal the femtometer scale structure of pion production, hence this field is often called femtoscopy.
2. Femtoscopy, Lévy Sources and the Critical End Point
The method of femtoscopy is based on the work of R. Hanbury Brown and R. Q. Twiss [
2] as well as Goldhaber and collaborators [
3]. The key relationship of this method shows that the spatial momentum correlations (
) are related to the properties of the particle emitting source (
, describing the probability density of particle creation) in the following way:
where
is the Fourier transform of
, and
q is the momentum difference of the pair (dependence on the average momentum
K is suppressed here). See more details e.g., in ref. [
4]. The usual assumption for the shape of the source is, based on the central limit theorem, Gaussian. A generalization of this assumption is to assume Lévy distributed sources. A possible reason is that due to the expanding medium, the mean free path may increase and thus anomalous diffusion and Lévy distributed sources may appear [
5,
6]. Alternatively, due to critical fluctuations and the appearance of large scale spatial correlations, similar power-law tailed sources may be present [
7]. Another reason for Lévy distributed sources may be the fractal structure of QCD jets, as discussed in ref. [
8]. Here we restrict our investigation to symmetric Lévy distributions, as they have proven to provide a suitable description of Bose–Einstein correlations in nucleus-nucleus collisions [
4]. Furthermore, we restrict ourselves to describe the spatial part of the source, and the time dependence is absorbed through the connection of momentum difference
q and average momentum
K in case of identical particles:
Then the symmetric Lévy distribution is characterized by two parameters: Lévy scale parameter
R and the Lévy exponent
. The distribution is defined as follows:
This distribution can be expressed analytically in two special cases. One is the already mentioned Gaussian distribution for
; furthermore,
leads to a Cauchy distribution. An important difference between Lévy distributions and Gaussians is the presence of a power-law tail in case of
, i.e., for large distances (
r), the following holds:
where
d represents the number of spatial dimensions. With Lévy sources, the Bose–Einstein or HBT correlation functions can be expressed in the following way:
where the
intercept parameter was introduced, which is defined as
where the
extrapolation is done in experimentally available
q regions (limited from below by the two-track resolution in the given measurement). The Core-Halo model [
9,
10] may be utilized to understand this
intercept parameter. The core–halo model splits the source into two pieces. The core part contains the primordially created pions (directly from hadronic freeze-out) or from short lived (strongly decaying) resonances. The halo consists of pions created from longer lived (compared to the usual source size of a few femtometers) resonances and general background. In this picture, the
parameter turns out to be connected to the ratio of the Core and the Halo as follows:
Finally, let us come back to the above mentioned point of connecting Lévy sources to the search for the Critical End Point of QCD. Critical points are characterized by critical exponents, one of which is the exponent of spatial correlations. This appears because due to the second order phase transition at the CEP, the spatial correlation functions becomes a power-law with an exponent of
(where
d is the dimension and
is the critical exponent of spatial correlations). We can see that the Lévy exponent
, given in Equation (
4), defines a similar power-law, and hence
may be regarded as identical to the critical exponent
. See further discussions in Ref. [
7]. Critical exponents are universal in the sense that they take the same values in case of physical systems belonging to the same universality class. It has been shown [
11] that the universality class of QCD is that of the 3D Ising model. The value of the critical exponent
has been calculated to be
[
12]. Alternatively, one may rely on the universality class of the 3D Ising model with a random external field, in which case an
value of
was calculated [
13]. Considering the previous statements, if we “scan” the phase diagram with different energies and systems and measure the values of the
exponent, we might be able to gain more information on the location and characteristics of the CEP.
3. Measurement Details
In this measurement we analyzed the 0–20% most central Be + Be collisions at 150A GeV/c. This dataset consists of about three million events, which after various event and track quality selections was reduced to around 300,000 events. The track acceptance in this analysis was as follows. The rapidity region of analyzed particles is (corresponding to in the center-of-mass frame), the azimuthal coverage is ; n this track sample, we identified pions based on their deposited energy in the TPC gas and charge obtained from the curvature of their trajectories in the magnetic field. We then analyzed negative pion pairs and positive pion pairs, as well as the combination of these two (i.e., created a dataset of identically charged pion pairs). These pairs were sorted into four (average pair transverse momentum) bins in the range of 0–600 MeV/c. In each momentum bin, we measured the pair distribution of pairs from the same event, let us call this the actual pair distribution. This contains quantumstatistical correlations, as well as many other residual effects related to kinematics and acceptance. To remove this undesirable effects, we created a mixed event for each actual event, by randomly selecting particles from other events of similar parameters, and making sure particles are each selected from a different event. Let us call the pair distribution from this sample , the background distribution. Then the correlation function is calculated as , provided a proper normalization is done in a q range where quantumstatistical correlations are not expected. Let us mention here, that our analysis was done with a one-dimensional momentum difference variable q, calculated in the longitudinally co-moving system (LCMS), as in this frame, an approximately spherically symmetric source can be expected, furthermore, the extrapolation of is equivalent with the three dimensional case.
Final state effects are still present in the
correlation function. Among these, for pion pairs, the most important is the Coulomb effect, responsible for the repulsion of same charged pairs. It is usually handled by a so-called Coulomb correction as follows. The Coulomb correction for Lévy type sources is a complicated numerical integral to calculate and fit, as discussed in Ref. [
4]. However it can be observed that the Coulomb correction does not strongly depend on the Lévy exponent
. Hence we can use then the approximate formula published by the CMS Collaboration in ref. [
14], valid for
:
where
is the so-called Sommerfeld parameter, and
is the fine-structure constant (neither should be confused with the above discussed exponents
and
). Utilizing the usual Bowler–Sinyukov method (i.e., Coulomb correcting only for the Core part of the source) as indicated in Refs. [
15,
16], one obtains:
with
N being a normalization parameter responsible for the proper normalization of the
ratio. This is the final fit function we are using to describe our data.
As also found in ref. [
4], the Lévy parameters (
) are highly correlated, and especially in a low statistics dataset, it is hard to determine them precisely. We might be able to reduce this correlation and the statistical uncertainty of the parameters, if we fix one of the three parameters to a well motivated value. The resulting statistical uncertainties and free parameter values are modified due to the additional physical assumptions used to fix one of the parameters. From a statistical point of view, a bootstrap type of method may also be used. However, our main aim with this is to see a more clear trend of the
dependence of the parameters, with additional physical assumptions. One assumption is that
(i.e., the shape of the pion emitting source) is independent of
, with that we may fix
to a weighted average of the four
values obtained in free parameter fits performed in each
bin. The other option is fixing
R with the following equation motivated by hydrodynamical predictions of the particle emission homogeneity length (essentially the HBT radii) in case of expanding fireballs [
10]. In this case, we fit the following equation to the
dependence (where
, i.e., the average transverse mass of the pair) of the
R Lévy scale:
Previous results with free parameter fits were shown in Ref. [
17], hence here we concentrate on the results of the above mentioned fixed parameter fits. We again note that our aim with fixing one of the parameters to a physically motivated value is to show the trend of the
dependence of the parameters.
4. Results
First, the measured correlation functions were fitted with the above mentioned (Equation (
10)) function with three free parameters (
), as shown in Ref. [
17]. Using the results from the free parameter fit, we fitted a constant function to the
values for all
bins, as well as the formula of Equation (
11) to the
R values in each bin. Then we analyzed first with one parameter fixed. All results are shown in
Figure 2,
Figure 3 and
Figure 4. Let us note here, that all the measurement settings (event selection, track selection, pair cuts, fitting interval) were varied systematically to obtain an estimate of systematic uncertainties. These, along with the statistical uncertainty of the parameters (obtained by the Minos algorithm) are shown in
Figure 2,
Figure 3 and
Figure 4.
The Lévy stability exponent
determines the source shape, and a value of 2 corresponds to a Gaussian source, a value of 1 to a Cauchy source, and 0.5 is the conjectured value at the critical point. Our results, along with these special cases (dotted yellow lines), are shown in
Figure 2. It is clear from the figure that the statistical uncertainties of
from the fixed
R fit are reduced by a factor 4–5; however, the values of
are similar in both cases. These values are far from the Gaussian case as well as the conjectured CEP value, motivating us to perform this measurement in different systems and at different energies as well.
The Lévy scale
R determines the correlation length of pion pairs from the given system. From a simple hydrodynamical picture one obtains a
type of affine linear dependence, as already mentioned above. This, i.e., Equation (
11), describes the free parameter fit data points of
Figure 3 well. The fixed
fits are within uncertainties compatible with the free parameter results, however, they suggest a more or less constant trend of
R versus
in the higher
bins. This motivated us to perform the same measurement in collision systems with larger multiplicities (where statistical uncertainties are expected to be reduced proportionally to the square of the mean multiplicity).
The last parameter to study was the correlation strength parameter
, given in Equation (
7). The transverse mass dependence of
is shown in
Figure 4. Comparing the three different fits (free parameter fit, fixed
fit, fixed
R fit), it was visible that
in free parameter fit was compatible (within statistical uncertainties) with both other cases. All fits showed a roughly constant
trend. This was in contrast to the findings at RHIC, see e.g., the compilation in ref. [
18]. This finding was, however, compatible with previous SPS measurements (in different systems), see e.g., ref. [
19].