The system created in non-central relativistic nucleus-nucleus collisions possesses large orbital angular momentum. Due to spin-orbit coupling, particles produced in such a system could become globally polarized along the direction of the system angular momentum. We present the results of Lambda and anti-Lambda hyperon global polarization measurements in Au+Au collisions at sqrt{s_NN}=62.4 GeV and 200 GeV performed with the STAR detector at RHIC. The observed global polarization of Lambda and anti-Lambda hyperons in the STAR acceptance is consistent with zero within the precision of the measurements. The obtained upper limit, |P_{Lambda,anti-Lambda}| <= 0.02, is compared to the theoretical values discussed recently in the literature.

Identified charged pion, kaon, and proton spectra are used to explore the system size dependence of bulk freeze-out properties in Cu+Cu collisions at $\sqrt{s_{NN}}$=200 and 62.4 GeV. The data are studied with hydrodynamically-motivated Blast-wave and statistical model frameworks in order to characterize the freeze-out properties of the system. The dependence of freeze-out parameters on beam energy and collision centrality is discussed. Using the existing results from Au+Au and $pp$ collisions, the dependence of freeze-out parameters on the system size is also explored. This multi-dimensional systematic study furthers our understanding of the QCD phase diagram revealing the importance of the initial geometrical overlap of the colliding ions. The analysis of Cu+Cu collisions, which expands the system size dependence studies from Au+Au data with detailed measurements in the smaller system, shows that the bulk freeze-out properties of charged particles studied here scale with the total charged particle multiplicity at mid-rapidity, suggesting the relevance of initial state effects.

We present a systematic analysis of two-pion interferometry in Au+Au collisions at $\sqrt{s_{\rm{NN}}}$ = 62.4 GeV and Cu+Cu collisions at $\sqrt{s_{\rm{NN}}}$ = 62.4 and 200 GeV using the STAR detector at RHIC. The multiplicity and transverse momentum dependences of the extracted femtoscopic radii are studied. The scaling of the apparent freeze-out volume with charged particle multiplicity is studied for the RHIC energy domain. The multiplicity scaling of the measured radii is found to be independent of colliding system and collision energy.

We present results for the measurement of $\phi$ meson production via its charged kaon decay channel $\phi \to K^+K^-$ in Au+Au collisions at $\sqrt{s_{_{NN}}}=62.4$, 130, 200 GeV, and in p+p and d+Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV from the STAR experiment at RHIC. The mid-rapidity ($|y|<0.5$) $\phi$ meson spectra in central Au+Au collisions are found to be well described by a single exponential distribution. On the other hand, the spectra from p+p, d+Au and peripheral Au+Au collisions show power-law tails at intermediate and high transverse momenta ($p_{T}$) and are described better by Levy distributions. The constant $\phi/K^-$ yield ratio vs. beam species, collision centrality and colliding energy is in contradiction with expectations from models having kaon coalescence as the dominant mechanism for $\phi$ production at RHIC. The $\Omega/\phi$ yield ratio as a function of $p_{T}$ is consistent with a model based on the recombination of thermal $s$ quarks up to $p_{T}\sim 4$ GeV/c, but disagrees at higher transverse momenta. The measured nuclear modification factor, $R_{dAu}$, for the $\phi$ meson increases above unity at intermediate $p_{T}$, similar to that for pions and protons, while $R_{AA}$ is suppressed due to jet quenching in central Au+Au collisions. Number of constituent quark scaling of both $R_{cp}$ and $v_{2}$ for the $\phi$ meson with respect to other hadrons in Au+Au collisions at $\sqrt{s_{_{NN}}}$=200 GeV at intermediate $p_{T}$ is observed. These observations support quark coalescence as being the dominant mechanism of hadronization in the intermediate $p_{T}$ region at RHIC.

