Abstract

We review topics related to the first moment of azimuthal anisotropy (), commonly known as directed flow, focusing on both charged particles and identified particles from heavy-ion collisions. Beam energies from the highest available, at the CERN LHC, down to projectile kinetic energies per nucleon of a few GeV per nucleon, as studied in experiments at the Brookhaven AGS, fall within our scope. We focus on experimental measurements and on theoretical work where direct comparisons with experiment have been emphasized. The physics addressed or potentially addressed by this review topic includes the study of Quark Gluon Plasma and, more generally, investigation of the Quantum Chromodynamics phase diagram and the equation of state describing the accessible phases.

1. Introduction

The purpose of relativistic nuclear collision experiments is the creation and study of nuclear matter at high energy densities. Experiments have established a new form of strongly interacting matter, called Quark Gluon Plasma (QGP) [17]. Collective motion of the particles emitted from such collisions is of special interest because it is sensitive to the equation of state in the early stages of the reaction [810]. Directed flow was the first type of collective motion identified among the fragments from nuclear collisions [1113] and, in current analyses, is characterized by the first harmonic coefficient in the Fourier expansion of the azimuthal distribution of the emitted particles with respect to each event’s reaction plane azimuth () [1416]:where is the azimuth of a charged particle, or more often, the azimuth of a particular particle species, and the angle brackets denote averaging over all such particles in all events. In some experimental analyses, is evaluated directly from (1); then a correction is applied for reaction plane resolution [15], whereas, in a typical modern analysis method, the directed flow correlation is extracted using cumulants [16]. In general, is of interest when plotted as a function of rapidity, , or sometimes pseudorapidity , where is the polar angle of the particle. The dependence of on collision centrality and on transverse momentum, , can offer additional insights.

Until relatively recently [17, 18], the rapidity-even component was always assumed to be zero or negligible in mass-symmetric collisions. In fact, fluctuations within the initial-state colliding nuclei, unrelated to the reaction plane, can generate a significant signal [17, 18]. This fluctuation effect falls beyond the scope of the present review, which focuses on fluid-like directed flow, , as per (1), and from here on, for mass-symmetric collisions implicitly signifies .

During the first decade of the study of in nuclear collisions, it was more commonly called sideward flow. It refers to a sideward collective motion of the emitted particles and is a repulsive collective deflection in the reaction plane. By convention, the positive direction of is taken to be the direction of “bounce-off” of projectile spectators in a fixed target experiment [8, 10]. Models imply that directed flow, especially the component closest to beam rapidity, is initiated during the passage time of the two colliding nuclei; the typical time-scale for this is [9, 10], where and are the nuclear radius and Lorentz factor, respectively. This is even earlier than the still-early time when elliptic flow, , is mostly imparted. Thus can probe the very early stages of the collision [19, 20], when the deconfined state of quarks and gluons is expected to dominate the collision dynamics [9, 10]. Both hydrodynamic [21, 22] and transport model [23, 24] calculations indicate that the directed flow of charged particles, especially baryons at midrapidity, is sensitive to the equation of state and can be used to explore the QCD phase diagram.

The theoretical work leading to the prediction of collective flow in nuclear collisions evolved gradually. In the mid-1950s, Belenkij and Landau [25] were the first to consider a hydrodynamic description of nuclear collisions. During the 1970s, as the Bevatron at Lawrence Berkeley National Lab was converted for use as the first accelerator of relativistic nuclear beams, the idea of hydrodynamic shock compression of nuclear matter emerged [2628], and these developments in turn led to increasingly realistic predictions [11, 12, 29] that paved the way for the first unambiguous measurement of directed flow at the Bevalac in the mid-1980s [13]. A frequent focus of theory papers during the subsequent years was the effort to use directed flow measurements to infer the incompressibility of the nuclear equation of state in the hadron gas phase and to infer properties of the relevant momentum-dependent potential [8, 10]. The observed directed flow at AGS energies [3038] and below is close to a linear function of rapidity throughout the populated region, and the slope can adequately quantify the strength of the signal. At SPS energies and above [3947], a more complex structure is observed in , with the slope in the midrapidity region being different from the slope in the regions closer to beam rapidity.

