The linear and mode-coupled contributions to higher-order anisotropic flow are presented for Au+Au collisions at $\sqrt{s_{\mathrm{NN}}}$ = 27, 39, 54.4, and 200 GeV and compared to similar measurements for Pb+Pb collisions at the Large Hadron Collider (LHC). The coefficients and the flow harmonics' correlations, which characterize the linear and mode-coupled response to the lower-order anisotropies, indicate a beam energy dependence consistent with an influence from the specific shear viscosity ($\eta/s$). In contrast, the dimensionless coefficients, mode-coupled response coefficients, and normalized symmetric cumulants are approximately beam-energy independent, consistent with a significant role from initial-state effects. These measurements could provide unique supplemental constraints to (i) distinguish between different initial-state models and (ii) delineate the temperature ($T$) and baryon chemical potential ($\mu_{B}$) dependence of the specific shear viscosity $\frac{\eta}{s} (T, \mu_B)$.

A linearly polarized photon can be quantized from the Lorentz-boosted electromagnetic field of a nucleus traveling at ultra-relativistic speed. When two relativistic heavy nuclei pass one another at a distance of a few nuclear radii, the photon from one nucleus may interact through a virtual quark-antiquark pair with gluons from the other nucleus forming a short-lived vector meson (e.g. ${\rho^0}$). In this experiment, the polarization was utilized in diffractive photoproduction to observe a unique spin interference pattern in the angular distribution of ${\rho^0\rightarrow\pi^+\pi^-}$ decays. The observed interference is a result of an overlap of two wave functions at a distance an order of magnitude larger than the ${\rho^0}$ travel distance within its lifetime. The strong-interaction nuclear radii were extracted from these diffractive interactions, and found to be $6.53\pm 0.06$ fm ($^{197} {\rm Au }$) and $7.29\pm 0.08$ fm ($^{238} {\rm U}$), larger than the nuclear charge radii. The observable is demonstrated to be sensitive to the nuclear geometry and quantum interference of non-identical particles.

The polarization of $\Lambda$ and $\bar{\Lambda}$ hyperons along the beam direction has been measured relative to the second and third harmonic event planes in isobar Ru+Ru and Zr+Zr collisions at $\sqrt{s_{NN}}$ = 200 GeV. This is the first experimental evidence of the hyperon polarization by the triangular flow originating from the initial density fluctuations. The amplitudes of the sine modulation for the second and third harmonic results are comparable in magnitude, increase from central to peripheral collisions, and show a mild $p_T$ dependence. The azimuthal angle dependence of the polarization follows the vorticity pattern expected due to elliptic and triangular anisotropic flow, and qualitatively disagree with most hydrodynamic model calculations based on thermal vorticity and shear induced contributions. The model results based on one of existing implementations of the shear contribution lead to a correct azimuthal angle dependence, but predict centrality and $p_T$ dependence that still disagree with experimental measurements. Thus, our results provide stringent constraints on the thermal vorticity and shear-induced contributions to hyperon polarization. Comparison to previous measurements at RHIC and the LHC for the second-order harmonic results shows little dependence on the collision system size and collision energy.

We report here the first observation of directed flow ($v_1$) of the hypernuclei $^3_{\Lambda}$H and $^4_{\Lambda}$H in mid-central Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 3 GeV at RHIC. These data are taken as part of the beam energy scan program carried out by the STAR experiment. From 165 $\times$ 10$^{6}$ events in 5%-40% centrality, about 8400 $^3_{\Lambda}$H and 5200 $^4_{\Lambda}$H candidates are reconstructed through two- and three-body decay channels. We observe that these hypernuclei exhibit significant directed flow. Comparing to that of light nuclei, it is found that the midrapidity $v_1$ slopes of $^3_{\Lambda}$H and $^4_{\Lambda}$H follow baryon number scaling, implying that the coalescence is the dominant mechanism for these hypernuclei production in such collisions.

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.

We report first measurements of $e^{+}e^{-}$ pair production in the mass region 0.4 $

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 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.

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.

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.

