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.

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.

The longitudinal and transverse spin transfers to $\Lambda$ ($\overline{\Lambda}$) hyperons in polarized proton-proton collisions are expected to be sensitive to the helicity and transversity distributions, respectively, of (anti-)strange quarks in the proton, and to the corresponding polarized fragmentation functions. We report improved measurements of the longitudinal spin transfer coefficient, $D_{LL}$, and the transverse spin transfer coefficient, $D_{TT}$, to $\Lambda$ and $\overline{\Lambda}$ in polarized proton-proton collisions at $\sqrt{s}$ = 200 GeV by the STAR experiment at RHIC. The data set includes longitudinally polarized proton-proton collisions with an integrated luminosity of 52 pb$^{-1}$, and transversely polarized proton-proton collisions with a similar integrated luminosity. Both data sets have about twice the statistics of previous results and cover a kinematic range of $|\eta_{\Lambda(\overline{\Lambda})}|$$<$ 1.2 and transverse momentum $p_{T,{\Lambda(\overline{\Lambda})}}$ up to 8 GeV/$c$. We also report the first measurements of the hyperon spin transfer coefficients $D_{LL}$ and $D_{TT}$ as a function of the fractional jet momentum $z$ carried by the hyperon, which can provide more direct constraints on the polarized fragmentation functions.

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.

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.

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

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.

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 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 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 precision measurements of hypernuclei ${}^3_\Lambda \rm{H}$ and ${}^4_\Lambda \rm{H}$ lifetimes obtained from Au+Au collisions at \snn = 3.0 GeV and 7.2 GeV collected by the STAR experiment at RHIC, and the first measurement of ${}^3_\Lambda \rm{H}$ and ${}^4_\Lambda \rm{H}$ mid-rapidity yields in Au+Au collisions at \snn = 3.0 GeV. ${}^3_\Lambda \rm{H}$ and ${}^4_\Lambda \rm{H}$, being the two simplest bound states composed of hyperons and nucleons, are cornerstones in the field of hypernuclear physics. Their lifetimes are measured to be $221\pm15(\rm stat.)\pm19(\rm syst.)$ ps for ${}^3_\Lambda \rm{H}$ and $218\pm6(\rm stat.)\pm13(\rm syst.)$ ps for ${}^4_\Lambda \rm{H}$. The $p_T$-integrated yields of ${}^3_\Lambda \rm{H}$ and ${}^4_\Lambda \rm{H}$ are presented in different centrality and rapidity intervals. It is observed that the shape of the rapidity distribution of ${}^4_\Lambda \rm{H}$ is different for 0--10% and 10--50% centrality collisions. Thermal model calculations, using the canonical ensemble for strangeness, describes the ${}^3_\Lambda \rm{H}$ yield well, while underestimating the ${}^4_\Lambda \rm{H}$ yield. Transport models, combining baryonic mean-field and coalescence (JAM) or utilizing dynamical cluster formation via baryonic interactions (PHQMD) for light nuclei and hypernuclei production, approximately describe the measured ${}^3_\Lambda \rm{H}$ and ${}^4_\Lambda \rm{H}$ yields. Our measurements provide means to precisely assess our understanding of the fundamental baryonic interactions with strange quarks, which can impact our understanding of more complicated systems involving hyperons, such as the interior of neutron stars or exotic hypernuclei.

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.

The STAR collaboration reports the first observation of exclusive rho^0 photo-production, AuAu->AuAu rho^0, and rho^0 production accompanied by mutual nuclear Coulomb excitation, AuAu->Au*Au*rho^0, in ultra-peripheral heavy-ion collisions. The rho^0 have low transverse momenta, consistent with coherent coupling to both nuclei. The cross sections at sqrt(s_NN)=130GeV agree with theoretical predictions treating rho^0 production and Coulomb excitation as independent processes.

We present the first measurement of midrapidity vector meson φ production in Au+Au collisions at RHIC (sNN=130 GeV) from the STAR detector. For the 11% highest multiplicity collisions, the slope parameter from an exponential fit to the transverse mass distribution is T=379±50(stat)±45(syst) MeV, the yield dN/dy=5.73±0.37(stat)±0.69(syst) per event, and the ratio Nφ/Nh− is found to be 0.021±0.001(stat)±0.004(syst). The measured ratio Nφ/Nh− and T for the φ meson at midrapidity do not change for the selected multiplicity bins.

Data from the first physics run at the Relativistic Heavy-Ion Collider at Brookhaven National Laboratory, Au+Au collisions at $\sqrt{s_{NN}}=130$ GeV, have been analyzed by the STAR Collaboration using three-pion correlations with charged pions to study whether pions are emitted independently at freezeout. We have made a high-statistics measurement of the three-pion correlation function and calculated the normalized three-particle correlator to obtain a quantitative measurement of the degree of chaoticity of the pion source. It is found that the degree of chaoticity seems to increase with increasing particle multiplicity.

We report the first measurement of inclusive antiproton production at mid-rapidity in Au+Au collisions at 130 GeV by the STAR experiment at RHIC. The antiproton transverse mass distributions in the measured transverse momentum range of 0.25 < pT < 0.95 GeV/c are found to fall less steeply for more central collisions. The extrapolated antiproton rapidity density is found to scale approximately with the negative hadron multiplicity density.