Various models also exhibit this kind of behavior. At these energies, both hydrodynamic and nuclear transport calculations predict a negative sign for charged particle near midrapidity, where pions are the dominant particle species. This negative near midrapidity has been given various names in the literature: “third flow component” [48], “antiflow” [49], or “wiggle” [10, 22, 50]. This phenomenon has been discussed as a possible QGP signature, and a negative for baryons has been argued [22] to be particularly significant. However, some aspects of antiflow can be explained in a model with only hadronic physics [50, 51] by assuming either incomplete baryon stopping with a positive space-momentum correlation [50], or full stopping with a tilted source [52].

A three-fluid hydrodynamic model [22] predicts a monotonic trend in net-baryon directed flow versus beam energy in the case of a purely hadronic equation of state, whereas a prominent minimum at AGS energies, dubbed the “softest point collapse,” is predicted when the equation of state incorporates a first-order phase transition between hadronic and quark-gluonic matter. Recent measurements of both proton and net-proton directed flow at RHIC [47] indeed indicate nonmonotonic directed flow as a function of beam energy, with the minimum lying between 11.5 and 19.6 GeV in . However, more recent hydrodynamic and nuclear transport calculations which incorporate significant theoretical improvements (see Section 5) do not reproduce the notable qualitative features of the data and therefore cast doubt on any overall conclusion about the inferred properties of the QCD phase diagram. Directed flow has also been measured at the LHC [53]. A negative slope of is observed for charged particles, but its magnitude is much smaller than at RHIC, which is thought to be a consequence of the smaller tilt of the participant zone at the LHC.

In this article, we review a representative set of directed flow results spanning AGS to LHC energies. In Sections 2 and 3, we discuss measurements of for charged particles in mass-symmetric and mass-asymmetric collisions, respectively. In Section 4, we cover measurements of for various identified particle species. Section 5 reviews some recent model calculations which lend themselves to direct comparisons with directed flow data. Section 6 presents a summary and future outlook.

2. Differential Measurements of Charged Particle Directed Flow

In this section, we review measurements of for all charged particles in cases where individual species were not identified. Studies of the dependence on transverse momentum , pseudorapidity , beam energy , system size, and centrality are included.

2.1. Dependence of on Transverse Momentum

The -dependence of for charged particles has been studied by the STAR experiment at RHIC [43, 44]. Figure 1(a) presents directed flow results for Au+Au collisions in two centrality intervals, 5–40% and 40–80%, and in two regions of pseudorapidity, and . In this case, because of the odd-functional property of , the backward pseudorapidity region by convention has its sign of reversed before summing over the indicated gate in pseudorapidity. These measurements represent the first instance of using only spectators to determine the estimated azimuth of the reaction plane [5456].

The measured in the midpseudorapidity interval crosses zero in the region of 1 to 2 GeV/. Extrapolations raise the possibility that a qualitatively similar zero crossing by forward pseudorapidity occurs at higher , but that region does not fall within the acceptance of the STAR detector.

The charged particle measurements by ALICE [53] in Pb+Pb collisions at 2.76 TeV are compared in Figure 1(b) with the corresponding data from 200 GeV Au+Au collisions at RHIC. The measurements at both LHC and RHIC show a similar trend, including a sign change around  GeV/ in central collisions and negative values at all for peripheral collisions.

There is an interest in the observation of zero crossing behavior in as it can be used to constrain hydrodynamic model calculations [57]. It has been pointed out that this sign change is an artifact of combining all species of charged particles together and can be explained [44] by the different sign of for pions and baryons, in conjunction with the enhanced production of baryons at higher [58]. This complication is one of the reasons why directed flow of identified particles, as reviewed in later sections, can be easier to interpret and offers additional insights.

2.2. Dependence of on Pseudorapidity

Figure 2(a) shows charged particle in Au+Au collisions at 19.6, 62.4, 130, and 200 GeV measured in the PHOBOS detector [42]. It is evident that charged particle within 3 to 4 units of on either side of has a sign opposite to that of the spectators on that side of (the antiflow phenomenon).