Global polarization of $\Lambda$ hyperons has been measured to be of the order of a few tenths of a percent in Au+Au collisions at $\sqrt{s_{_{NN}}}$ = 200 GeV, with no significant difference between $\Lambda$ and $\bar{\Lambda}$. These new results reveal the collision energy dependence of the global polarization together with the results previously observed at $\sqrt{s_{_{NN}}}$ = 7.7 -- 62.4 GeV and indicate noticeable vorticity of the medium created in non-central heavy-ion collisions at the highest RHIC collision energy. The signal is in rough quantitative agreement with the theoretical predictions from a hydrodynamic model and from the AMPT (A Multi-Phase Transport) model. The polarization is larger in more peripheral collisions, and depends weakly on the hyperon's transverse momentum and pseudorapidity $\eta^H$ within $|\eta^H|<1$. An indication of the polarization dependence on the event-by-event charge asymmetry is observed at the $2\sigma$ level, suggesting a possible contribution to the polarization from the axial current induced by the initial magnetic field.

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 results on the total and elastic cross sections in proton-proton collisions at $\sqrt{s}=200$ GeV obtained with the Roman Pot setup of the STAR experiment at the Relativistic Heavy Ion Collider (RHIC). The elastic differential cross section was measured in the squared four-momentum transfer range $0.045 \leq -t \leq 0.135$ GeV$^2$. The value of the exponential slope parameter $B$ of the elastic differential cross section $d\sigma/dt \sim e^{-Bt}$ in the measured $-t$ range was found to be $B = 14.32 \pm 0.09 (stat.)^{\scriptstyle +0.13}_{\scriptstyle -0.28} (syst.)$ GeV$^{-2}$. The total cross section $\sigma_{tot}$, obtained from extrapolation of the $d\sigma/dt$ to the optical point at $-t = 0$, is $\sigma_{tot} = 54.67 \pm 0.21 (stat.) ^{\scriptstyle +1.28}_{\scriptstyle -1.38} (syst.)$ mb. We also present the values of the elastic cross section $\sigma_{el} = 10.85 \pm 0.03 (stat.) ^{\scriptstyle +0.49}_{\scriptstyle -0.41}(syst.)$ mb, the elastic cross section integrated within the STAR $t$-range $\sigma^{det}_{el} = 4.05 \pm 0.01 (stat.) ^{\scriptstyle+0.18}_{\scriptstyle -0.17}(syst.)$ mb, and the inelastic cross section $\sigma_{inel} = 43.82 \pm 0.21 (stat.) ^{\scriptstyle +1.37}_{\scriptstyle -1.44} (syst.)$ mb. The results are compared with the world data.

We report on the $W$ and $Z/\gamma^*$ differential and total cross sections as well as the $W^+$/$W^-$ and $(W^+ + W^-)$/$(Z/\gamma^*)$ cross-section ratios measured by the STAR experiment at RHIC in $p+p$ collisions at $\sqrt{s} = 500$ GeV and $510$ GeV. The cross sections and their ratios are sensitive to quark and antiquark parton distribution functions. In particular, at leading order, the $W$ cross-section ratio is sensitive to the $\bar{d}/\bar{u}$ ratio. These measurements were taken at high $Q^2 \sim M_W^2,M_Z^2$ and can serve as input into global analyses to provide constraints on the sea quark distributions. The results presented here combine three STAR data sets from 2011, 2012, and 2013, accumulating an integrated luminosity of 350 pb$^{-1}$. We also assess the expected impact that our $W^+/W^-$ cross-section ratios will have on various quark distributions, and find sensitivity to the $\bar{u}-\bar{d}$ and $\bar{d}/\bar{u}$ distributions.

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 STAR Collaboration at the Relativistic Heavy Ion Collider reports the first measurement of inclusive jet production in peripheral and central Au+Au collisions at $\sqrt{s_{NN}}$=200 GeV. Jets are reconstructed with the anti-k$_{T}$ algorithm using charged tracks with pseudorapidity $|\eta|<1.0$ and transverse momentum $0.2<p_{T,jet}^{ch}<30$ GeV/$c$, with jet resolution parameter $R$=0.2, 0.3, and 0.4. The large background yield uncorrelated with the jet signal is observed to be dominated by statistical phase space, consistent with a previous coincidence measurement. This background is suppressed by requiring a high-transverse-momentum (high-$p_T$) leading hadron in accepted jet candidates. The bias imposed by this requirement is assessed, and the $p_T$ region in which the bias is small is identified. Inclusive charged-particle jet distributions are reported in peripheral and central Au+Au collisions for $5<p_{T,jet}^{ch}<25$ GeV/$c$ and $5<p_{T,jet}^{ch}<30$ GeV/$c$, respectively. The charged-particle jet inclusive yield is suppressed for central Au+Au collisions, compared to both the peripheral Au+Au yield from this measurement and to the $pp$ yield calculated using the PYTHIA event generator. The magnitude of the suppression is consistent with that of inclusive hadron production at high $p_T$, and that of semi-inclusive recoil jet yield when expressed in terms of energy loss due to medium-induced energy transport. Comparison of inclusive charged-particle jet yields for different values of $R$ exhibits no significant evidence for medium-induced broadening of the transverse jet profile for $R<0.4$ in central Au+Au collisions. The measured distributions are consistent with theoretical model calculations that incorporate jet quenching.