New measurements of directed flow for charged hadrons, characterized by the Fourier coefficient \vone, are presented for transverse momenta $\mathrm{p_T}$, and centrality intervals in Au+Au collisions recorded by the STAR experiment for the center-of-mass energy range $\mathrm{\sqrt{s_{_{NN}}}} = 7.7 - 200$ GeV. The measurements underscore the importance of momentum conservation and the characteristic dependencies on $\mathrm{\sqrt{s_{_{NN}}}}$, centrality and $\mathrm{p_T}$ are consistent with the expectations of geometric fluctuations generated in the initial stages of the collision, acting in concert with a hydrodynamic-like expansion. The centrality and $\mathrm{p_T}$ dependencies of $\mathrm{v^{even}_{1}}$, as well as an observed similarity between its excitation function and that for $\mathrm{v_3}$, could serve as constraints for initial-state models. The $\mathrm{v^{even}_{1}}$ excitation function could also provide an important supplement to the flow measurements employed for precision extraction of the temperature dependence of the specific shear viscosity.

We report the energy dependence of mid-rapidity (anti-)deuteron production in Au+Au collisions at $\sqrt{s_\text{NN}} =\ $7.7, 11.5, 14.5, 19.6, 27, 39, 62.4, and 200 GeV, measured by the STAR experiment at RHIC. The yield of deuterons is found to be well described by the thermal model. The collision energy, centrality, and transverse momentum dependence of the coalescence parameter $B_2$ are discussed. We find that the values of $B_2$ for anti-deuterons are systematically lower than those for deuterons, indicating that the correlation volume of anti-baryons is larger than that of baryons at $\sqrt{s_\text{NN}}$ from 19.6 to 39 GeV. In addition, values of $B_2$ are found to vary with collision energy and show a broad minimum around $\sqrt{s_\text{NN}}=\ $20 to 40 GeV, which might imply a change of the equation of state of the medium in these collisions.

The extreme temperatures and energy densities generated by ultra-relativistic collisions between heavy nuclei produce a state of matter with surprising fluid properties. Non-central collisions have angular momentum on the order of 1000$\hbar$, and the resulting fluid may have a strong vortical structure that must be understood to properly describe the fluid. It is also of particular interest because the restoration of fundamental symmetries of quantum chromodynamics is expected to produce novel physical effects in the presence of strong vorticity. However, no experimental indications of fluid vorticity in heavy ion collisions have so far been found. Here we present the first measurement of an alignment between the angular momentum of a non-central collision and the spin of emitted particles, revealing that the fluid produced in heavy ion collisions is by far the most vortical system ever observed. We find that $\Lambda$ and $\overline{\Lambda}$ hyperons show a positive polarization of the order of a few percent, consistent with some hydrodynamic predictions. A previous measurement that reported a null result at higher collision energies is seen to be consistent with the trend of our new observations, though with larger statistical uncertainties. These data provide the first experimental access to the vortical structure of the "perfect fluid" created in a heavy ion collision. They should prove valuable in the development of hydrodynamic models that quantitatively connect observations to the theory of the Strong Force. Our results extend the recent discovery of hydrodynamic spin alignment to the subatomic realm.

Balance functions have been measured in terms of relative pseudorapidity ($\Delta \eta$) for charged particle pairs at the Relativistic Heavy-Ion Collider (RHIC) from Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 7.7 GeV to 200 GeV using the STAR detector. These results are compared with balance functions measured at the Large Hadron Collider (LHC) from Pb+Pb collisions at $\sqrt{s_{\rm NN}}$ = 2.76 TeV by the ALICE Collaboration. The width of the balance function decreases as the collisions become more central and as the beam energy is increased. In contrast, the widths of the balance functions calculated using shuffled events show little dependence on centrality or beam energy and are larger than the observed widths. Balance function widths calculated using events generated by UrQMD are wider than the measured widths in central collisions and show little centrality dependence. The measured widths of the balance functions in central collisions are consistent with the delayed hadronization of a deconfined quark gluon plasma (QGP). The narrowing of the balance function in central collisions at $\sqrt{s_{\rm NN}}$ = 7.7 GeV implies that a QGP is still being created at this relatively low energy.