Figure 2(b) shows the dependence of for charged particles in 2.76 TeV Pb+Pb collisions, as measured by the ALICE collaboration [53] at the LHC. The ALICE results are compared in this panel with RHIC measurements from STAR [41, 44]. The slope at the LHC and at the top RHIC energies has the same negative sign, but the slope magnitude at the LHC is a factor ~3 smaller than at top RHIC energy. This pattern is consistent with the participant zone at the LHC having a smaller tilt, as predicted [52], and does not support a proposed picture at LHC energies in which a strong rotation is imparted to the central fireball [59, 60].

2.3. System Size and Beam Energy Dependence of

The beam energy and system size dependence of have been studied at RHIC using data from two colliding species at two beam energies. Figure 3(a) shows charged particle for midcentrality (30–60%) Au+Au and Cu+Cu collisions at and 200 GeV, measured by STAR [44]. A trend of decreasing is observed as beam energy increases for both Au+Au and Cu+Cu collisions. Across the reported pseudorapidity range, is independent of system size within errors at each beam energy. This is a remarkable finding, given that the Au+Au system mass is three times that of the Cu+Cu system, and given that neither the AMPT [6163] nor the UrQMD [23, 24] models exhibit such a scaling behavior.

A different scaling behavior is presented in Figure 3(b). Here, the data in Figure 3(a) are transformed into the rest frame of the beam nucleus; that is, zero on the -axis corresponds to for each of the two collision energies involved. Within errors, the measurements lie on a universal curve across about three units of pseudorapidity. This behavior is known in the heavy-ion literature as limiting fragmentation, and had previously been observed in Au+Au collisions as a function of beam energy by STAR [43] and PHOBOS [42]. The term “limiting fragmentation” was originally employed by Feynman [64] and Benecke et al. [65] to describe the analogous phenomenon of the measured ratio in cosmic ray showers at sea level being almost independent of muon energy.

2.4. Centrality Dependence of

Figure 4 shows as a function of collision centrality in 62.4 and 200 GeV Au+Au [44] from STAR and in 2.76 TeV Pb+Pb [53] from ALICE, for the midpseudorapidity region . Figure 4 also shows (divided by 6 in order to fit conveniently on a common scale) as a function of collision centrality in 200 GeV Au+Au [44] from STAR at the forward pseudorapidity region .

Directed flow magnitude at midpseudorapidity increases monotonically going from central to peripheral collisions at all three beam energies, and there is a strong trend for this magnitude to decrease with increasing beam energy. A similar but stronger centrality dependence is observed at forward pseudorapidity, and the magnitude also increases strongly from midpseudorapidity to forward pseudorapidity.

It has been pointed out by Caines [67] that many aspects of soft physics (defined as GeV/) in heavy-ion collisions at relativistic energies depend only on event multiplicity per unit rapidity; in other words, for a fixed value of , there is no significant dependence on beam energy, or on centrality, or on the mass of the colliding system. This type of scaling is called “entropy-driven” soft physics. The directed flow results for two beam energies and two colliding systems reported by STAR in [44] (see Figure 3) represented one of the first (and still few) violations of entropy-driven multiplicity scaling. In contrast, this scaling is observed to hold, with caveats, for homogeneity lengths from femtoscopy [67, 68], for elliptic flow per average participant eccentricity [67, 69], and for various strangeness yields [67].

3. Directed Flow of Charged Particles in Mass-Asymmetric Collisions

In mass-asymmetric collisions like Cu+Au, the well-defined distinction that exists in mass-symmetric collisions between the odd component (a hydrodynamic effect correlated with the reaction plane) and the even component (an initial-state fluctuation effect unrelated to the reaction plane) no longer holds. In a recent paper from the PHENIX collaboration, they report midrapidity charged hadron in Cu+Au collisions at GeV for centralities of 10–20%, 20–30%, 30–40%, and 40–50%, using spectator neutrons from the Au side of the collision to determine the event plane; see Figure 5 [66]. However, they preserve the standard convention for the sign of by defining the direction of bounce-off by remnants of the first nucleus in the A+A system (Cu) to be positive. An even more recent paper from STAR reports distributions for the same system and centrality that are consistent within errors [70].

