Hydro-kinetic approach to relativistic heavy ion collisions

We develop a combined hydro-kinetic approach which incorporates a hydrodynamical expansion of the systems formed in A + A collisions and their dynamical decoupling described by escape probabilities. The method corresponds to a generalized relaxation time (tau(rel)) approximation for the Boltzmann equation applied to inhomogeneous expanding systems; at small tau(rel) it also allows one to catch the viscous effects in hadronic component-hadron-resonance gas. We demonstrate how the approximation of sudden freeze-out can be obtained within this dynamical picture of continuous emission and find that hypersurfaces, corresponding to a sharp freeze-out limit, are momentum dependent. The pion m(T) spectra are computed in the developed hydro-kinetic model, and compared with those obtained from ideal hydrodynamics with the Cooper-Frye isothermal prescription. Our results indicate that there does not exist a universal freeze-out temperature for pions with different momenta, and support an earlier decoupling of higher p(T) particles. By performing numerical simulations for various initial conditions and equations of state we identify several characteristic features of the bulk QCD matter evolution preferred in view of the current analysis of heavy ion collisions at RHIC energies.

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.

gamma cross section in p plus p collisions at root s=200 GeV

We report on a measurement of the gamma(1S + 2S + 3S) -> e(+)e(-) cross section at midrapidity in p + p collisions at root s = 200 GeV. We find the cross section to be 114 +/- 38(stat + fit)(-24)(+23)(syst) pb. Perturbative QCD calculations at next-to-leading order in the color evaporation model are in agreement with our measurement, while calculations in the color singlet model underestimate it by 2 sigma. Our result is consistent with the trend seen in world data as a function of the center-of-mass energy of the collision and extends the availability of gamma data to RHIC energies. The dielectron continuum in the invariant-mass range near the gamma is also studied to obtain a combined yield of e(+)e(-) pairs from the sum of the Drell-Yan process and b-(b) over bar production.

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.

Charged and strange hadron elliptic flow in Cu plus Cu collisions at root s(NN)=62.4 and 200 GeV

We present the results of an elliptic flow, v(2), analysis of Cu + Cu collisions recorded with the solenoidal tracker detector (STAR) at the BNL Relativistic Heavy Ion Collider at root s(NN) = 62.4 and 200 GeV. Elliptic flow as a function of transverse momentum, v(2)(p(T)), is reported for different collision centralities for charged hadrons h(+/-) and strangeness-ontaining hadrons K(S)(0), Lambda, Xi, and phi in the midrapidity region vertical bar eta vertical bar < 1.0. Significant reduction in systematic uncertainty of the measurement due to nonflow effects has been achieved by correlating particles at midrapidity, vertical bar eta vertical bar < 1.0, with those at forward rapidity, 2.5 < vertical bar eta vertical bar < 4.0. We also present azimuthal correlations in p + p collisions at root s = 200 GeV to help in estimating nonflow effects. To study the system-size dependence of elliptic flow, we present a detailed comparison with previously published results from Au + Au collisions at root s(NN) = 200 GeV. We observe that v(2)(p(T)) of strange hadrons has similar scaling properties as were first observed in Au + Au collisions, that is, (i) at low transverse momenta, p(T) < 2 GeV/c, v(2) scales with transverse kinetic energy, m(T) - m, and (ii) at intermediate p(T), 2 < p(T) < 4 GeV/c, it scales with the number of constituent quarks, n(q.) We have found that ideal hydrodynamic calculations fail to reproduce the centrality dependence of v(2)(p(T)) for K(S)(0) and Lambda. Eccentricity scaled v(2) values, v(2)/epsilon, are larger in more central collisions, suggesting stronger collective flow develops in more central collisions. The comparison with Au + Au collisions, which go further in density, shows that v(2)/epsilon depends on the system size, that is, the number of participants N(part). This indicates that the ideal hydrodynamic limit is not reached in Cu + Cu collisions, presumably because the assumption of thermalization is not attained.

Identified particle production, azimuthal anisotropy, and interferometry measurements in Au plus Au collisions at root s(NN)=9.2 GeV

We present the first measurements of identified hadron production, azimuthal anisotropy, and pion interferometry from Au + Au collisions below the nominal injection energy at the BNL Relativistic Heavy-Ion Collider (RHIC) facility. The data were collected using the large acceptance solenoidal tracker at RHIC (STAR) detector at root s(NN) = 9.2 GeV from a test run of the collider in the year 2008. Midrapidity results on multiplicity density dN/dy in rapidity y, average transverse momentum < p(T)>, particle ratios, elliptic flow, and Hanbury-Brown-Twiss (HBT) radii are consistent with the corresponding results at similar root s(NN) from fixed-target experiments. Directed flow measurements are presented for both midrapidity and forward-rapidity regions. Furthermore the collision centrality dependence of identified particle dN/dy, < p(T)>, and particle ratios are discussed. These results also demonstrate that the capabilities of the STAR detector, although optimized for root s(NN) = 200 GeV, are suitable for the proposed QCD critical-point search and exploration of the QCD phase diagram at RHIC.

Long range rapidity correlations and jet production in high energy nuclear collisions