Density fluctuations near the QCD critical point can be probed via an intermittency analysis in relativistic heavy-ion collisions. We report the first measurement of intermittency in Au$+$Au collisions at $\sqrt{s_\mathrm{_{NN}}}$ = 7.7-200 GeV measured by the STAR experiment at the Relativistic Heavy Ion Collider (RHIC). The scaled factorial moments of identified charged hadrons are analyzed at mid-rapidity and within the transverse momentum phase space. We observe a power-law behavior of scaled factorial moments in Au$+$Au collisions and a decrease in the extracted scaling exponent ($\nu$) from peripheral to central collisions. The $\nu$ is consistent with a constant for different collisions energies in the mid-central (10-40%) collisions. Moreover, the $\nu$ in the 0-5% most central Au$+$Au collisions exhibits a non-monotonic energy dependence that reaches a possible minimum around $\sqrt{s_\mathrm{_{NN}}}$ = 27 GeV. The physics implications on the QCD phase structure are discussed.

We present the directed flow ($v_1$) measured in Au+Au collisions at $\sqrt{s_{_{NN}}}$ = 62.4 GeV in the mid-pseudorapidity region $|\eta|<1.3$ and in the forward pseudorapidity region $2.5 < |\eta| < 4.0$. The results are obtained using the three-particle cumulant method, the event plane method with mixed harmonics, and for the first time at the Relativistic Heavy Ion Collider (RHIC), the standard method with the event plane reconstructed from spectator neutrons. Results from all three methods are in good agreement. Over the pseudorapidity range studied, charged particle directed flow is in the direction opposite to that of fragmentation neutrons.

Vector mesons may be photoproduced in relativistic heavy-ion collisions when a virtual photon emitted by one nucleus scatters from the other nucleus, emerging as a vector meson. The STAR Collaboration has previously presented measurements of coherent $\rho^0$ photoproduction at center of mass energies of 130 GeV and 200 GeV in AuAu collisions. Here, we present a measurement of the cross section at 62.4 GeV; we find that the cross section for coherent $\rho^0$ photoproduction with nuclear breakup is $10.5\pm1.5\pm 1.6$ mb at 62.4 GeV. The cross-section ratio between 200 GeV and 62.4 GeV is $2.8\pm0.6$, less than is predicted by most theoretical models. It is, however, proportionally much larger than the previously observed $15\pm 55$% increase between 130 GeV and 200 GeV.

We present two-particle $p_{\rm t}$ correlations as a function of event centrality for Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 7.7, 11.5, 14.5, 19.6, 27, 39, 62.4, and 200 GeV at the Relativistic Heavy Ion Collider using the STAR detector. These results are compared to previous measurements from CERES at the Super Proton Synchrotron and from ALICE at the Large Hadron Collider. The data are compared with UrQMD model calculations and with a model based on a Boltzmann-Langevin approach incorporating effects from thermalization. The relative dynamical correlations for Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 200 GeV show a power law dependence on the number of participant nucleons and agree with the results for Pb+Pb collisions at $\sqrt{s_{\rm NN}} = 2.76~ {\rm TeV}$ from ALICE. As the collision energy is lowered from $\sqrt{s_{\rm NN}}$ = 200 GeV to 7.7 GeV, the centrality dependence of the relative dynamical correlations departs from the power law behavior observed at the higher collision energies. In central collisions, the relative dynamical correlations increase with collision energy up to $\sqrt{s_{\rm NN}}$ = 200 GeV in contrast to previous measurements that showed little dependence on the collision energy.

The two-particle angular correlation functions, $R_2$, of pions, kaons, and protons in Au+Au collisions at $\sqrt{s_{NN}}=$ 7.7, 11.5, 14.5, 19.6, 27, 39, 62.4, and 200 GeV were measured by the STAR experiment at RHIC. These correlations were measured for both like-sign and unlike-sign charge combinations and versus the centrality. The correlations of pions and kaons show the expected near-side ({\it i.e.}, at small relative angles) peak resulting from short-range mechanisms. The amplitudes of these short-range correlations decrease with increasing beam energy. However, the proton correlation functions exhibit strong anticorrelations in the near-side region. This behavior is observed for the first time in an A+A collision system. The observed anticorrelation is $p_{T}$-independent and decreases with increasing beam energy and centrality. The experimental results are also compared to the Monte Carlo models UrQMD, Hijing, and AMPT.