The PHENIX results in Figure 5 [66] reveal that the higher particles at midrapidity (above 1 or 1.5 GeV/ in ) and at all the studied centralities have negative and so are preferentially emitted with azimuths parallel to the Au fragment bounce-off direction (and antiparallel to the Cu fragment bounce-off direction). Whether or not the more abundant particles below 1 GeV/ are preferentially emitted with opposite azimuths, as might be expected based on momentum conservation, cannot be answered within the systematic uncertainty of the measurements [66].

In the STAR collaboration’s analysis of charged particle directed flow in Cu+Au collisions at GeV, a particular focus is the difference between positive and negative charges. This difference has the potential to be sensitive to the strong electric field between the two incident ions, whose electric charges differ by units; this field has a lifetime on the order of a fraction of a fm/. Figure 6 [70] presents at medium (GeV/) and intermediate centrality (10–40%). Like the PHENIX result, this measurement was made relative to the event plane from spectator neutrons, dominated by the Au side. It is evident from Figure 6 that both even and odd components are present, and most interestingly, there is a significant pattern showing a larger magnitude for negative particles.

The Parton-Hadron String Dynamics (PHSD) model [72, 73], when the initial electric field is explicitly modeled, predicts a difference signal that is an order of magnitude larger than observed [70]. On the other hand, parton distribution functions [74] can be used to estimate the number of quarks and antiquarks at very early times in relation to the number created in the collision; then given certain plausible assumptions, as set out in [70], it can be inferred that only a small fraction of the total quarks created in the collision are produced during the lifetime of the initial electric field. In addition to this important insight, the charged-dependent directed flow measurements in Cu+Au collisions offer new and valuable quantitative information with relevance to the Chiral Magnetic Effect [75, 76] and the Chiral Magnetic Wave [77, 78].

4. Differential Measurements of Identified Particle Directed Flow

The charged particle measurements reviewed in Section 2 are an admixture of all emitted particle species. Measurements of directed flow for identified particles offer more insights into the underlying physics that controls this observable. In this section, we discuss the dependence of on , , and centrality for several identified particle species.

4.1. Dependence of on Transverse Momentum

Measurements of for protons, antiprotons, and charged pions have been reported by the E877 collaboration at the AGS (11 GeV/) [3033]. For antiprotons, large negative values of are observed for GeV/ but with large statistical errors. For protons and charged pions, in various rapidity gates have been published, and these results are also divided into various bins of transverse energy, , which is a proxy for centrality. The E877 collaboration also provides the information needed to convert from their intervals of into percent centrality [31].

The NA49 collaboration [39] measured proton and pion in Pb+Pb collisions at a projectile kinetic energy of 158 GeV, as shown in Figure 7. The rapidity gate for these measurements is , which corresponds to a forward region (midrapidity is ). The NA49 collaboration describes the behavior as “peculiar,” especially for pions. However, they point out that negative at low has been predicted by Voloshin [79] and is explained by the interaction of radial and directed flow. Various types of nonflow effect were also mentioned as possible contributors to the observed pion behavior [39].

4.2. Dependence of on Rapidity

Various models suggest that the structure of near midrapidity, especially the pattern for baryons, is sensitive to the QCD equation of state and therefore this signal can be used to investigate QGP production and changes of phase [810, 22, 80, 81]. Figure 8 presents proton and pion in central, midcentral, and peripheral Pb+Pb collisions at projectile kinetic energies of 40 and 158 GeV [40], as reported by the NA49 collaboration at the CERN SPS. The data points at negative rapidity are mirrored from the positive side. The for pions is similar in magnitude and shape at both 40 and 158 GeV. The proton measurements suggest that proton antiflow (also called wiggle) [10, 22, 4850], which is not observed at the AGS (see Figure 9), might begin to happen at SPS energies. However, based on detailed studies using a variety of methods, including the approach introduced by Borghini et al. [82], the NA49 collaboration reports that the observed pattern of proton at SPS energies could be influenced by nonflow effects [40].

