Korteveg-de Vries solitons in a cold quark-gluon plasma

The relativistic heavy ion program developed at RHIC and now at LHC motivated a deeper study of the properties of the quark-gluon plasma (QGP) and, in particular, the study of perturbations in this kind of plasma. We are interested on the time evolution of perturbations in the baryon and energy densities. If a localized pulse in baryon density could propagate throughout the QGP for long distances preserving its shape and without loosing localization, this could have interesting consequences for relativistic heavy ion physics and for astrophysics. A mathematical way to prove that this can happen is to derive (under certain conditions) from the hydrodynamical equations of the QGP a Korteveg-de Vries (KdV) equation. The solution of this equation describes the propagation of a KdV soliton. The derivation of the KdV equation depends crucially on the equation of state (EOS) of the QGP. The use of the simple MIT bag model EOS does not lead to KdV solitons. Recently we have developed an EOS for the QGP which includes both perturbative and nonperturbative corrections to the MIT one and is still simple enough to allow for analytical manipulations. With this EOS we were able to derive a KdV equation for the cold QGP.

Nonlinear waves in a quark gluon plasma

We study the propagation of perturbations in the quark gluon plasma. This subject has been addressed in other works and in most of the theoretical descriptions of this phenomenon the hydrodynamic equations have been linearized for simplicity. We propose an alternative approach, also based on hydrodynamics but taking into account the nonlinear terms of the equations. We show that these terms may lead to localized waves or even solitons. We use a simple equation of state for the QGP and expand the hydrodynamic equations around equilibrium configurations. The resulting differential equations describe the propagation of perturbations in the energy density. We solve them numerically and find that localized perturbations can propagate for long distances in the plasma. Under certain conditions our solutions mimic the propagation of Korteweg-de Vries solitons.

BREAKING WAVES IN AN IDEAL QUARK GLUON PLASMA

We study the propagation of perturbations in the energy density in a quark gluon plasma. Expanding the Euler and continuity equations of relativistic hydrodynamics around equilibrium configurations we obtain a nonlinear differential equation called the breaking wave equation. We solve it numerically and follow the time-evolution of initially localized pulses. We find that, quite unexpectedly, these pulses live for a very long time (compared to the reaction time-scales) before breaking. In practice, they mimick the Korteweg-de Vries solitons. Their existence may have some observable consequences.

Beam-energy and system-size dependence of dynamical net charge fluctuations

We present measurements of net charge fluctuations in Au+Au collisions at s(NN)=19.6, 62.4, 130, and 200 GeV, Cu+Cu collisions at s(NN)=62.4 and 200 GeV, and p+p collisions at s=200 GeV using the dynamical net charge fluctuations measure nu(+-,dyn). We observe that the dynamical fluctuations are nonzero at all energies and exhibit a modest dependence on beam energy. A weak system size dependence is also observed. We examine the collision centrality dependence of the net charge fluctuations and find that dynamical net charge fluctuations violate 1/N(ch) scaling but display approximate 1/N(part) scaling. We also study the azimuthal and rapidity dependence of the net charge correlation strength and observe strong dependence on the azimuthal angular range and pseudorapidity widths integrated to measure the correlation.

Spin alignment measurements of the K*(0)(892) and phi(1020) vector mesons in heavy ion collisions at root S(NN)=200 GeV

We present the first spin alignment measurements for the K*(0)(892) and phi(1020) vector mesons produced at midrapidity with transverse momenta up to 5 GeV/c at root s(NN) = 200 GeV at RHIC. The diagonal spin-density matrix elements with respect to the reaction plane in Au+Au collisions are rho(00) = 0.32 +/- 0.04 (stat) +/- 0.09 (syst) for the K*(0) (0.8 < p(T) < 5.0 GeV/c) and rho(00) = 0.34 +/- 0.02 (stat) +/- 0.03 (syst) for the phi (0.4 < p(T) < 5.0 GeV/c) and are constant with transverse momentum and collision centrality. The data are consistent with the unpolarized expectation of 1/3 and thus no evidence is found for the transfer of the orbital angular momentum of the colliding system to the vector-meson spins. Spin alignments for K(*0) and phi in Au+Au collisions were also measured with respect to the particle's production plane. The phi result, rho(00) = 0.41 +/- 0.02 (stat) +/- 0.04 (syst), is consistent with that in p+p collisions, rho(00) = 0.39 +/- 0.03 (stat) +/- 0.06 (syst), also measured in this work. The measurements thus constrain the possible size of polarization phenomena in the production dynamics of vector mesons.

