We present measurements of $e^+e^-$ production at midrapidity in Au$+$Au collisions at $\sqrt{s_{_{NN}}}$ = 200 GeV. The invariant yield is studied within the PHENIX detector acceptance over a wide range of mass ($m_{ee} <$ 5 GeV/$c^2$) and pair transverse momentum ($p_T$ $<$ 5 GeV/$c$), for minimum bias and for five centrality classes. The \ee yield is compared to the expectations from known sources. In the low-mass region ($m_{ee}=0.30$--0.76 GeV/$c^2$) there is an enhancement that increases with centrality and is distributed over the entire pair \pt range measured. It is significantly smaller than previously reported by the PHENIX experiment and amounts to $2.3\pm0.4({\rm stat})\pm0.4({\rm syst})\pm0.2^{\rm model}$ or to $1.7\pm0.3({\rm stat})\pm0.3({\rm syst})\pm0.2^{\rm model}$ for minimum bias collisions when the open-heavy-flavor contribution is calculated with {\sc pythia} or {\sc mc@nlo}, respectively. The inclusive mass and $p_T$ distributions as well as the centrality dependence are well reproduced by model calculations where the enhancement mainly originates from the melting of the $\rho$ meson resonance as the system approaches chiral symmetry restoration. In the intermediate-mass region ($m_{ee}$ = 1.2--2.8 GeV/$c^2$), the data hint at a significant contribution in addition to the yield from the semileptonic decays of heavy-flavor mesons.

The PHENIX experiment has measured $\phi$ meson production in $d$$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV using the dimuon and dielectron decay channels. The $\phi$ meson is measured in the forward (backward) $d$-going (Au-going) direction, $1.2<y<2.2$ ($-2.2<y<-1.2$) in the transverse-momentum ($p_T$) range from 1--7 GeV/$c$, and at midrapidity $|y|<0.35$ in the $p_T$ range below 7 GeV/$c$. The $\phi$ meson invariant yields and nuclear-modification factors as a function of $p_T$, rapidity, and centrality are reported. An enhancement of $\phi$ meson production is observed in the Au-going direction, while suppression is seen in the $d$-going direction, and no modification is observed at midrapidity relative to the yield in $p$$+$$p$ collisions scaled by the number of binary collisions. Similar behavior was previously observed for inclusive charged hadrons and open heavy flavor indicating similar cold-nuclear-matter effects.

The measurement of the direct-photon spectrum from Au$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV is presented by the PHENIX collaboration using the external-photon-conversion technique for 0%--93% central collisions in a transverse-momentum ($p_T$) range of 0.8--10 GeV/$c$. An excess of direct photons, above prompt-photon production from hard-scattering processes, is observed for $p_T<6$ GeV/$c$. Nonprompt direct photons are measured by subtracting the prompt component, which is estimated as $N_{\rm coll}$-scaled direct photons from $p$ $+$ $p$ collisions at 200 GeV, from the direct-photon spectrum. Results are obtained for $0.8<p_T<6.0$ GeV/$c$ and suggest that the spectrum has an increasing inverse slope from ${\approx}0.2$ to 0.4 GeV/$c$ with increasing $p_T$, which indicates a possible sensitivity of the measurement to photons from earlier stages of the evolution of the collision. In addition, like the direct-photon production, the $p_T$-integrated nonprompt direct-photon yields also follow a power-law scaling behavior as a function of collision-system size. The exponent, $\alpha$, for the nonprompt component is found to be consistent with 1.1 with no apparent $p_T$ dependence.

