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 report a new measurement of the production of electrons from open heavy-flavor hadron decays (HFEs) at mid-rapidity ($|y|<$ 0.7) in Au+Au collisions at $\sqrt{s_{\rm NN}}=200$ GeV. Invariant yields of HFEs are measured for the transverse momentum range of $3.5 < p_{\rm T} < 9$ GeV/$c$ in various configurations of the collision geometry. The HFE yields in head-on Au+Au collisions are suppressed by approximately a factor of 2 compared to that in $p$+$p$ collisions scaled by the average number of binary collisions, indicating strong interactions between heavy quarks and the hot and dense medium created in heavy-ion collisions. Comparison of these results with models provides additional tests of theoretical calculations of heavy quark energy loss in the quark-gluon plasma.

We report a new measurement of $D^0$-meson production at mid-rapidity ($|y|$\,$<$\,1) in Au+Au collisions at ${\sqrt{s_{\rm NN}} = \rm{200\,GeV}}$ utilizing the Heavy Flavor Tracker, a high resolution silicon detector at the STAR experiment. Invariant yields of $D^0$-mesons with transverse momentum $p_{T}$ $\lesssim 9$\,GeV/$c$ are reported in various centrality bins (0--10\%, 10--20\%, 20--40\%, 40--60\% and 60--80\%). Blast-Wave thermal models are used to fit the $D^0$-meson $p_{T}$ spectra to study $D^0$ hadron kinetic freeze-out properties. The average radial flow velocity extracted from the fit is considerably smaller than that of light hadrons ($\pi,K$ and $p$), but comparable to that of hadrons containing multiple strange quarks ($\phi,\Xi^-$), indicating that $D^0$ mesons kinetically decouple from the system earlier than light hadrons. The calculated $D^0$ nuclear modification factors re-affirm that charm quarks suffer large amount of energy loss in the medium, similar to those of light quarks for $p_{T}$\,$>$\,4\,GeV/$c$ in central 0--10\% Au+Au collisions. At low $p_{T}$, the nuclear modification factors show a characteristic structure qualitatively consistent with the expectation from model predictions that charm quarks gain sizable collective motion during the medium evolution. The improved measurements are expected to offer new constraints to model calculations and help gain further insights into the hot and dense medium created in these collisions.

A study of multiplicity and pseudorapidity distributions of inclusive photons measured in pp and p$-$Pb collisions at a center-of-mass energy per nucleon$-$nucleon collision of $\sqrt{s_{\rm NN}} = 5.02$ TeV using the ALICE detector in the forward pseudorapidity region $2.3 < \eta_{\rm lab} < 3.9$ is presented. Measurements in p$-$Pb collisions are reported for two beam configurations in which the directions of the proton and lead ion beam were reversed. The pseudorapidity distributions in p$-$Pb collisions are obtained for seven centrality classes which are defined based on different event activity estimators, i.e., the charged-particle multiplicity measured at midrapidity as well as the energy deposited in a calorimeter at beam rapidity. The inclusive photon multiplicity distributions for both pp and p$-$Pb collisions are described by double negative binomial distributions. The pseudorapidity distributions of inclusive photons are compared to those of charged particles at midrapidity in \pp collisions and for different centrality classes in p$-$Pb collisions. The results are compared to predictions from various Monte Carlo event generators. None of the generators considered in this paper reproduces the inclusive photon multiplicity distributions in the reported multiplicity range. The pseudorapidity distributions are, however, better described by the same generators.

We report the triton ($t$) production in mid-rapidity ($|y| <$ 0.5) Au+Au collisions at $\sqrt{s_\mathrm{NN}}$= 7.7--200 GeV measured by the STAR experiment from the first phase of the beam energy scan at the Relativistic Heavy Ion Collider (RHIC). The nuclear compound yield ratio ($\mathrm{N}_t \times \mathrm{N}_p/\mathrm{N}_d^2$), which is predicted to be sensitive to the fluctuation of local neutron density, is observed to decrease monotonically with increasing charged-particle multiplicity ($dN_{ch}/d\eta$) and follows a scaling behavior. The $dN_{ch}/d\eta$ dependence of the yield ratio is compared to calculations from coalescence and thermal models. Enhancements in the yield ratios relative to the coalescence baseline are observed in the 0%-10% most central collisions at 19.6 and 27 GeV, with a significance of 2.3$\sigma$ and 3.4$\sigma$, respectively, giving a combined significance of 4.1$\sigma$. The enhancements are not observed in peripheral collisions or model calculations without critical fluctuation, and decreases with a smaller $p_{T}$ acceptance. The physics implications of these results on the QCD phase structure and the production mechanism of light nuclei in heavy-ion collisions are discussed.