Notwithstanding decades of progress since Yukawa first developed a description of the force between nucleons in terms of meson exchange, a full understanding of the strong interaction remains a major challenge in modern science. One remaining difficulty arises from the non-perturbative nature of the strong force, which leads to the phenomenon of quark confinement at distances on the order of the size of the proton. Here we show that in relativistic heavy-ion collisions, where quarks and gluons are set free over an extended volume, two species of produced vector (spin-1) mesons, namely $\phi$ and $K^{*0}$, emerge with a surprising pattern of global spin alignment. In particular, the global spin alignment for $\phi$ is unexpectedly large, while that for $K^{*0}$ is consistent with zero. The observed spin-alignment pattern and magnitude for the $\phi$ cannot be explained by conventional mechanisms, while a model with a connection to strong force fields, i.e. an effective proxy description within the Standard Model and Quantum Chromodynamics, accommodates the current data. This connection, if fully established, will open a potential new avenue for studying the behaviour of strong force fields.

We present results from a harmonic decomposition of two-particle azimuthal correlations measured with the STAR detector in Au+Au collisions for energies ranging from $\sqrt{s_{NN}}=7.7$ GeV to 200 GeV. The third harmonic $v_3^2\{2\}=\langle \cos3(\phi_1-\phi_2)\rangle$, where $\phi_1-\phi_2$ is the angular difference in azimuth, is studied as a function of the pseudorapidity difference between particle pairs $\Delta\eta = \eta_1-\eta_2$. Non-zero {\vthree} is directly related to the previously observed large-$\Delta\eta$ narrow-$\Delta\phi$ ridge correlations and has been shown in models to be sensitive to the existence of a low viscosity Quark Gluon Plasma (QGP) phase. For sufficiently central collisions, $v_3^2\{2\}$ persist down to an energy of 7.7 GeV suggesting that QGP may be created even in these low energy collisions. In peripheral collisions at these low energies however, $v_3^2\{2\}$ is consistent with zero. When scaled by pseudorapidity density of charged particle multiplicity per participating nucleon pair, $v_3^2\{2\}$ for central collisions shows a minimum near {\snn}$=20$ GeV.

Local parity-odd domains are theorized to form inside a Quark-Gluon-Plasma (QGP) which has been produced in high-energy heavy-ion collisions. The local parity-odd domains manifest themselves as charge separation along the magnetic field axis via the chiral magnetic effect (CME). The experimental observation of charge separation has previously been reported for heavy-ion collisions at the top RHIC energies. In this paper, we present the results of the beam-energy dependence of the charge correlations in Au+Au collisions at midrapidity for center-of-mass energies of 7.7, 11.5, 19.6, 27, 39 and 62.4 GeV from the STAR experiment. After background subtraction, the signal gradually reduces with decreased beam energy, and tends to vanish by 7.7 GeV. The implications of these results for the CME will be discussed.

We study the beam-energy and system-size dependence of \phi meson production (using the hadronic decay mode \phi -- K+K-) by comparing the new results from Cu+Cu collisions and previously reported Au+Au collisions at \sqrt{s_NN} = 62.4 and 200 GeV measured in the STAR experiment at RHIC. Data presented are from mid-rapidity (|y|&lt;0.5) for 0.4 &lt; pT &lt; 5 GeV/c. At a given beam energy, the transverse momentum distributions for \phi mesons are observed to be similar in yield and shape for Cu+Cu and Au+Au colliding systems with similar average numbers of participating nucleons. The \phi meson yields in nucleus-nucleus collisions, normalised by the average number of participating nucleons, are found to be enhanced relative to those from p+p collisions with a different trend compared to strange baryons. The enhancement for \phi mesons is observed to be higher at \sqrt{s_NN} = 200 GeV compared to 62.4 GeV. These observations for the produced \phi(s\bar{s}) mesons clearly suggest that, at these collision energies, the source of enhancement of strange hadrons is related to the formation of a dense partonic medium in high energy nucleus-nucleus collisions and cannot be alone due to canonical suppression of their production in smaller systems.