The invariant yields for $J/\psi$ production at forward rapidity $(1.2<|y|<2.2)$ in U$+$U collisions at $\sqrt{s_{_{NN}}}$=193 GeV have been measured as a function of collision centrality. The invariant yields and nuclear-modification factor $R_{AA}$ are presented and compared with those from Au$+$Au collisions in the same rapidity range. Additionally, the direct ratio of the invariant yields from U$+$U and Au$+$Au collisions within the same centrality class is presented, and used to investigate the role of $c\bar{c}$ coalescence. Two different parameterizations of the deformed Woods-Saxon distribution were used in Glauber calculations to determine the values of the number of nucleon-nucleon collisions in each centrality class, $N_{\rm coll}$, and these were found to give significantly different $N_{\rm coll}$ values. Results using $N_{\rm coll}$ values from both deformed Woods-Saxon distributions are presented. The measured ratios show that the $J/\psi$ suppression, relative to binary collision scaling, is similar in U$+$U and Au$+$Au for peripheral and midcentral collisions, but that $J/\psi$ show less suppression for the most central U$+$U collisions. The results are consistent with a picture in which, for central collisions, increase in the $J/\psi$ yield due to $c\bar{c}$ coalescence becomes more important than the decrease in yield due to increased energy density. For midcentral collisions, the conclusions about the balance between $c\bar{c}$ coalescence and suppression depend on which deformed Woods-Saxon distribution is used to determine $N_{\rm coll}$.

The PHENIX collaboration presents first measurements of low-momentum ($0.4<p_T<3$ GeV/$c$) direct-photon yields from Au$+$Au collisions at $\sqrt{s_{_{NN}}}$=39 and 62.4 GeV. For both beam energies the direct-photon yields are substantially enhanced with respect to expectations from prompt processes, similar to the yields observed in Au$+$Au collisions at $\sqrt{s_{_{NN}}}$=200. Analyzing the photon yield as a function of the experimental observable $dN_{\rm ch}/d\eta$ reveals that the low-momentum ($>$1\,GeV/$c$) direct-photon yield $dN_{\gamma}^{\rm dir}/d\eta$ is a smooth function of $dN_{\rm ch}/d\eta$ and can be well described as proportional to $(dN_{\rm ch}/d\eta)^\alpha$ with $\alpha{\sim}$1.25. This new scaling behavior holds for a wide range of beam energies at the Relativistic Heavy Ion Collider and Large Hadron Collider, for centrality selected samples, as well as for different, $A$$+$$A$ collision systems. At a given beam energy the scaling also holds for high $p_T$ ($>5$\,GeV/$c$) but when results from different collision energies are compared, an additional $\sqrt{s_{_{NN}}}$-dependent multiplicative factor is needed to describe the integrated-direct-photon yield.

We have measured direct photons for $p_T<5~$GeV/$c$ in minimum bias and 0\%--40\% most central events at midrapidity for Cu$+$Cu collisions at $\sqrt{s_{_{NN}}}=200$ GeV. The $e^{+}e^{-}$ contribution from quasi-real direct virtual photons has been determined as an excess over the known hadronic contributions in the $e^{+}e^{-}$ mass distribution. A clear enhancement of photons over the binary scaled $p$$+$$p$ fit is observed for $p_T<4$ GeV/$c$ in Cu$+$Cu data. The $p_T$ spectra are consistent with the Au$+$Au data covering a similar number of participants. The inverse slopes of the exponential fits to the excess after subtraction of the $p$$+$$p$ baseline are 285$\pm$53(stat)$\pm$57(syst)~MeV/$c$ and 333$\pm$72(stat)$\pm$45(syst)~MeV/$c$ for minimum bias and 0\%--40\% most central events, respectively. The rapidity density, $dN/dy$, of photons demonstrates the same power law as a function of $dN_{\rm ch}/d\eta$ observed in Au$+$Au at the same collision energy.

We present a detailed measurement of charged two-pion correlation functions in 0%-30% centrality $\sqrt{s_{_{NN}}}=200$ GeV Au$+$Au collisions by the PHENIX experiment at the Relativistic Heavy Ion Collider. The data are well described by Bose-Einstein correlation functions stemming from L\'evy-stable source distributions. Using a fine transverse momentum binning, we extract the correlation strength parameter $\lambda$, the L\'evy index of stability $\alpha$ and the L\'evy length scale parameter $R$ as a function of average transverse mass of the pair $m_T$. We find that the positively and the negatively charged pion pairs yield consistent results, and their correlation functions are represented, within uncertainties, by the same L\'evy-stable source functions. The $\lambda(m_T)$ measurements indicate a decrease of the strength of the correlations at low $m_T$. The L\'evy length scale parameter $R(m_T)$ decreases with increasing $m_T$, following a hydrodynamically predicted type of scaling behavior. The values of the L\'evy index of stability $\alpha$ are found to be significantly lower than the Gaussian case of $\alpha=2$, but also significantly larger than the conjectured value that may characterize the critical point of a second-order quark-hadron phase transition.