Centrality dependence of charged hadron and strange hadron elliptic flow from root s(NN)=200 GeVAu+Au collisions

We present STAR results on the elliptic flow upsilon(2) Of charged hadrons, strange and multistrange particles from,root s(NN) = 200 GeV Au+Au collisions at the BNL Relativistic Heavy Ion Collider (RHIC). The detailed study of the centrality dependence of upsilon(2) over a broad transverse momentum range is presented. Comparisons of different analysis methods are made in order to estimate systematic uncertainties. To discuss the nonflow effect, we have performed the first analysis Of upsilon(2) with the Lee-Yang zero method for K(S)(0) and A. In the relatively low PT region, P(T) <= 2 GeV/c, a scaling with m(T) - m is observed for identified hadrons in each centrality bin studied. However, we do not observe nu 2(p(T))) scaled by the participant eccentricity to be independent of centrality. At higher PT, 2 1 <= PT <= 6 GeV/c, V2 scales with quark number for all hadrons studied. For the multistrange hadron Omega, which does not suffer appreciable hadronic interactions, the values of upsilon(2) are consistent with both m(T) - m scaling at low p(T) and number-of-quark scaling at intermediate p(T). As a function ofcollision centrality, an increase of p(T)-integrated upsilon(2) scaled by the participant eccentricity has been observed, indicating a stronger collective flow in more central Au+Au collisions.

High p(T) direct photon and pi(0) triggered azimuthal jet correlations and measurement of k(T) for isolated direct photons in p plus p collisions at root s=200 GeV

Correlations of charged hadrons of 1< p(T) < 10 Gev/c with high pT direct photons and pi(0) mesons in the range 5< p(T) < 15 Gev/c are used to study jet fragmentation in the gamma + jet and dijet channels, respectively. The magnitude of the partonic transverse momentum, k(T), is obtained by comparing to a model incorporating a Gaussian kT smearing. The sensitivity of the associated charged hadron spectra to the underlying fragmentation function is tested and the data are compared to calculations using recent global fit results. The shape of the direct photon-associated hadron spectrum as well as its charge asymmetry are found to be consistent with a sample dominated by quark-gluon Compton scattering. No significant evidence of fragmentation photon correlated production is observed within experimental uncertainties.

Azimuthal Anisotropy of pi(0) Production in Au plus Au Collisions at root s(NN)=200 GeV: Path-Length Dependence of Jet Quenching and the Role of Initial Geometry

We have measured the azimuthal anisotropy of pi(0) production for 1 < p(T) < 18 GeV/c for Au + Au collisions at root s(NN) = 200 GeV. The observed anisotropy shows a gradual decrease for 3 less than or similar to p(T) less than or similar to 7-10 GeV/c, but remains positive beyond 10 GeV/c. The magnitude of this anisotropy is underpredicted, up to at least similar to 10 GeV/c, by current perturbative QCD (PQCD) energy-loss model calculations. An estimate of the increase in anisotropy expected from initial-geometry modification due to gluon saturation effects and fluctuations is insufficient to account for this discrepancy. Calculations that implement a path-length dependence steeper than what is implied by current PQCD energy-loss models show reasonable agreement with the data.

Elliptic and Hexadecapole Flow of Charged Hadrons in Au plus Au Collisions at root s(NN)=200 GeV