Measurements of the elliptic flow, $v_{2}$, of identified hadrons ($\pi^{\pm}$, $K^{\pm}$, $K_{s}^{0}$, $p$, $\bar{p}$, $\phi$, $\Lambda$, $\bar{\Lambda}$, $\Xi^{-}$, $\bar{\Xi}^{+}$, $\Omega^{-}$, $\bar{\Omega}^{+}$) in Au+Au collisions at $\sqrt{s_{NN}}=$ 7.7, 11.5, 19.6, 27, 39 and 62.4 GeV are presented. The measurements were done at mid-rapidity using the Time Projection Chamber and the Time-of-Flight detectors of the STAR experiment during the Beam Energy Scan program at RHIC. A significant difference in the $v_{2}$ values for particles and the corresponding anti-particles was observed at all transverse momenta for the first time. The difference increases with decreasing center-of-mass energy, $\sqrt{s_{NN}}$ (or increasing baryon chemical potential, $\mu_{B}$) and is larger for the baryons as compared to the mesons. This implies that particles and anti-particles are no longer consistent with the universal number-of-constituent quark (NCQ) scaling of $v_{2}$ that was observed at $\sqrt{s_{NN}}=$ 200 GeV. However, for the group of particles NCQ scaling at $(m_{T}-m_{0})/n_{q}>$ 0.4 GeV/$c^{2}$ is not violated within $\pm$10%. The $v_{2}$ values for $\phi$ mesons at 7.7 and 11.5 GeV are approximately two standard deviations from the trend defined by the other hadrons at the highest measured $p_{T}$ values.

A systematic study is presented for centrality, transverse momentum ($p_T$) and pseudorapidity ($\eta$) dependence of the inclusive charged hadron elliptic flow ($v_2$) at midrapidity($|\eta| < 1.0$) in Au+Au collisions at $\sqrt{s_{NN}}$ = 7.7, 11.5, 19.6, 27 and 39 GeV. The results obtained with different methods, including correlations with the event plane reconstructed in a region separated by a large pseudorapidity gap and 4-particle cumulants ($v_2{4}$), are presented in order to investigate non-flow correlations and $v_2$ fluctuations. We observe that the difference between $v_2{2}$ and $v_2{4}$ is smaller at the lower collision energies. Values of $v_2$, scaled by the initial coordinate space eccentricity, $v_{2}/\varepsilon$, as a function of $p_T$ are larger in more central collisions, suggesting stronger collective flow develops in more central collisions, similar to the results at higher collision energies. These results are compared to measurements at higher energies at the Relativistic Heavy Ion Collider ($\sqrt{s_{NN}}$ = 62.4 and 200 GeV) and at the Large Hadron Collider (Pb + Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV). The $v_2(p_T)$ values for fixed $p_T$ rise with increasing collision energy within the $p_T$ range studied ($< 2 {\rm GeV}/c$). A comparison to viscous hydrodynamic simulations is made to potentially help understand the energy dependence of $v_{2}(p_{T})$. We also compare the $v_2$ results to UrQMD and AMPT transport model calculations, and physics implications on the dominance of partonic versus hadronic phases in the system created at Beam Energy Scan (BES) energies are discussed.

Identified charged particle spectra of $\pi^{\pm}$, $K^{\pm}$, $p$ and $\pbar$ at mid-rapidity ($|y|<0.1$) measured by the $\dedx$ method in the STAR-TPC are reported for $pp$ and d+Au collisions at $\snn = 200$ GeV and for Au+Au collisions at 62.4 GeV, 130 GeV, and 200 GeV. ... [Shortened for arXiv list. Full abstract in manuscript.]

We present results on strange and multi-strange particle production in Au+Au collisions at $\sqrt{s_{NN}}=62.4$ GeV as measured with the STAR detector at RHIC. Mid-rapidity transverse momentum spectra and integrated yields of $K^{0}_{S}$, $\Lambda$, $\Xi$, $\Omega$ and their anti-particles are presented for different centrality classes. The particle yields and ratios follow a smooth energy dependence. Chemical freeze-out parameters, temperature, baryon chemical potential and strangeness saturation factor obtained from the particle yields are presented. Intermediate transverse momentum ($p_T$) phenomena are discussed based on the ratio of the measured baryon-to-meson spectra and nuclear modification factor. The centrality dependence of various measurements presented show a similar behavior as seen in Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