The STAR Collaboration at the Relativistic Heavy Ion Collider presents a systematic study of high-transverse-momentum charged-di-hadron correlations at small azimuthal pair separation Delta phi in d+Au and central Au+Au collisions at s(NN)=200 GeV. Significant correlated yield for pairs with large longitudinal separation Delta eta is observed in central Au+Au collisions, in contrast to d+Au collisions. The associated yield distribution in Delta eta x Delta phi can be decomposed into a narrow jet-like peak at small angular separation which has a similar shape to that found in d+Au collisions, and a component that is narrow in Delta phi and depends only weakly on Delta eta, the ""ridge."" Using two systematically independent determinations of the background normalization and shape, finite ridge yield is found to persist for trigger p(t)>6 GeV/c, indicating that it is correlated with jet production. The transverse-momentum spectrum of hadrons comprising the ridge is found to be similar to that of bulk particle production in the measured range (2 < p(t)< 4 GeV/c).

J/psi production at high transverse momenta in p plus p and Cu plus Cu collisions at root s(NN)=200 GeV

The STAR Collaboration at the Relativistic Heavy Ion Collider presents measurements of J/psi e(+) e(-) at midrapidity and high transverse momentum (pT > 5 GeV/c) in p + p and central Cu + Cu collisions at root s(NN) = 200 GeV. The inclusive J/psi production cross section for Cu + Cu collisions is found to be consistent at high p(T) with the binary collision-scaled cross section for p + p collisions. At a confidence level of 97%, this is in contrast to a suppression of J/psi production observed at lower p(T). Azimuthal correlations of J/psi with charged hadrons in p + p collisions provide an estimate of the contribution of B-hadron decays to J/psi production of 13% +/- 5%.

Pion interferometry in Au plus Au and Cu plus Cu collisions at s(NN)=62.4 and 200 GeV

We present a systematic analysis of two-pion interferometry in Au+Au collisions at s(NN)=62.4 GeV and Cu+Cu collisions at s(NN)=62.4 and 200 GeV using the STAR detector at the Relativistic Heavy Ion Collider (RHIC). The multiplicity and transverse momentum dependences of the extracted correlation lengths (radii) are studied. The scaling with charged particle multiplicity of the apparent system volume at final interaction is studied for the RHIC energy domain. The multiplicity scaling of the measured correlation radii is found to be independent of colliding system and collision energy.

Pion femtoscopy in p plus p collisions at root s=200 GeV

The STAR Collaboration at the BNL Relativistic Heavy Ion Collider has measured two-pion correlation functions from p + p collisions at root s = 200 GeV. Spatial scales are extracted via a femtoscopic analysis of the correlations, though this analysis is complicated by the presence of strong nonfemtoscopic effects. Our results are put into the context of the world data set of femtoscopy in hadron-hadron collisions. We present the first direct comparison of femtoscopy in p + p and heavy ion collisions, under identical analysis and detector conditions.

EXACT AND QUASI-EXACT MODELS OF STRANGE STARS

We construct and compare in this work a variety of simple models for strange stars, namely, hypothetical self-bound objects made of a cold stable version of the quark-gluon plasma. Exact, quasi-exact and numerical models are examined to find the most economical description for these objects. A simple and successful parametrization of them is given in terms of the central density, and the differences among the models are explicitly shown and discussed. In particular, we present a model starting with a Gaussian ansatz for the density profile that provides a very accurate and almost complete analytical integration of the problem, modulo a small difference for one of the metric potentials.

Thermal photon production in Au plus Au collisions: Viscous corrections in two different hydrodynamic formalisms

We calculate the spectra of produced thermal photons in Au + Au collisions taking into account the nonequilibrium contribution to photon production due to finite shear viscosity. The evolution of the fireball is modeled by second-order as well as by divergence-type 2 + 1 dissipative hydrodynamics, both with an ideal equation of state and with one based on Lattice QCD that includes an analytical crossover. The spectrum calculated in the divergence-type theory is considerably enhanced with respect to the one calculated in the second-order theory, the difference being entirely due to differences in the viscous corrections to photon production. Our results show that the differences in hydrodynamic formalisms are an important source of uncertainty in the extraction of the value of eta/s from measured photon spectra. The uncertainty in the value of eta/s associated with different hydrodynamic models used to compute thermal photon spectra is larger than the one occurring in matching hadron elliptic flow to RHIC data. (C) 2010 Elsevier B.V. All rights reserved.

Sound waves and solitons in hot and dense nuclear matter

Assuming that nuclear matter can be treated as a perfect fluid, we study the propagation of perturbations in the baryon density. The equation of state is derived from a relativistic mean field model, which is a variant of the non-linear Walecka model. The expansion of the Euler and continuity equations of relativistic hydrodynamics around equilibrium configurations leads to differential equations for the density perturbation. We solve them numerically for linear and spherical perturbations and follow the propagation of the initial pulses. For linear perturbations we find single soliton solutions and solutions with one or more solitons followed by ""radiation"". Depending on the equation of state a strong damping may occur. We consider also the evolution of perturbations in a medium without dispersive effects. In this case we observe the formation and breaking of shock waves. We study all these equations also for matter at finite temperature. Our results may be relevant for the analysis of RHIC data. They suggest that the shock waves formed in the quark gluon plasma phase may survive and propagate in the hadronic phase. (C) 2009 Elseiver. B.V. All rights reserved.