Differential measurements of the elliptic (upsilon(2)) and hexadecapole (upsilon(4)) Fourier flow coefficients are reported for charged hadrons as a function of transverse momentum (p(T)) and collision centrality or number of participant nucleons (N(part)) for Au + Au collisions at root s(NN) = 200 GeV/ The upsilon(2,4) measurements at pseudorapidity vertical bar eta vertical bar <= 0.35, obtained with four separate reaction-plane detectors positioned in the range 1.0 < vertical bar eta vertical bar < 3.9, show good agreement, indicating the absence of significant Delta eta-dependent nonflow correlations. Sizable values for upsilon(4)(p(T)) are observed with a ratio upsilon(4)(p(T), N(part))/upsilon(2)(2)(p(T), N(part)) approximate to 0.8 for 50 less than or similar to N(part) less than or similar to 200, which is compatible with the combined effects of a finite viscosity and initial eccentricity fluctuations. For N(part) greater than or similar to 200 this ratio increases up to 1.7 in the most central collisions.

Enhanced Production of Direct Photons in Au plus Au Collisions at root s(NN)=200 GeV and Implications for the Initial Temperature

The production of e(+)e(-) pairs for m(e+e-) < 0.3 GeV/c(2) and 1< p(T) < 5 GeV/c is measured in p + p and Au + Au collisions at root s(NN) = 200 GeV. An enhanced yield above hadronic sources is observed. Treating the excess as photon internal conversions, the invariant yield of direct photons is deduced. In central Au + Au collisions, the excess of the direct photon yield over p + p is exponential in transverse momentum, with an inverse slope T = 221 +/- 19(stat) +/- 19(syst) MeV. Hydrodynamical models with initial temperatures ranging from T(init) similar to 300-600 MeV at times of similar to 0.6-0.15 fm/c after the collision are in qualitative agreement with the data. Lattice QCD predicts a phase transition to quark gluon plasma at similar to 170 MeV.

Detailed measurement of the e(+)e(-) pair continuum in p plus p and Au plus Au collisions at root s(NN)=200 GeV and implications for direct photon production

PHENIX has measured the e(+)e(-) pair continuum in root s(NN) = 200 GeV Au+Au and p+p collisions over a wide range of mass and transverse momenta. The e(+)e(-) yield is compared to the expectations from hadronic sources, based on PHENIX measurements. In the intermediate-mass region, between the masses of the phi and the J/psi meson, the yield is consistent with expectations from correlated c (c) over bar production, although other mechanisms are not ruled out. In the low-mass region, below the phi, the p+p inclusive mass spectrum is well described by known contributions from light meson decays. In contrast, the Au+Au minimum bias inclusive mass spectrum in this region shows an enhancement by a factor of 4.7 +/- 0.4(stat) +/- 1.5(syst) +/- 0.9(model). At low mass (m(ee) < 0.3 GeV/c(2)) and high p(T) (1 < p(T) < 5 GeV/c) an enhanced e(+)e(-) pair yield is observed that is consistent with production of virtual direct photons. This excess is used to infer the yield of real direct photons. In central Au+Au collisions, the excess of the direct photon yield over the p+p is exponential in p(T), with inverse slope T = 221 +/- 19(stat) +/- 19(syst) MeV. Hydrodynamical models with initial temperatures ranging from T(init) similar or equal to 300-600 MeV at times of 0.6-0.15 fm/c after the collision are in qualitative agreement with the direct photon data in Au+Au. For low p(T) < 1 GeV/c the low-mass region shows a further significant enhancement that increases with centrality and has an inverse slope of T similar or equal to 100 MeV. Theoretical models underpredict the low-mass, low-p(T) enhancement.

K*(0) production in Cu plus Cu and Au plus Au collisions at root s(NN)=62.4 GeV and 200 GeV

We report on K*(0) production at midrapidity in Au + Au and Cu + Cu collisions at root s(NN) = 62.4 and 200 GeV collected by the Solenoid Tracker at the Relativistic Heavy Ion Collider detector. The K*(0) is reconstructed via the hadronic decays K*(0) -> K(+)pi(-) and (K*(0)) over bar -> K(+)pi(-). Transverse momentum, p(T), spectra are measured over a range of p(T) extending from 0.2 GeV/c up to 5 GeV/c. The center-of-mass energy and system size dependence of the rapidity density, dN/dy, and the average transverse momentum, < p(T)>, are presented. The measured N(K*(0))/N(K) and N(phi)/N(K*(0)) ratios favor the dominance of rescattering of decay daughters of K*(0) over the hadronic regeneration for the K*(0) production. In the intermediate p(T) region (2.0 < p(T) < 4.0 GeV/c), the elliptic flow parameter, v(2), and the nuclear modification factor, R(CP), agree with the expectations from the quark coalescence model of particle production.