The EOS-E895 experiment carried out a beam energy scan at the Brookhaven AGS during 1996. E895 featured the first Time Projection Chamber with pad readout and reported directed flow in the form of , as well as in the form of the older observable [83] (in-plane transverse momentum) for several identified particle species: , , , , and , in Au+Au collisions at projectile kinetic energies 2, 4, 6, and 8 GeV [3438]. Figure 9(a) shows for protons at the four E895 beam energies, where denotes normalized rapidity such that the target and projectile are always at and , respectively. The slope of remains positive for all four beam energies.

The slope of the rapidity dependence of directed flow was extracted by E895 using a cubic fit . Figure 9(b) shows the fitted proton slope (not using normalized rapidity) and, unlike in [34], the horizontal axis in Figure 9(b) uses the now-conventional beam energy scale . Two additional points for proton in Au+Au collisions are also plotted here: a measurement at the top energy of the Berkeley Bevalac (1.2 GeV) [71] using the same EOS TPC detector as in E895 and an E877 measurement at the top AGS energy of 11 GeV [31]. The plotted data show a peak near  GeV beam kinetic energy (GeV) and, thereafter, a smooth decrease with beam energy across the AGS range.

Three-fluid hydrodynamic calculations [22] predict a “softest point” in the Equation of State in the AGS energy range, but the E895 beam energy scan did not reveal any nonmonotonic behavior in the energy dependence. Overall, hadronic models with a momentum-dependent mean field [84] show better agreement with data at AGS and SPS energies.

Phase-I of the beam energy scan (BES) [85, 86] at RHIC took data in 2010, 2011, and 2014, spanning the range from 200 GeV down to 7.7 GeV. Lattice QCD calculations [8789] imply a smooth crossover from hadronic matter to QGP at beam energies near and above the top energy of RHIC, whereas phase diagram features like a first-order phase transition and a critical point may become evident in nuclear collisions as the beam energy is scanned across RHIC’s BES region. The STAR experiment initially reported measurements of directed flow for protons, antiprotons, and charged pions in the energy range of 7.7 to 200 GeV [47]. At the 2015 Quark Matter meeting, preliminary directed flow data at BES energies for nine particle species were presented, along with slopes for protons, , and in 9 bins of centrality [90]. Figure 10 shows at intermediate centrality (10–40%) for , , , , , , and at 7.7 to 39 GeV.

The overall strength of directed flow has been characterized by the slope from a linear fit over the range [90], and the beam energy dependence of these slopes for the same nine particle species (, , , , , , and ) in 10–40% centrality Au+Au collisions is presented in Figure 11.

The most noteworthy feature of the data is a minimum within the range of 10 to 20 GeV in the slope for protons at intermediate centrality. The same quantity for hyperons is consistent with the proton result, but the larger statistical errors for do not allow an independent determination of a possible minimum in the beam energy dependence for this species. The proton and directed flow slope change from positive to negative close to 11.5 GeV and remain negative at all the remaining energies (up to and including 200 GeV in the case of protons). The remaining species have negative slope at all the studied beam energies. At 7.7 GeV, for is closer to zero than for , which supports the inference that and experience nuclear potentials that are repulsive and attractive, respectively, under the conditions created at this beam energy [91].

4.3. Net-Particle Directed Flow

There are two separate contributions to the energy dependence of proton directed flow in the vicinity of midrapidity: one part arises from baryon number transported from the initial-state at beam rapidity towards by the stopping process of the collision, while the other part arises from baryon-antibaryon pairs produced in the fireball near midrapidity. Clearly, these two contributions have very different dependence on beam energy, and disentangling them has a good potential to generate new insights. Towards this end, the STAR collaboration has defined [47] net-proton directed flow according to where is the observed rapidity dependence of the antiproton to proton ratio. Net-kaon is defined analogously, with and substituted for and , respectively [90].