Presented are the first measurements of the transverse single-spin asymmetries ($A_N$) for neutral pions and eta mesons in $p$+Au and $p$+Al collisions at $\sqrt{s_{_{NN}}}=200$ GeV in the pseudorapidity range $|\eta|<$0.35 with the PHENIX detector at the Relativistic Heavy Ion Collider. The asymmetries are consistent with zero, similar to those for midrapidity neutral pions and eta mesons produced in $p$+$p$ collisions. These measurements show no evidence of additional effects that could potentially arise from the more complex partonic environment present in proton-nucleus collisions.

Recently, the PHENIX Collaboration has published second- and third-harmonic Fourier coefficients $v_2$ and $v_3$ for midrapidity ($|\eta|<0.35$) charged hadrons in 0%--5% central $p$$+$Au, $d$ $+$Au, and $^3$He$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV utilizing three sets of two-particle correlations for two detector combinations with different pseudorapidity acceptance [Phys. Rev. C {\bf 105}, 024901 (2022)]. This paper extends these measurements of $v_2$ to all centralities in $p$ $+$Au, $d$ $+$Au, and $^3$He$+$Au collisions, as well as $p$$+$$p$ collisions, as a function of transverse momentum ($p_T$) and event multiplicity. The kinematic dependence of $v_2$ is quantified as the ratio $R$ of $v_2$ between the two detector combinations as a function of event multiplicity for $0.5$ $<$ $p_T$ $<$ $1$ and $2$ $<$ $p_T$ $<$ $2.5$ GeV/$c$. A multiphase-transport (AMPT) model can reproduce the observed $v_2$ in most-central to midcentral $d$$+$Au and $^3$He$+$Au collisions. However, the AMPT model systematically overestimates the measurements in $p$ $+$ $p$, $p$ $+$Au, and peripheral $d$$+$Au and $^3$He$+$Au collisions, indicating a higher nonflow contribution in AMPT than in the experimental data. The AMPT model fails to describe the observed $R$ for $0.5$ $<$ $p_T$$<$ $1$ GeV/$c$, but there is qualitative agreement with the measurements for $2$ $<$ $p_T$ $<$ $2.5$ GeV/$c$.