We measure triangular flow relative to the reaction plane at 3 GeV center-of-mass energy in Au+Au collisions at the BNL Relativistic Heavy Ion Collider. A significant $v_3$ signal for protons is observed, which increases for higher rapidity, higher transverse momentum, and more peripheral collisions. The triangular flow is essentially rapidity-odd with a slope at mid-rapidity, $dv_3/dy|_{(y=0)}$, opposite in sign compared to the slope for directed flow. No significant $v_3$ signal is observed for charged pions and kaons. Comparisons with models suggest that a mean field potential is required to describe these results, and that the triangular shape of the participant nucleons is the result of stopping and nuclear geometry.

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.

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.

Azimuthal anisotropy of produced particles is one of the most important observables used to access the collective properties of the expanding medium created in relativistic heavy-ion collisions. In this paper, we present second ($v_{2}$) and third ($v_{3}$) order azimuthal anisotropies of $K_{S}^{0}$, $\phi$, $\Lambda$, $\Xi$ and $\Omega$ at mid-rapidity ($|y|<$1) in Au+Au collisions at $\sqrt{s_{\text{NN}}}$ = 54.4 GeV measured by the STAR detector. The $v_{2}$ and $v_{3}$ are measured as a function of transverse momentum and centrality. Their energy dependence is also studied. $v_{3}$ is found to be more sensitive to the change in the center-of-mass energy than $v_{2}$. Scaling by constituent quark number is found to hold for $v_{2}$ within 10%. This observation could be evidence for the development of partonic collectivity in 54.4 GeV Au+Au collisions. Differences in $v_{2}$ and $v_{3}$ between baryons and anti-baryons are presented, and ratios of $v_{3}$/$v_{2}^{3/2}$ are studied and motivated by hydrodynamical calculations. The ratio of $v_{2}$ of $\phi$ mesons to that of anti-protons ($v_{2}(\phi)/v_{2}(\bar{p})$) shows centrality dependence at low transverse momentum, presumably resulting from the larger effects from hadronic interactions on anti-proton $v_{2}$.

We report the beam energy and collision centrality dependence of fifth and sixth order cumulants ($C_{5}$, $C_{6}$) and factorial cumulants ($\kappa_{5}$, $\kappa_{6}$) of net-proton and proton distributions, from $\sqrt{s_{NN}} = 3 - 200$ GeV Au+Au collisions at RHIC. The net-proton cumulant ratios generally follow the hierarchy expected from QCD thermodynamics, except for the case of collisions at $\sqrt{s_{NN}}$ = 3 GeV. $C_{6}/C_{2}$ for 0-40% centrality collisions is increasingly negative with decreasing $\sqrt{s_{NN}}$, while it is positive for the lowest $\sqrt{s_{NN}}$ studied. These observed negative signs are consistent with QCD calculations (at baryon chemical potential, $\mu_{B} \leq$ 110 MeV) that include a crossover quark-hadron transition. In addition, for $\sqrt{s_{NN}} \geq$ 11.5 GeV, the measured proton $\kappa_{n}$, within uncertainties, does not support the two-component shape of proton distributions that would be expected from a first-order phase transition. Taken in combination, the hyper-order proton number fluctuations suggest that the structure of QCD matter at high baryon density, $\mu_{B}\sim 750$ MeV ($\sqrt{s_{NN}}$ = 3 GeV) is starkly different from those at vanishing $\mu_{B}\sim 20$MeV ($\sqrt{s_{NN}}$ = 200 GeV and higher).

We report on new measurements of inclusive J/$\psi$ polarization at mid-rapidity in p+p collisions at $\sqrt{s}$ = 200 GeV by the STAR experiment at RHIC. The polarization parameters, $\lambda_\theta$, $\lambda_\phi$, and $\lambda_{\theta\phi}$, are measured as a function of transverse momentum ($p_T$) in both the Helicity and Collins-Soper (CS) reference frames within $p_T< 10$ GeV/$C$. Except for $\lambda_\theta$ in the CS frame at the highest measured $p_T$, all three polarization parameters are consistent with 0 in both reference frames without any strong $p_T$ dependence. Several model calculations are compared with data, and the one using the Color Glass Condensate effective field theory coupled with non-relativistic QCD gives the best overall description of the experimental results, even though other models cannot be ruled out due to experimental uncertainties.