STAR’s measurements of net-proton and net-kaon directed flow slope as a function of beam energy are reproduced in Figure 12 [90]. The net-proton directed flow shows a double sign change, and a clear minimum around the same beam energy where the proton directed flow has its minimum. This coincidence is not surprising, since antibaryon production is low at and below the beam energy of the minimum, and therefore net-proton and proton observables only begin to deviate from each other at energies above the minimum.

The observed minimum in directed flow for protons and net protons resembles the predicted “softest point collapse” of directed flow [19, 22], and these authors are open to an interpretation in terms of a first-order phase transition. Other theorists (see Section 5) point out that mechanisms other than a first-order phase transition can cause a drop in pressure (a softening of the equation of state) and therefore a definitive conclusion requires more research.

The net-kaon directed flow reproduced in Figure 12 [90] shows close agreement with the net-proton result near and above 11.5 GeV but deviates very strongly at 7.7 GeV. This deviation is not understood. To date, all of the several model comparisons with STAR measurements at BES energies have considered only particle as opposed to net-particle .

5. Recent Model Calculations of Directed Flow

Models that explicitly incorporate properties of the Quantum Chromodynamics phase diagram and its equation of state (typically hydrodynamic models) suggest that the magnitude of directed flow is an excellent indicator of the relative pressure during the early, high-density stage of the collision. Therefore, directed flow at sub-AGS beam energies can reveal information about hadron gas incompressibility, while at the higher energies that are a focus of the present review, it can flag the softening, or drop in pressure, that may accompany a transition to a different phase, notably Quark Gluon Plasma. For example, there may be a spinodal decomposition associated with a first-order phase transition [92, 93], which would cause a large softening effect. However, interpretation of flow measurements is not straightforward, and it is known that directed flow can also be sensitive to poorly understood model inputs like momentum-dependent potentials [91] in the nuclear medium. More theoretical work is needed to elucidate the quantitative connection between softening signatures, like the beam energy dependence directed flow, and QCD phase changes.

During the period since publication of the STAR BES directed flow results in 2014, there have been several theoretical papers [19, 20, 91, 94, 95] aimed towards interpretation of these measurements. The Frankfurt hybrid model [96] used for the data comparison by Steinheimer et al. [95] is based on a Boltzmann transport approach similar to UrQMD for the initial and late stages of the collision process, while a hydrodynamic evolution is employed for the intermediate hot and dense stage. The equation of state for the hydro stage includes crossover and first-order phase transition options. The data comparison by Konchakovski et al. [20] uses the Parton-Hadron String Dynamics (PHSD) model [72] of the Giessen group, a microscopic approach with a crossover equation of state having properties similar to the crossover of lattice QCD [8789]. The PHSD code also has a mode named Hadron String Dynamics (HSD), which features purely hadronic physics throughout the collision evolution and which yields directed flow predictions in close agreement with those [47] of the UrQMD model. The data comparison by Ivanov and Soldatov [94] uses a relativistic 3-fluid hydrodynamic model (3FD) [97] with equations of state that include a crossover option and a first-order phase transition option. The most recent comparison to the STAR BES data, by Nara et al. [19], uses the Jet AA Microscopic (JAM) model [98]. JAM is a purely hadronic Boltzmann transport code, but the authors of [19] introduce an option to switch from the normal stochastic binary scattering style to a modified style where the elementary 2-body scatterings are always oriented like attractive orbits [99, 100]. They argue that the switch-over from random to attractive binary orbits mimics the softening effect of a first-order phase transition.

Figure 13 focuses on the most promising directed flow measurement from the RHIC Beam Energy Scan, namely, the for protons at 10–40% centrality, and summarizes recent model comparisons [19, 20, 94, 95] with these data. These authors are largely in agreement that the data disfavor models with purely hadronic physics. However, some conclude that a crossover deconfinement transition is favored [20, 94], while others conclude that a first-order phase transition is still a possible explanation [19]. Note that the argument of Nara et al. [19] is that a more sophisticated implementation of a first-order phase transition would transition from the “JAM” curve at low BES energies to the “JAM-attractive” curve at higher BES energies.