The STAR Collaboration reports measurements of the transverse single-spin asymmetries, $A_N$, for inclusive jets and identified `hadrons within jets' production at midrapidity from transversely polarized $pp$ collisions at $\sqrt{s}$ = 200 GeV, based on data recorded in 2012 and 2015. The inclusive jet asymmetry measurements include $A_N$ for inclusive jets and $A_N$ for jets containing a charged pion carrying a momentum fraction $z>0.3$ of the jet momentum. The identified hadron within jet asymmetry measurements include the Collins effect for charged pions, kaons and protons, and the Collins-like effect for charged pions. The measured asymmetries are determined for several distinct kinematic regions, characterized by the jet transverse momentum $p_{T}$ and pseudorapidity $\eta$, as well as the hadron momentum fraction $z$ and momentum transverse to the jet axis $j_{T}$. These results probe higher momentum scales ($Q^{2}$ up to $\sim$ 900 GeV$^{2}$) than current, semi-inclusive deep inelastic scattering measurements, and they provide new constraints on quark transversity in the proton and enable tests of evolution, universality and factorization breaking in the transverse-momentum-dependent formalism.

We present the first inclusive measurements of the invariant and SoftDrop jet mass in proton-proton collisions at $\sqrt{s}=200$ GeV at STAR. The measurements are fully corrected for detector effects, and reported differentially in both the jet transverse momentum and jet radius parameter. We compare the measurements to established leading-order Monte Carlo event generators and find that STAR-tuned PYTHIA-6 reproduces the data, while LHC tunes of PYTHIA-8 and HERWIG-7 do not agree with the data, providing further constraints on parameter tuning. Finally, we observe that SoftDrop grooming, for which the contribution of wide-angle non-perturbative radiation is suppressed, shifts the jet mass distributions into closer agreement with the partonic jet mass as determined by both PYTHIA-8 and a next-to-leading-logarithmic accuracy perturbative QCD calculation. These measurements complement recent LHC measurements in a different kinematic region, as well as establish a baseline for future jet mass measurements in heavy-ion collisions at RHIC.

Two-particle correlation measurements projected onto two-dimensional, transverse rapidity coordinates ($y_{T1},y_{T2}$), allow access to dynamical properties of the QCD medium produced in relativistic heavy-ion collisions that angular correlation measurements are not sensitive to. We report non-identified charged-particle correlations for Au + Au minimum-bias collisions at $\sqrt{s_{\rm NN}}$ = 200 GeV taken by the STAR experiment at the Relativistic Heavy-Ion Collider (RHIC). Correlations are presented as 2D functions of transverse rapidity for like-sign, unlike-sign and all charged-particle pairs, as well as for particle pairs whose relative azimuthal angles lie on the near-side, the away-side, or at all relative azimuth. The correlations are constructed using charged particles with transverse momentum $p_T \geq 0.15$ GeV/$c$, pseudorapidity from $-$1 to 1, and azimuthal angles from $-\pi$ to $\pi$. The significant correlation structures that are observed evolve smoothly with collision centrality. The major correlation features include a saddle shape plus a broad peak with maximum near $y_T \approx 3$, corresponding to $p_T \approx$ 1.5 GeV/$c$. The broad peak is observed in both like- and unlike-sign charge combinations and in near- and away-side relative azimuthal angles. The all-charge, all-azimuth correlation measurements are compared with the theoretical predictions of {\sc hijing} and {\sc epos}. The results indicate that the correlations for peripheral to mid-central collisions can be approximately described as a superposition of nucleon + nucleon collisions with minimal effects from the QCD medium. Strong medium effects are indicated in mid- to most-central collisions.

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.

We present high-precision measurements of elliptic, triangular, and quadrangular flow $v_{2}$, $v_{3}$, and $v_{4}$, respectively, at midrapidity ($|\eta|<1.0$) for identified hadrons $\pi$, $p$, $K$, $\varphi$, $K_s$, $\Lambda$ as a function of centrality and transverse momentum in Au+Au collisions at the center-of-mass energy $\sqrt{s_{\rm NN}}=$ 200 GeV. We observe similar $v_{n}$ trends between light and strange mesons which indicates that the heavier strange quarks flow as strongly as the lighter up and down quarks. The number-of-constituent-quark scaling for $v_{2}$, $v_{3}$, and $v_{4}$ is found to hold within statistical uncertainty for 0-10$\%$, 10-40$\%$ and 40-80$\%$ collision centrality intervals. The results are compared to several viscous hydrodynamic calculations with varying initial conditions, and could serve as an additional constraint to the development of hydrodynamic models.

We present the first measurements of the longitudinal double-spin asymmetry $A_{LL}$ for dijets with at least one jet reconstructed within the pseudorapidity range $0.8 < \eta < 1.8$. The dijets were measured in polarized $pp$ collisions at a center-of-mass energy $\sqrt{s}$ = 200 GeV. Values for $A_{LL}$ are determined for several distinct event topologies, defined by the jet pseudorapidities, and span a range of parton momentum fraction $x$ down to $x \sim$ 0.01. The measured asymmetries are found to be consistent with the predictions of global analyses that incorporate the results of previous RHIC measurements. They will provide new constraints on $\Delta g(x)$ in this poorly constrained region when included in future global analyses.

Polarized proton-proton collisions provide leading-order access to gluons, presenting an opportunity to constrain gluon spin-momentum correlations within transversely polarized protons and enhance our understanding of the three-dimensional structure of the proton. Midrapidity open-heavy-flavor production at $\sqrt{s}=200$ GeV is dominated by gluon-gluon fusion, providing heightened sensitivity to gluon dynamics relative to other production channels. Transverse single-spin asymmetries of positrons and electrons from heavy-flavor hadron decays are measured at midrapidity using the PHENIX detector at the Relativistic Heavy Ion Collider. These charge-separated measurements are sensitive to gluon correlators that can in principle be related to gluon orbital angular momentum via model calculations. Explicit constraints on gluon correlators are extracted for two separate models, one of which had not been constrained previously.

We present measurements of the cross section and double-helicity asymmetry $A_{LL}$ of direct-photon production in $\vec{p}+\vec{p}$ collisions at $\sqrt{s}=510$ GeV. The measurements have been performed at midrapidity ($|\eta|<0.25$) with the PHENIX detector at the Relativistic Heavy Ion Collider. At relativistic energies, direct photons are dominantly produced from the initial quark-gluon hard scattering and do not interact via the strong force at leading order. Therefore, at $\sqrt{s}=510$ GeV, where leading-order-effects dominate, these measurements provide clean and direct access to the gluon helicity in the polarized proton in the gluon-momentum-fraction range $0.02<x<0.08$, with direct sensitivity to the sign of the gluon contribution.

Rapidity-odd directed flow measurements at midrapidity are presented for $\Lambda$, $\bar{\Lambda}$, $K^\pm$, $K^0_s$ and $\phi$ at $\sqrt{s_{NN}} =$ 7.7, 11.5, 14.5, 19.6, 27, 39, 62.4 and 200 GeV in Au+Au collisions recorded by the STAR detector at the Relativistic Heavy Ion Collider. These measurements greatly expand the scope of data available to constrain models with differing prescriptions for the equation of state of quantum chromodynamics. Results show good sensitivity for testing a picture where flow is assumed to be imposed before hadron formation and the observed particles are assumed to form via coalescence of constituent quarks. The pattern of departure from a coalescence-inspired sum-rule can be a valuable new tool for probing the collision dynamics.

The chiral magnetic effect (CME) is predicted to occur as a consequence of a local violation of $\cal P$ and $\cal CP$ symmetries of the strong interaction amidst a strong electro-magnetic field generated in relativistic heavy-ion collisions. Experimental manifestation of the CME involves a separation of positively and negatively charged hadrons along the direction of the magnetic field. Previous measurements of the CME-sensitive charge-separation observables remain inconclusive because of large background contributions. In order to better control the influence of signal and backgrounds, the STAR Collaboration performed a blind analysis of a large data sample of approximately 3.8 billion isobar collisions of $^{96}_{44}$Ru+$^{96}_{44}$Ru and $^{96}_{40}$Zr+$^{96}_{40}$Zr at $\sqrt{s_{\rm NN}}=200$ GeV. Prior to the blind analysis, the CME signatures are predefined as a significant excess of the CME-sensitive observables in Ru+Ru collisions over those in Zr+Zr collisions, owing to a larger magnetic field in the former. A precision down to 0.4% is achieved, as anticipated, in the relative magnitudes of the pertinent observables between the two isobar systems. Observed differences in the multiplicity and flow harmonics at the matching centrality indicate that the magnitude of the CME background is different between the two species. No CME signature that satisfies the predefined criteria has been observed in isobar collisions in this blind analysis.

Understanding gluon density distributions and how they are modified in nuclei are among the most important goals in nuclear physics. In recent years, diffractive vector meson production measured in ultra-peripheral collisions (UPCs) at heavy-ion colliders has provided a new tool for probing the gluon density. In this Letter, we report the first measurement of $J/\psi$ photoproduction off the deuteron in UPCs at the center-of-mass energy $\sqrt{s_{_{\rm NN}}}=200~\rm GeV$ in d$+$Au collisions. The differential cross section as a function of momentum transfer $-t$ is measured. In addition, data with a neutron tagged in the deuteron-going Zero-Degree Calorimeter is investigated for the first time, which is found to be consistent with the expectation of incoherent diffractive scattering at low momentum transfer. Theoretical predictions based on the Color Glass Condensate saturation model and the gluon shadowing model are compared with the data quantitatively. A better agreement with the saturation model has been observed. With the current measurement, the results are found to be directly sensitive to the gluon density distribution of the deuteron and the deuteron breakup, which provides insights into the nuclear gluonic structure.

The ALICE collaboration at the LHC reports measurement of the inclusive production cross section of electrons from semi-leptonic decays of beauty hadrons with rapidity $|y|<0.8$ and transverse momentum $1<p_{\mathrm{T}}<10$ GeV/$c$, in pp collisions at $\sqrt{s} = $ 2.76 TeV. Electrons not originating from semi-electronic decay of beauty hadrons are suppressed using the impact parameter of the corresponding tracks. The production cross section of beauty decay electrons is compared to the result obtained with an alternative method which uses the distribution of the azimuthal angle between heavy-flavour decay electrons and charged hadrons. Perturbative QCD calculations agree with the measured cross section within the experimental and theoretical uncertainties. The integrated visible cross section, $\sigma_{\mathrm{b} \rightarrow \mathrm{e}} = 3.47\pm0.40(\mathrm{stat})^{+1.12}_{-1.33}(\mathrm{sys})\pm0.07(\mathrm{norm}) \mu$b, was extrapolated to full phase space using Fixed Order plus Next-to-Leading Log (FONLL) predictions to obtain the total b$\bar{\mathrm{b}}$ production cross section, $\sigma_{\mathrm{b\bar{b}}} = 130\pm15.1(\mathrm{stat})^{+42.1}_{-49.8}(\mathrm{sys})^{+3.4}_{-3.1}(\mathrm{extr})\pm2.5(\mathrm{norm})\pm4.4(\mathrm{BR}) \mu$b.

We report a new measurement of the production cross section for inclusive electrons from open heavy-flavor hadron decays as a function of transverse momentum ($p_{\rm T}$) at mid-rapidity ($|y|<$ 0.7) in $p$+$p$ collisions at $\sqrt{s} = 200$ GeV. The result is presented for 2.5 $<p_{\rm T}<$ 10 GeV/$c$ with an improved precision above 6 GeV/$c$ with respect to the previous measurements, providing more constraints on perturbative QCD calculations. Moreover, this measurement also provides a high-precision reference for measurements of nuclear modification factors for inclusive electrons from open-charm and -bottom hadron decays in heavy-ion collisions.

We report the first measurements of a complete second-order cumulant matrix of net-charge, net-proton, and net-kaon multiplicity distributions for the first phase of the beam energy scan program at RHIC. This includes the centrality and, for the first time, the pseudorapidity window dependence of both diagonal and off-diagonal cumulants in Au+Au collisions at \sNN~= 7.7-200 GeV. Within the available acceptance of $|\eta|<0.5$, the cumulants grow linearly with the pseudorapidity window. Relative to the corresponding measurements in peripheral collisions, the ratio of off-diagonal over diagonal cumulants in central collisions indicates an excess correlation between net-charge and net-kaon, as well as between net-charge and net-proton. The strength of such excess correlation increases with the collision energy. The correlation between net-proton and net-kaon multiplicity distributions is observed to be negative at \sNN~= 200 GeV and change to positive at the lowest collision energy. Model calculations based on non-thermal (UrQMD) and thermal (HRG) production of hadrons cannot explain the data. These measurements will help map the QCD phase diagram, constrain hadron resonance gas model calculations, and provide new insights on the energy dependence of baryon-strangeness correlations. An erratum has been added to address the issue of self-correlation in the previously considered efficiency correction for off-diagonal cumulant measurement. Previously considered unidentified (net-)charge correlation results ($\sigma^{11}_{Q,p}$ and $\sigma^{11}_{Q,k})$ are now replaced with identified (net-)charge correlation ($\sigma^{11}_{Q^{PID},p}$ and $\sigma^{11}_{Q^{PID},k}$)

According to first-principle lattice QCD calculations, the transition from quark-gluon plasma to hadronic matter is a smooth crossover in the region μB ≤ T c. In this range the ratio, C6=C2, of net-baryon distributions are predicted to be negative. In this Letter, we report the first measurement of the midrapidity net-proton C6=C2 from 27, 54.4, and 200 GeV Au þ Au collisions at the Relativistic Heavy Ion Collider (RHIC). The dependence on collision centrality and kinematic acceptance in (p T , y) are analyzed. While for 27 and 54.4 GeV collisions the C6=C2 values are close to zero within uncertainties, it is observed that for 200 GeV collisions, the C6=C2 ratio becomes progressively negative from peripheral to central collisions. Transport model calculations without critical dynamics predict mostly positive values except for the most central collisions within uncertainties. These observations seem to favor a smooth crossover in the high-energy nuclear collisions at top RHIC 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.

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.

The STAR Collaboration reports measurements of the transverse single-spin asymmetry (TSSA) of inclusive $\pi^0$ at center-of-mass energies ($\sqrt s$) of 200 GeV and 500 GeV in transversely polarized proton-proton collisions in the pseudo-rapidity region 2.7 to 4.0. The results at the two different energies show a continuous increase of the TSSA with Feynman-$x$, and, when compared to previous measurements, no dependence on $\sqrt s$ from 19.4 GeV to 500 GeV is found. To investigate the underlying physics leading to this large TSSA, different topologies have been studied. $\pi^0$ with no nearby particles tend to have a higher TSSA than inclusive $\pi^0$. The TSSA for inclusive electromagnetic jets, sensitive to the Sivers effect in the initial state, is substantially smaller, but shows the same behavior as the inclusive $\pi^0$ asymmetry as a function of Feynman-$x$. To investigate final-state effects, the Collins asymmetry of $\pi^0$ inside electromagnetic jets has been measured. The Collins asymmetry is analyzed for its dependence on the $\pi^0$ momentum transverse to the jet thrust axis and its dependence on the fraction of jet energy carried by the $\pi^0$. The asymmetry was found to be small in each case for both center-of-mass energies. All the measurements are compared to QCD-based theoretical calculations for transverse-momentum-dependent parton distribution functions and fragmentation functions. Some discrepancies are found, which indicates new mechanisms might be involved.