The PHENIX experiment at the Relativistic Heavy Ion Collider has measured low mass vector meson, $\omega$, $\rho$, and $\phi$, production through the dimuon decay channel at forward rapidity ($1.2<|y|<2.2$) in $p$$+$$p$ collisions at $\sqrt{s}=200$ GeV. The differential cross sections for these mesons are measured as a function of both $p_T$ and rapidity. We also report the integrated differential cross sections over $1<p_T<7$ GeV/$c$ and $1.2<|y|<2.2$: $d\sigma/dy(\omega+\rho\rightarrow\mu\mu) = 80 \pm 6 \mbox{(stat)} \pm 12 \mbox{(syst)}$ nb and $d\sigma/dy(\phi\rightarrow\mu\mu) = 27 \pm 3 \mbox{(stat)} \pm 4 \mbox{(syst)}$ nb. These results are compared with midrapidity measurements and calculations.

The measurement of direct photons from Au$+$Au collisions at $\sqrt{s_{_{NN}}}=39$ and 62.4 GeV in the transverse-momentum range $0.4<p_T<3$ Gev/$c$ is presented by the PHENIX collaboration at the Relativistic Heavy Ion Collider. A significant direct-photon yield is observed in both collision systems. A universal scaling is observed when the direct-photon $p_T$ spectra for different center-of-mass energies and for different centrality selections at $\sqrt{s_{_{NN}}}=62.4$ GeV is scaled with $(dN_{\rm ch}/d\eta)^{\alpha}$ for $\alpha=1.21{\pm}0.04$. This scaling also holds true for direct-photon spectra from Au$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV measured earlier by PHENIX, as well as the spectra from Pb$+$Pb at $\sqrt{s_{_{NN}}}=2760$ GeV published by ALICE. The scaling power $\alpha$ seems to be independent of $p_T$, center of mass energy, and collision centrality. The spectra from different collision energies have a similar shape up to $p_T$ of 2 GeV/$c$. The spectra have a local inverse slope $T_{\rm eff}$ increasing with $p_T$ of $0.174\pm0.018$ GeV/$c$ in the range $0.4<p_T<1.3$ GeV/$c$ and increasing to $0.289\pm0.024$ GeV/$c$ for $0.9<p_T<2.1$ GeV/$c$. The observed similarity of low-$p_T$ direct-photon production from $\sqrt{s_{_{NN}}}= 39$ to 2760 GeV suggests a common source of direct photons for the different collision energies and event centrality selections, and suggests a comparable space-time evolution of direct-photon emission.

Production of $\pi^0$ and $\eta$ mesons has been measured at midrapidity in Cu$+$Au collisions at $\sqrt{s_{_{NN}}}$=200 GeV. Measurements were performed in $\pi^0(\eta)\rightarrow\gamma\gamma$ decay channel in the 1(2)-20 GeV/$c$ transverse momentum range. A strong suppression is observed for $\pi^0$ and $\eta$ meson production at high transverse momentum in central Cu$+$Au collisions relative to the $p$$+$$p$ results scaled by the number of nucleon-nucleon collisions. In central collisions the suppression is similar to Au$+$Au with comparable nuclear overlap. The $\eta/\pi^0$ ratio measured as a function of transverse momentum is consistent with $m_T$-scaling parameterization down to $p_T=$2 GeV/$c$, its asymptotic value is constant and consistent with Au$+$Au and $p$$+$$p$ and does not show any significant dependence on collision centrality. Similar results were obtained in hadron-hadron, hadron-nucleus, and nucleus-nucleus collisions as well as in $e^+e^-$ collisions in a range of collision energies $\sqrt{s_{_{NN}}}=$3--1800 GeV. This suggests that the quark-gluon-plasma medium produced in Cu$+$Cu collisions either does not affect the jet fragmentation into light mesons or it affects the $\pi^0$ and $\eta$ the same way.