Overall, Figure 13 underlines the fact that no option in any of the model calculations to date reproduces, even qualitatively, the most striking feature of the data, namely, the minimum in proton directed flow in the region of GeV. It is also noteworthy that the difference between nominally equivalent equations of state implementations in different models is very large. For example, the differences in between the 1st-order phase transition in the hybrid model [95] and a similar nominal quantity for the 3FD model [94] (see Figure 14) are currently more than an order of magnitude larger than the experimental measurement being interpreted and are larger still than the error on the measured data.

6. Summary and Outlook

In this review, we discuss heavy-ion directed flow results for charged particles and for identified particle types, covering beam energies from the Brookhaven AGS to the CERN LHC. Charged particle directed flow measurements have been published as a function of transverse momentum, pseudorapidity, and collision centrality, while Cu+Cu and Au+Au have also been compared. The charged particle directed flow magnitude at the LHC is a factor of three smaller than that at top RHIC energy. The observations from RHIC suggest that the charged particle directed flow is independent of system size but depend on the incident beam energy. Limiting fragmentation scaling is observed for at RHIC energies, but entropy-driven multiplicity scaling in terms of is not seen at RHIC. In mass-asymmetric collisions, specifically Cu+Au, recent directed flow measurements at GeV have opened a new window into quark and antiquark formation at the very earliest times of the collision evolution ( fm/) and could clarify theoretical and experimental questions related to the Chiral Magnetic Effect and the Chiral Magnetic Wave.

Measurements of for identified species offer deeper insights into the development of hydrodynamic flow. Opposite for pions and protons at AGS/SPS and at lower RHIC energies suggests an important role for nuclear shadowing. Signals of antiflow of neutral kaons in AGS/E895 together with kaon measurements in the RHIC Beam Energy Scan region point to kaon-nucleon potential effects. The single sign change in proton slope and a double sign change in net-proton slope, with a clear minimum around ~11.5–19.6 GeV, show a qualitative resemblance to a hydrodynamic model prediction called “softest point collapse of flow.” This original prediction assumed a first-order phase transition, but a crossover from hadron gas to a deconfined phase can also cause softening (a drop in pressure). None of the current state-of-the-art models can explain the main features of the STAR directed flow measurements, and different models with nominally similar equations of state diverge from each other very widely over the BES range.

Looking ahead to likely developments during the period 2017-2018 in the area of directed flow at of a few GeV and above, we can expect new BES Phase-I results for the meson, as well as final publication of current preliminary RHIC Beam Energy Scan Phase-I results, like those in [90]. We can also expect parallel theoretical work on related physics and interpretation of the newest data. The preliminary results for new particle species like s and charged and neutral kaons, as well as BES for protons, s, and pions in narrow bins of centrality (nine bins’ spanning 0–80% centrality) amount to very stringent constraints on the next round of theoretical interpretation in terms of the QCD phase diagram and equation of state.

More comprehensive measurements of the many phenomenological aspects of directed flow outlined in this review will be possible beginning in the year 2019, as a consequence of the much increased statistics of Phase-II of the RHIC Beam Energy Scan (to take data in 2019 and 2020), in conjunction with the anticipated upgrades to the performance of the STAR detector [101]. Thereafter, new facilities like FAIR [102] in Germany, NICA [103] in Russia, and J-PARC-HI [104, 105] in Japan, will begin coming online. These new dedicated facilities will further strengthen worldwide research on the QCD phase diagram at high baryon chemical potential.

Competing Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work was supported in part by the Office of Nuclear Physics within the US DOE Office of Science, under Grant DE-FG02-89ER40531. The authors are grateful to the following colleagues for insightful advice and discussions: J. Auvinen, H. Caines, D. Cebra, V. Dexheimer, M. Lisa, W. Llope, G. Odyniec, Y. Pandit, H. Petersen, A. Poskanzer, H. G. Ritter, J. Steinheimer, H. Stöcker, A. H. Tang, S. Voloshin, G. Wang, N. Xu, and Z. Xu.