For central heavy ion collisions at the RHIC energy, the entropy per unit rapidity dS/dy at freeze-out is extracted with minimal model dependence from available experimental measurements of particle yields, spectra, and source sizes estimated from two-particle interferometry. The extracted entropy rapidity density is consistent with lattice gauge theory results for a thermalized quark-gluon plasma with an energy density estimated from transverse energy production at RHIC. PACS numbers: 12.38.Mh, 24.85.+p, Lattice quantum chromodynamics (QCD) calculations indicate [1] that at zero net baryon density at a critical temperature of T c ∼ 170 MeV and energy density ǫ c ∼ 1 GeV/fm 3 a color deconfined and a chirally restored quark-gluon plasma (QGP) phase is formed. These energy densities are attained in relativistic heavy ion collisions, where it is believed that a dense system of quarks and gluons is created which undergoes rapid collective expansion before the partons hadronize and eventually decouple.The ongoing experimental program at RHIC is devoted to the search for the signals of QGP. Unfortunately, measurements are confined to final-state particles which are mostly hadrons. One interesting quantity that may provide valuable insight into the state of the matter in the early stages of the collision is the net entropy which is roughly conserved between initial thermalization and freeze-out [2,3]. After freeze-out, when particles are freestreaming, the entropy is fixed. Entropy would also be conserved during the expansion stage if the mean free paths were very short. Viscous effects, shock waves, and the decoupling process can only increase the entropy. Thus, the entropy from the final state provides an upper bound for the entropy of the initial state.Compared to a pion gas at the same energy density, the QGP is a high-entropy state as can be understood by considering the equation,where ∆S is the change in entropy resulting from the addition of an amount of heat Q into a fixed volume at temperature T . Lattice calculations show that the transition is nearly a first order with a distinct peak in the specific heat near T c . In this region, the temperature stays relatively fixed while the energy density changes substantially. Thus, the existence of a latent heat keeps the temperature lower than a pion gas at the same energy density, which results in a higher entropy as can be seen from Eq. (1). At high energy densities, the temperature of a QGP remains lower than that of the pion gas due to the disparity in the number of degrees of freedom.Asymptotically, the energy density of a QGP approaches that of a massless gas of quarks and hadrons,where N B and N F are the number of Bosonic and Fermionic degrees of freedom. Whereas a pion gas has 3 bosonic degrees of freedom, a QGP has 47.5 effective degrees of freedom if one considers gluons and the three light quark flavors. Thus, relative to a pion gas, the QGP has a lower temperature and higher entropy. If the partonic degrees of freedom were not fully liberated at RHIC, o...