Elliptic flow (v_2) values for identified particles at midrapidity in Au + Au collisions measured by the STAR experiment in the Beam Energy Scan at the Relativistic Heavy Ion Collider at sqrt{s_{NN}}= 7.7--62.4 GeV are presented for three centrality classes. The centrality dependence and the data at sqrt{s_{NN}}= 14.5 GeV are new. Except at the lowest beam energies we observe a similar relative v_2 baryon-meson splitting for all centrality classes which is in agreement within 15% with the number-of-constituent quark scaling. The larger v_2 for most particles relative to antiparticles, already observed for minimum bias collisions, shows a clear centrality dependence, with the largest difference for the most central collisions. Also, the results are compared with A Multiphase Transport Model and fit with a Blast Wave model.

We present the first measurement of multiplicity and pseudorapidity distributions of photons in the pseudorapidity region 2.3 $\le$ $\eta$ $\le$ 3.7 for different centralities in Au + Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV. We find that the photon yield in this pseudorapidity range scales with the number of participating nucleons at all collision centralities studied. The pseudorapidity distribution of photons, dominated by neutral pion decays, has been compared to those of identified charged pions, photons, and inclusive charged particles from heavy ion and nucleon-nucleon collisions at various energies. The photon production in the measured pseudorapidity region has been shown to be consistent with the energy and centrality independent limiting fragmentation scenario.

We present the results of an elliptic flow analysis of Cu+Cu collisions recorded with the STAR detector at 62.4 and 200GeV. Elliptic flow as a function of transverse momentum is reported for different collision centralities for charged hadrons and strangeness containing hadrons $K_{S}^{0}$, $\Lambda$, $\Xi$, $\phi$ in the midrapidity region $|eta|<1.0$. Significant reduction in systematic uncertainty of the measurement due to non-flow effects has been achieved by correlating particles at midrapidity, $|\eta|<1.0$, with those at forward rapidity, $2.5<|\eta|<4.0$. We also present azimuthal correlations in p+p collisions at 200 GeV to help estimating non-flow effects. To study the system-size dependence of elliptic flow, we present a detailed comparison with previously published results from Au+Au collisions at 200 GeV. We observe that $v_{2}$($p_{T}$) of strange hadrons has similar scaling properties as were first observed in Au+Au collisions, i.e.: (i) at low transverse momenta, $p_T<2GeV/c$, $v_{2}$ scales with transverse kinetic energy, $m_{T}-m$, and (ii) at intermediate $p_T$, $2<p_T<4GeV/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 $v_{2}/\epsilon$ depend on the system size, 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.

We report a systematic measurement of cumulants, $C_{n}$, for net-proton, proton and antiproton multiplicity distributions, and correlation functions, $\kappa_n$, for proton and antiproton multiplicity distributions up to the fourth order in Au+Au collisions at $\sqrt{s_{\mathrm {NN}}}$ = 7.7, 11.5, 14.5, 19.6, 27, 39, 54.4, 62.4 and 200 GeV. The $C_{n}$ and $\kappa_n$ are presented as a function of collision energy, centrality and kinematic acceptance in rapidity, $y$, and transverse momentum, $p_{T}$. The data were taken during the first phase of the Beam Energy Scan (BES) program (2010 -- 2017) at the BNL Relativistic Heavy Ion Collider (RHIC) facility. The measurements are carried out at midrapidity ($|y| <$ 0.5) and transverse momentum 0.4 $<$$p_{\rm T}$$<$ 2.0 GeV/$c$, using the STAR detector at RHIC. We observe a non-monotonic energy dependence ($\sqrt{s_{\mathrm {NN}}}$ = 7.7 -- 62.4 GeV) of the net-proton $C_{4}$/$C_{2}$ with the significance of 3.1$\sigma$ for the 0-5% central Au+Au collisions. This is consistent with the expectations of critical fluctuations in a QCD-inspired model. Thermal and transport model calculations show a monotonic variation with $\sqrt{s_{\mathrm {NN}}}$. For the multiparticle correlation functions, we observe significant negative values for a two-particle correlation function, $\kappa_2$, of protons and antiprotons, which are mainly due to the effects of baryon number conservation. Furthermore, it is found that the four-particle correlation function, $\kappa_4$, of protons plays a role in determining the energy dependence of proton $C_4/C_1$ below 19.6 GeV, which cannot be understood by the effect of baryon number conservation.