Parton energy loss in heavy-ion collisions via direct-photon and charged-particle azimuthal correlations

Charged-particle spectra associated with direct photon (gamma(dir)) and pi(0) are measured in p + p and Au + Au collisions at center-of-mass energy root(S)(NN) = 200 GeV with the STAR detector at the Relativistic Heavy Ion Collider. A shower-shape analysis is used to partially discriminate between gamma(dir) and pi(0). Assuming no associated charged particles in the gamma(dir) direction ( near side) and small contribution from fragmentation photons (gamma(frag)), the associated charged-particle yields opposite to gamma(dir) (away side) are extracted. In central Au + Au collisions, the charged-particle yields at midrapidity (vertical bar eta vertical bar < 1) and high transverse momentum (3 < (assoc)(PT) < 16 GeV/c) associated with gamma(dir) and pi(0) (vertical bar eta vertical bar < 0.9, 8 < (trig)(PT) < 16 GeV/c) are suppressed by a factor of 3-5 compared with p + p collisions. The observed suppression of the associated charged particles is similar for gamma(dir) and pi(0) and independent of the gamma(dir) energy within uncertainties. These measurements indicate that, in the kinematic range covered and within our current experimental uncertainties, the parton energy loss shows no sensitivity to the parton initial energy, path length, or color charge.

Balance functions from Au+Au, d+Au, and p+p collisions at root s(NN)=200 GeV

Balance functions have been measured for charged-particle pairs, identified charged-pion pairs, and identified charged-kaon pairs in Au + Au, d + Au, and p + p collisions at root s(NN) = 200 GeV at the Relativistic Heavy Ion Collider using the STAR detector. These balance functions are presented in terms of relative pseudorapidity, Delta eta, relative rapidity, Delta y, relative azimuthal angle, Delta phi, and invariant relative momentum, q(inv). For charged-particle pairs, the width of the balance function in terms of Delta eta scales smoothly with the number of participating nucleons, while HIJING and UrQMD model calculations show no dependence on centrality or system size. For charged-particle and charged-pion pairs, the balance functions widths in terms of Delta eta and Delta y are narrower in central Au + Au collisions than in peripheral collisions. The width for central collisions is consistent with thermal blast-wave models where the balancing charges are highly correlated in coordinate space at breakup. This strong correlation might be explained by either delayed hadronization or limited diffusion during the reaction. Furthermore, the narrowing trend is consistent with the lower kinetic temperatures inherent to more central collisions. In contrast, the width of the balance function for charged-kaon pairs in terms of Delta y shows little centrality dependence, which may signal a different production mechanism for kaons. The widths of the balance functions for charged pions and kaons in terms of q(inv) narrow in central collisions compared to peripheral collisions, which may be driven by the change in the kinetic temperature.

Observation of charge-dependent azimuthal correlations and possible local strong parity violation in heavy-ion collisions

Parity (P)-odd domains, corresponding to nontrivial topological solutions of the QCD vacuum, might be created during relativistic heavy-ion collisions. These domains are predicted to lead to charge separation of quarks along the orbital momentum of the system created in noncentral collisions. To study this effect, we investigate a three-particle mixed-harmonics azimuthal correlator which is a P-even observable, but directly sensitive to the charge-separation effect. We report measurements of this observable using the STAR detector in Au + Au and Cu + Cu collisions at root s(NN) = 200 and 62 GeV. The results are presented as a function of collision centrality, particle separation in rapidity, and particle transverse momentum. A signal consistent with several of the theoretical expectations is detected in all four data sets. We compare our results to the predictions of existing event generators and discuss in detail possible contributions from other effects that are not related to P violation.