We report the first observation of $K^{\star}(892)^{0}\to\pi K$ in relativistic heavy ion collisions. The transverse momentum spectrum of $(K^{\star0}+\bar{K}^{\star0})/2$ from central Au+Au collisions at $\sqrt{s_{_{NN}}}=130$ GeV is presented. The ratios of the $K^{\star0}$ yield derived from these data to the yields of negative hadrons, charged kaons, and $\phi$ mesons have been measured in central and minimum bias collisions and compared with model predictions and comparable $e^{+}e^{-}$, $pp$, and $\bar{p}p$ results. The data indicate no dramatic reduction of $K^{\star0}$ production in relativistic heavy ion collisions despite expected losses due to rescattering effects.

The first measurements of light antinucleus production in Au+Au collisions at RHIC are reported. The observed production rates for antideuterons and antihelions are much larger than in lower energy nucleus-nucleus collisions. A coalescence model analysis of the yields indicates that there is little or no increase in the antinucleon freeze-out volume compared to collisions at SPS energy. These analyses also indicate that the antihelion freeze-out volume is smaller than the antideuteron freeze-out volume.

We report the first measurement of strange ($\Lambda$) and anti-strange ($\bar{\Lambda}$) baryon production from $\sqrt{s_{_{NN}}}=130$ GeV Au+Au collisions at the Relativistic Heavy Ion Collider (RHIC). Rapidity density and transverse mass distributions at mid-rapidity are presented as a function of centrality. The yield of $\Lambda$ and $\bar{\Lambda}$ hyperons is found to be approximately proportional to the number of negative hadrons. The production of $\bar{\Lambda}$ hyperons relative to negative hadrons increases very rapidly with transverse momentum. The magnitude of the increase cannot be described by existing hadronic string fragmentation models.

The balance function is a new observable based on the principle that charge is locally conserved when particles are pair produced. Balance functions have been measured for charged particle pairs and identified charged pion pairs in Au+Au collisions at $\sqrt{s_{NN}}$ = 130 GeV at the Relativistic Heavy Ion Collider using STAR. Balance functions for peripheral collisions have widths consistent with model predictions based on a superposition of nucleon-nucleon scattering. Widths in central collisions are smaller, consistent with trends predicted by models incorporating late hadronization.

Mid-rapidity transverse mass spectra and multiplicity densities of charged and neutral kaons are reported for Au+Au collisions at $\snn$=130 GeV at RHIC. The spectra are exponential in transverse mass, with an inverse slope of about 280 MeV in central collisions. The multiplicity densities for these particles scale with the negative hadron pseudo-rapidity density. The charged kaon to pion ratios are $K^+/\pi^- = 0.161 \pm 0.002 {\rm (stat)} \pm 0.024 {\rm (syst)}$ and $K^-/\pi^- = 0.146 \pm 0.002 {\rm (stat)} \pm 0.022 {\rm (syst)}$ for the most central collisions. The $K^+/\pi^-$ ratio is lower than the same ratio observed at the SPS while the $K^-/\pi^-$ is higher than the SPS result. Both ratios are enhanced by about 50% relative to p+p and $\bar{\rm p}$+p collision data at similar energies.

Values of the ratios in the mid-rapidity yields of anti-Lambda/Lambda = 0.71 +/- 0.01(stat.) +/- 0.04(sys.), anti-Xi+/Xi- = 0.83 +/- 0.04(stat.) +/- 0.05 (sys.), anti-Omega+/Omega- = 0.95 +/- 0.15(stat) +/- 0.05(sys.) and K+/K- 1.092 +/- 0.023(combined) were obtained in central sqrt(s_NN) = 130 GeV Au+Au collisions using the STAR detector. The ratios indicate that a fraction of the net-baryon number from the initial system is present in the excess of hyperons over anti-hyperons at mid-rapidity. The trend in the progression of the baryon ratios, with increasing strange quark content, is similar to that observed in heavy-ion collisions at lower energies. The value of these ratios may be related to the charged kaon ratio in the framework of simple quark-counting and thermal models.

Elliptic flow holds much promise for studying the early-time thermalization attained in ultrarelativistic nuclear collisions. Flow measurements also provide a means of distinguishing between hydrodynamic models and calculations which approach the low density (dilute gas) limit. Among the effects that can complicate the interpretation of elliptic flow measurements are azimuthal correlations that are unrelated to the reaction plane (non-flow correlations). Using data for Au + Au collisions at sqrt{s_{NN}} = 130 GeV from the STAR TPC, it is found that four-particle correlation analyses can reliably separate flow and non-flow correlation signals. The latter account for on average about 15% of the observed second-harmonic azimuthal correlation, with the largest relative contribution for the most peripheral and the most central collisions. The results are also corrected for the effect of flow variations within centrality bins. This effect is negligible for all but the most central bin, where the correction to the elliptic flow is about a factor of two. A simple new method for two-particle flow analysis based on scalar products is described. An analysis based on the distribution of the magnitude of the flow vector is also described.