We report on measurements of sequential $\Upsilon$ suppression in Au+Au collisions at $\sqrt{s_{_\mathrm{NN}}}$ = 200 GeV with the STAR detector at the Relativistic Heavy Ion Collider (RHIC) through both the dielectron and dimuon decay channels. In the 0-60% centrality class, the nuclear modification factors ($R_{\mathrm{AA}}$), which quantify the level of yield suppression in heavy-ion collisions compared to $p$+$p$ collisions, for $\Upsilon$(1S) and $\Upsilon$(2S) are $0.40 \pm 0.03~\textrm{(stat.)} \pm 0.03~\textrm{(sys.)} \pm 0.09~\textrm{(norm.)}$ and $0.26 \pm 0.08~\textrm{(stat.)} \pm 0.02~\textrm{(sys.)} \pm 0.06~\textrm{(norm.)}$, respectively, while the upper limit of the $\Upsilon$(3S) $R_{\mathrm{AA}}$ is 0.17 at a 95% confidence level. This provides experimental evidence that the $\Upsilon$(3S) is significantly more suppressed than the $\Upsilon$(1S) at RHIC. The level of suppression for $\Upsilon$(1S) is comparable to that observed at the much higher collision energy at the Large Hadron Collider. These results point to the creation of a medium at RHIC whose temperature is sufficiently high to strongly suppress excited $\Upsilon$ states.

Measurement by the STAR experiment at RHIC of the cold nuclear matter (CNM) effects experienced by inclusive $J/\psi$ at mid-rapidity in 0-100%$p$+Au collisions at $\sqrt{s_{_{\mathrm{NN}}}}$ = 200 GeV is presented. Such effects are quantified utilizing the nuclear modification factor, $R_{p\mathrm{Au}}$, obtained by taking a ratio of $J/\psi$ yield in $p$+Au collisions to that in $p$+$p$ collisions scaled by the number of binary nucleon-nucleon collisions. The differential $J/\psi$ yield in both $p$+$p$ and $p$+Au collisions is measured through the dimuon decay channel, taking advantage of the trigger capability provided by the Muon Telescope Detector in the RHIC 2015 run. Consequently, the $J/\psi$$R_{p\mathrm{Au}}$ is derived within the transverse momentum ($p_{\mathrm{T}}$) range of 0 to 10 GeV/$c$. A suppression of approximately 30% is observed for $p_{\mathrm{T}}<2$ GeV/$c$, while $J/\psi$ $R_{p\mathrm{Au}}$ becomes compatible with unity for $p_{\mathrm{T}}$ greater than 3 GeV/$c$, indicating the $J/\psi$ yield is minimally affected by the CNM effects at high $p_{\mathrm{T}}$. Comparison to a similar measurement from 0-20% central Au+Au collisions reveals that the observed strong $J/\psi$ suppression above 3 Gev/$c$ is mostly due to the hot medium effects, providing strong evidence for the formation of the quark-gluon plasma in these collisions. Several model calculations show qualitative agreement with the measured $J/\psi$ $R_{p\mathrm{Au}}$, while their agreement with the $J/\psi$ yields in $p$+$p$ and $p$+Au collisions is worse.

We present systematic measurements of azimuthal anisotropy for strange and multistrange hadrons ($K^{0}_{s}$, $\Lambda$, $\Xi$, and $\Omega$) and $\phi$ mesons at midrapidity ($|y| <$ 1.0) in collisions of U + U nuclei at $\sqrt{s_{NN}} = 193$ GeV, recorded by the STAR detector at the Relativistic Heavy Ion Collider. Transverse momentum ($p_{\text{T}}$) dependence of flow coefficients ($v_{2}$, $v_{3}$, and $v_{4}$) is presented for minimum bias collisions and three different centrality intervals. Number of constituent quark scaling of the measured flow coefficients in U + U collisions is discussed. We also present the ratio of $v_{n}$ scaled by the participant eccentricity ($\varepsilon_{n}\left\lbrace 2 \right\rbrace$) to explore system size dependence and collectivity in U + U collisions. The magnitude of $v_{2}/\varepsilon_{2}$ is found to be smaller in U + U collisions than that in central Au + Au collisions contradicting naive eccentricity scaling. Furthermore, the ratios between various flow harmonics ($v_{3}/v_{2}^{3/2}$, $v_{4}/v_{2}^{4/2}$) are studied and compared with hydrodynamic and transport model calculations.