We report on a polarization measurement of inclusive $J/\psi$ mesons in the di-electron decay channel at mid-rapidity at 2 $<p_{T}<$ 6 GeV/$c$ in $p+p$ collisions at $\sqrt{s}$ = 200 GeV. Data were taken with the STAR detector at RHIC. The $J/\psi$ polarization measurement should help to distinguish between different models of the $J/\psi$ production mechanism since they predict different $p_{T}$ dependences of the $J/\psi$ polarization. In this analysis, $J/\psi$ polarization is studied in the helicity frame. The polarization parameter $\lambda_{\theta}$ measured at RHIC becomes smaller towards high $p_{T}$, indicating more longitudinal $J/\psi$ polarization as $p_{T}$ increases. The result is compared with predictions of presently available models.

The $\jpsi$ $\pt$ spectrum and nuclear modification factor ($\raa$) are reported for $\pt < 5 \ \gevc$ and $|y|<1$ from 0\% to 60\% central Au+Au and Cu+Cu collisions at $\snn = 200 \ \gev$ at STAR. A significant suppression of $\pt$-integrated $\jpsi$ production is observed in central Au+Au events. The Cu+Cu data are consistent with no suppression, although the precision is limited by the available statistics. $\raa$ in Au+Au collisions exhibits a strong suppression at low transverse momentum and gradually increases with $\pt$. The data are compared to high-$\pt$ STAR results and previously published BNL Relativistic Heavy Ion Collider results. Comparing with model calculations, it is found that the invariant yields at low $\pt$ are significantly above hydrodynamic flow predictions but are consistent with models that include color screening and regeneration.

We present a measurement of inclusive $J/\psi$ production at mid-rapidity ($|y|<1$) in $p+p$ collisions at a center-of-mass energy of $\sqrt{s}$ = 200 GeV with the STAR experiment at the Relativistic Heavy Ion Collider (RHIC). The differential production cross section for $J/\psi$ as a function of transverse momentum ($p_T$) for $0

We report the measurement of $K^{*0}$ meson at midrapidity ($|y|<$ 1.0) in Au+Au collisions at $\sqrt{s_{\rm NN}}$~=~7.7, 11.5, 14.5, 19.6, 27 and 39 GeV collected by the STAR experiment during the RHIC beam energy scan (BES) program. The transverse momentum spectra, yield, and average transverse momentum of $K^{*0}$ are presented as functions of collision centrality and beam energy. The $K^{*0}/K$ yield ratios are presented for different collision centrality intervals and beam energies. The $K^{*0}/K$ ratio in heavy-ion collisions are observed to be smaller than that in small system collisions (e+e and p+p). The $K^{*0}/K$ ratio follows a similar centrality dependence to that observed in previous RHIC and LHC measurements. The data favor the scenario of the dominance of hadronic re-scattering over regeneration for $K^{*0}$ production in the hadronic phase of the medium.

The polarizations of the Y(1S), Y(2S), and Y(3S) mesons are measured as a function of the charged particle multiplicity in proton-proton collisions at sqrt(s) = 7 TeV. The measurements are performed with a dimuon data sample collected in 2011 by the CMS experiment, corresponding to an integrated luminosity of 4.9 inverse femtobarns. The results are extracted from the dimuon decay angular distributions, in two ranges of Y(nS) transverse momentum (10-15 and 15-35 GeV), and in the rapidity interval abs(y) < 1.2. The results do not show significant changes from low- to high-multiplicity pp collisions, although large uncertainties preclude definite statements in the Y(2S) and Y(3S) cases.

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 report on the measurement of $\rm{J}/\psi$ production in the dielectron channel at mid-rapidity (|y|<1) in p+p and d+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from the STAR experiment at the Relativistic Heavy Ion Collider. The transverse momentum $p_{T}$ spectra in p+p for $p_{T}$ < 4 GeV/c and d+Au collisions for $p_{T}$ < 3 GeV/c are presented. These measurements extend the STAR coverage for $\rm{J}/\psi$ production in p+p collisions to low $p_{T}$. The $<p_{T}^{2}>$ from the measured $\rm{J}/\psi$ invariant cross section in p+p and d+Au collisions are evaluated and compared to similar measurements at other collision energies. The nuclear modification factor for $\rm{J}/\psi$ is extracted as a function of $p_{T}$ and collision centrality in d+Au and compared to model calculations using the modified nuclear Parton Distribution Function and a final-state $\rm{J}/\psi$ nuclear absorption cross section.

A search is described for the production of a pair of bottom-type vector-like quarks (B VLQs) with mass greater than 1000 GeV. Each B VLQ decays into a b quark and a Higgs boson, a b quark and a Z boson, or a t quark and a W boson. This analysis considers both fully hadronic final states and those containing a charged lepton pair from a Z boson decay. The products of the H $to$ bb boson decay and of the hadronic Z or W boson decays can be resolved as two distinct jets or merged into a single jet, so the final states are classified by the number of reconstructed jets. The analysis uses data corresponding to an integrated luminosity of 138 fb$^{-1}$ collected in proton-proton collisions at $\sqrt{s}$ = 13 TeV with the CMS detector at the LHC from 2016 to 2018. No excess over the expected background is observed. Lower limits are set on the B VLQ mass at 95% confidence level. These depend on the B VLQ branching fractions and are 1570 and 1540 GeV for 100% B $\to$ bH and 100% B $\to$ bZ, respectively. In most cases, the mass limits obtained exceed previous limits by at least 100 GeV.

A search is described for the production of a pair of bottom-type vector-like quarks (VLQs), each decaying into a b or $\mathrm{\bar{b}}$ quark and either a Higgs or a Z boson, with a mass greater than 1000 GeV. The analysis is based on data from proton-proton collisions at a 13 TeV center-of-mass energy recorded at the CERN LHC, corresponding to a total integrated luminosity of 137 fb$^{-1}$. As the predominant decay modes of the Higgs and Z bosons are to a pair of quarks, the analysis focuses on final states consisting of jets resulting from the six quarks produced in the events. Since the two jets produced in the decay of a highly Lorentz-boosted Higgs or Z boson can merge to form a single jet, nine independent analyses are performed, categorized by the number of observed jets and the reconstructed event mode. No signal in excess of the expected background is observed. Lower limits are set on the VLQ mass at 95% confidence level equal to 1570 GeV in the case where the VLQ decays exclusively to a b quark and a Higgs boson, 1390 GeV for when it decays exclusively to a b quark and a Z boson, and 1450 GeV for when it decays equally in these two modes. These limits represent significant improvements over the previously published VLQ limits.

A search for decays to invisible particles of Higgs bosons produced in association with a top-antitop quark pair or a vector boson, which both decay to a fully hadronic final state, has been performed using proton-proton collision data collected at $\sqrt{s}$ = 13 TeV by the CMS experiment at the LHC, corresponding to an integrated luminosity of 138 fb$^{-1}$. The 95% confidence level upper limit set on the branching fraction of the 125 GeV Higgs boson to invisible particles, $\mathcal{B}$(H $\to$ inv), is 0.54 (0.39 expected), assuming standard model production cross sections. The results of this analysis are combined with previous $\mathcal{B}$(H $\to$ inv) searches carried out at $\sqrt{s}$ = 7, 8, and 13 TeV in complementary production modes. The combined upper limit at 95% confidence level on $\mathcal{B}$(H $\to$ inv) is 0.15 (0.08 expected).

A search for new phenomena is performed in final states containing one or more jets and an imbalance in transverse momentum in pp collisions at a centre-of-mass energy of 13 TeV. The analysed data sample, recorded with the CMS detector at the CERN LHC, corresponds to an integrated luminosity of 2.3 inverse femtobarns. Several kinematic variables are employed to suppress the dominant background, multijet production, as well as to discriminate between other standard model and new physics processes. The search provides sensitivity to a broad range of new-physics models that yield a stable weakly interacting massive particle. The number of observed candidate events is found to agree with the expected contributions from standard model processes, and the result is interpreted in the mass parameter space of fourteen simplified supersymmetric models that assume the pair production of gluinos or squarks and a range of decay modes. For models that assume gluino pair production, masses up to 1575 and 975 GeV are excluded for gluinos and neutralinos, respectively. For models involving the pair production of top squarks and compressed mass spectra, top squark masses up to 400 GeV are excluded.

A search for new light bosons decaying into muon pairs is presented using a data sample corresponding to an integrated luminosity of 35.9 fb$^{-1}$ of proton-proton collisions at a center-of-mass energy $\sqrt{s} =$ 13 TeV, collected with the CMS detector at the CERN LHC. The search is model independent, only requiring the pair production of a new light boson and its subsequent decay to a pair of muons. No significant deviation from the predicted background is observed. A model independent limit is set on the product of the production cross section times branching fraction to dimuons squared times acceptance as a function of new light boson mass. This limit varies between 0.16 and 0.45 fb over a range of new light boson masses from 0.25 to 8.5 GeV. It is then interpreted in the context of the next-to-minimal supersymmetric standard model and a dark supersymmetry model that allows for nonnegligible light boson lifetimes. In both cases, there is significant improvement over previously published limits.

The differential cross sections for the production of photons in Z to mu+ mu- gamma decays are presented as a function of the transverse energy of the photon and its separation from the nearest muon. The data for these measurements were collected with the CMS detector and correspond to an integrated luminosity of 4.7 inverse femtobarns of pp collisions at sqrt(s) = 7 TeV delivered by the CERN LHC. The cross sections are compared to simulations with POWHEG and PYTHIA, where PYTHIA is used to simulate parton showers and final-state photons. These simulations match the data to better than 5%.

The first measurement of the CP structure of the Yukawa coupling between the Higgs boson and $\tau$ leptons is presented. The measurement is based on data collected in proton-proton collisions at $\sqrt{s} =$ 13 TeV by the CMS detector at the LHC, corresponding to an integrated luminosity of 137 fb$^{-1}$. The analysis uses the angular correlation between the decay planes of $\tau$ leptons produced in Higgs boson decays. The effective mixing angle between CP-even and CP-odd $\tau$ Yukawa couplings is found to be $-$1 $\pm$ 19$^\circ$, compared to an expected value of 0 $\pm$ 21$^\circ$ at the 68.3% confidence level. The data disfavour the pure CP-odd scenario at 3.0 standard deviations. The results are compatible with predictions for the standard model Higgs boson.

Angular distributions of the decay B$^+$$\to$ K$^*$(892)$^+\mu^+\mu^-$ are studied using events collected with the CMS detector in $\sqrt{s} =$ 8 TeV proton-proton collisions at the LHC, corresponding to an integrated luminosity of 20.0 fb$^{-1}$. The forward-backward asymmetry of the muons and the longitudinal polarization of the K$^*$(892)$^+$ meson are determined as a function of the square of the dimuon invariant mass. These are the first results from this exclusive decay mode and are in agreement with a standard model prediction.

The angular distribution of the flavor-changing neutral current decay B$^+$$\to$ K$^+\mu^+\mu^-$ is studied in proton-proton collisions at a center-of-mass energy of 8 TeV. The analysis is based on data collected with the CMS detector at the LHC, corresponding to an integrated luminosity of 20.5 fb$^{-1}$. The forward-backward asymmetry $A_{\mathrm{FB}}$ of the dimuon system and the contribution $F_{\mathrm{H}}$ from the pseudoscalar, scalar, and tensor amplitudes to the decay width are measured as a function of the dimuon mass squared. The measurements are consistent with the standard model expectations.

The angular distributions and the differential branching fraction of the decay B0 to K*0(892) mu mu are studied using data corresponding to an integrated luminosity of 20.5 inverse femtobarns collected with the CMS detector at the LHC in pp collisions at sqrt(s) = 8 TeV. From 1430 signal decays, the forward-backward asymmetry of the muons, the K*0(892) longitudinal polarization fraction, and the differential branching fraction are determined as a function of the dimuon invariant mass squared. The measurements are among the most precise to date and are in good agreement with standard model predictions.

Measurements of the five most significant angular coefficients, A[0] through A[4], for Z bosons produced in pp collisions at $\sqrt{s}$ = 8 TeV and decaying to $\mu^+ \mu^-$ are presented as a function of the transverse momentum and rapidity of Z boson. The integrated luminosity of the dataset collected with the CMS detector at the LHC corresponds to 19.7 inverse femtobarns. These measurements provide comprehensive information about Z boson production mechanisms, and are compared to QCD predictions at leading order, next-to-leading order, and next-to-next-to-leading order in perturbation theory.

We present two-dimensional (2D) two-particle angular correlations on relative pseudorapidity $\eta$ and azimuth $\phi$ for charged particles from Au-Au collisions at $\sqrt{s_{\rm NN}} = 62$ and 200 GeV with transverse momentum $p_t \geq 0.15$ GeV/$c$, $|\eta| \leq 1$ and $2\pi$ azimuth. Observed correlations include a {same-side} (relative azimuth $< \pi/2$) 2D peak, a closely-related away-side azimuth dipole, and an azimuth quadrupole conventionally associated with elliptic flow. The same-side 2D peak and away-side dipole are explained by semihard parton scattering and fragmentation (minijets) in proton-proton and peripheral nucleus-nucleus collisions. Those structures follow N-N binary-collision scaling in Au-Au collisions until mid-centrality where a transition to a qualitatively different centrality trend occurs within a small centrality interval. Above the transition point the number of same-side and away-side correlated pairs increases rapidly {relative to} binary-collision scaling, the $\eta$ width of the same-side 2D peak also increases rapidly ($\eta$ elongation) and the $\phi$ width actually decreases significantly. Those centrality trends are more remarkable when contrasted with expectations of jet quenching in a dense medium. Observed centrality trends are compared to {\sc hijing} predictions and to the expected trends for semihard parton scattering and fragmentation in a thermalized opaque medium. We are unable to reconcile a semihard parton scattering and fragmentation origin for the observed correlation structure and centrality trends with heavy ion collision scenarios which invoke rapid parton thermalization. On the other hand, if the collision system is effectively opaque to few-GeV partons the observations reported here would be inconsistent with a minijet picture.

Results on high transverse momentum charged particle emission with respect to the reaction plane are presented for Au+Au collisions at $\sqrt{s_{_{NN}}}$= 200 GeV. Two- and four-particle correlations results are presented as well as a comparison of azimuthal correlations in Au+Au collisions to those in $p+p$ at the same energy. Elliptic anisotropy, $v_2$, is found to reach its maximum at $p_t \sim 3$ GeV/c, then decrease slowly and remain significant up to $p_t\approx 7$ -- 10 GeV/c. Stronger suppression is found in the back-to-back high-$p_t$ particle correlations for particles emitted out-of-plane compared to those emitted in-plane. The centrality dependence of $v_2$ at intermediate $p_t$ is compared to simple models based on jet quenching.

The results from the STAR Collaboration on directed flow (v_1), elliptic flow (v_2), and the fourth harmonic (v_4) in the anisotropic azimuthal distribution of particles from Au+Au collisions at sqrtsNN = 200 GeV are summarized and compared with results from other experiments and theoretical models. Results for identified particles are presented and fit with a Blast Wave model. Different anisotropic flow analysis methods are compared and nonflow effects are extracted from the data. For v_2, scaling with the number of constituent quarks and parton coalescence is discussed. For v_4, scaling with v_2^2 and quark coalescence is discussed.

The azimuthal anisotropic flow of identified and unidentified charged particles has been systematically studied in Cu+Au collisions at $\sqrt{s_{_{NN}}}$ = 200 GeV for harmonics $n=$ 1-4 in the pseudorapidity range $|\eta|<1$. The directed flow in Cu+Au collisions is compared with the rapidity-odd and, for the first time, the rapidity-even components of charged particle directed flow in Au+Au collisions at $\sqrt{s_{_{NN}}}$ = 200~GeV. The slope of the directed flow pseudorapidity dependence in Cu+Au collisions is found to be similar to that in Au+Au collisions, with the intercept shifted toward positive $\eta$ values, i.e., the Cu-going direction. The mean transverse momentum projected onto the spectator plane, $\langle p_x\rangle$, in Cu+Au collision also exhibits approximately linear dependence on $\eta$ with the intercept at about $\eta\approx-0.4$, closer to the rapidity of the Cu+Au system center-of-mass. The observed dependencies find natural explanation in a picture of the directed flow originating partly due the "tilted source" and partly due to the rapidity dependent asymmetry in the initial density distribution. Charge-dependence of the $\langle p_x\rangle$ was also observed in Cu+Au collisions, indicating an effect of the initial electric field created by charge difference of the spectator protons in two colliding nuclei. The rapidity-even component of directed flow in Au+Au collisions is close to that in Pb+Pb collisions at $\sqrt{s_{_{NN}}}$ = 2.76 TeV, indicating a similar magnitude of dipole-like fluctuations in the initial-state density distribution. Higher harmonic flow in Cu+Au collisions exhibits similar trends to those observed in Au+Au and Pb+Pb collisions and is qualitatively reproduced by a viscous hydrodynamic model and a multi-phase transport model. For all harmonics with $n\ge2$ we observe an approximate scaling of $v_n$ with the number of constituent quarks.

Collisions between prolate uranium nuclei are used to study how particle production and azimuthal anisotropies depend on initial geometry in heavy-ion collisions. We report the two- and four-particle cumulants, $v_2\{2\}$ and $v_2\{4\}$, for charged hadrons from U+U collisions at $\sqrt{s_{\rm NN}}$ = 193 GeV and Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 200 GeV. Nearly fully overlapping collisions are selected based on the amount of energy deposited by spectators in the STAR Zero Degree Calorimeters (ZDCs). Within this sample, the observed dependence of $v_2\{2\}$ on multiplicity demonstrates that ZDC information combined with multiplicity can preferentially select different overlap configurations in U+U collisions. An initial-state model with gluon saturation describes the slope of $v_2\{2\}$ as a function of multiplicity in central collisions better than one based on Glauber with a two-component multiplicity model.

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

We report STAR results on the azimuthal anisotropy parameter v2 for strange particles K0S, L and Lbar at midrapidity in Au+Au collisions at sNN = 130 GeV at RHIC. The value of v2 as a function of transverse momentum of the produced particles pt and collision centrality is presented for both particles up to pt 3.0 GeV/c. A strong pt dependence in v2 is observed up to 2.0 GeV/c. The v2 measurement is compared with hydrodynamic model calculations. The physics implications of the pt integrated v2 magnitude as a function of particle mass are also discussed.

The Fourier coefficients v[2] and v[3] characterizing the anisotropy of the azimuthal distribution of charged particles produced in PbPb collisions at sqrt(s[NN]) = 5.02 TeV are measured with data collected by the CMS experiment. The measurements cover a broad transverse momentum range, 1 < pT < 100 GeV. The analysis focuses on pT > 10 GeV range, where anisotropic azimuthal distributions should reflect the path-length dependence of parton energy loss in the created medium. Results are presented in several bins of PbPb collision centrality, spanning the 60% most central events. The v[2] coefficient is measured with the scalar product and the multiparticle cumulant methods, which have different sensitivities to the initial-state fluctuations. The values of both methods remain positive up to pT of about 60-80 GeV, in all examined centrality classes. The v[3] coefficient, only measured with the scalar product method, tends to zero for pT greater than or equal to 20 GeV. Comparisons between theoretical calculations and data provide new constraints on the path-length dependence of parton energy loss in heavy ion collisions and highlight the importance of the initial-state fluctuations.

The path-length dependent parton energy loss within the dense partonic medium created in lead-lead collisions at a nucleon-nucleon center-of-mass energy of $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV is studied by determining the azimuthal anisotropies for dijets with high transverse momentum. The data were collected by the CMS experiment in 2018 and correspond to an integrated luminosity of 1.69 nb$^{-1}$. For events containing back-to-back jets, correlations in relative azimuthal angle and pseudorapidity ($\eta$) between jets and hadrons, and between two hadrons, are constructed. The anisotropies are expressed as the Fourier expansion coefficients $v_n$, $n = $ 2-4 of these azimuthal distributions. The dijet $v_n$ values are extracted from long-range (1.5 $\lt$$\vert\Delta\eta\vert$$\lt$ 2.5) components of these correlations, which suppresses the background contributions from jet fragmentation processes. Positive dijet $v_2$ values are observed which increase from central to more peripheral events, while the $v_3$ and $v_4$ values are consistent with zero within experimental uncertainties.

Azimuthal correlations between the two jets with the largest transverse momenta $ {p_{\mathrm{T}}} $ in inclusive 2-, 3-, and 4-jet events are presented for several regions of the leading jet $ {p_{\mathrm{T}}} $ up to 4 TeV. For 3- and 4-jet scenarios, measurements of the minimum azimuthal angles between any two of the three or four leading $ {p_{\mathrm{T}}} $ jets are also presented. The analysis is based on data from proton-proton collisions collected by the CMS Collaboration at a centre-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 35.9 fb$^{-1}$. Calculations based on leading-order matrix elements supplemented with parton showering and hadronization do not fully describe the data, so next-to-leading-order calculations matched with parton shower and hadronization models are needed to better describe the measured distributions. Furthermore, we show that azimuthal jet correlations are sensitive to details of the parton showering, hadronization, and multiparton interactions. A next-to-leading-order calculation matched with parton showers in the MC@NLO method, as implemented in HERWIG 7, gives a better overall description of the measurements than the POWHEG method.

The production of Z bosons associated with jets is measured in pp collisions at $\sqrt{s}$ = 13 TeV with data recorded with the CMS experiment at the LHC corresponding to an integrated luminosity of 36.3 fb$^{-1}$. The multiplicity of jets with transverse momentum $p_\mathrm{T}$$\gt$ 30 GeV is measured for different regions of the Z boson's $p_\mathrm{T}$(Z), from lower than 10 GeV to higher than 100 GeV. The azimuthal correlation $\Delta \phi$ between the Z boson and the leading jet, as well as the correlations between the two leading jets are measured in three regions of $p_\mathrm{T}$(Z). The measurements are compared with several predictions at leading and next-to-leading orders, interfaced with parton showers. Predictions based on transverse-momentum dependent parton distributions and corresponding parton showers give a good description of the measurement in the regions where multiple parton interactions and higher jet multiplicities are not important. The effects of multiple parton interactions are shown to be important to correctly describe the measured spectra in the low $p_\mathrm{T}$(Z) regions.

The structure of nucleons is multidimensional and depends on the transverse momenta, spatial geometry, and polarization of the constituent partons. Such a structure can be studied using high-energy photons produced in ultraperipheral heavy-ion collisions. The first measurement of the azimuthal angular correlations of exclusively produced events with two jets in photon-lead interactions at large momentum transfer is presented, a process that is considered to be sensitive to the underlying nuclear gluon polarization. This study uses a data sample of ultraperipheral lead-lead collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV, corresponding to an integrated luminosity of 0.38 nb$^{-1}$, collected with the CMS experiment at the LHC. The measured second harmonic of the correlation between the sum and difference of the two jet momenta is found to be positive, and rising, as the dijet momentum increases. A well-tuned model that has been successful at describing a wide range of proton scattering data from the HERA experiments fails to describe the observed correlations, suggesting the presence of gluon polarization effects.

Yields, correlation shapes, and mean transverse momenta \pt{} of charged particles associated with intermediate to high-\pt{} trigger particles ($2.5 < \pt < 10$ \GeVc) in d+Au and Au+Au collisions at $\snn=200$ GeV are presented. For associated particles at higher $\pt \gtrsim 2.5$ \GeVc, narrow correlation peaks are seen in d+Au and Au+Au, indicating that the main production mechanism is jet fragmentation. At lower associated particle $\pt < 2$ \GeVc, a large enhancement of the near- ($\dphi \sim 0$) and away-side ($\dphi \sim \pi$) associated yields is found, together with a strong broadening of the away-side azimuthal distributions in Au+Au collisions compared to d+Au measurements, suggesting that other particle production mechanisms play a role. This is further supported by the observed significant softening of the away-side associated particle yield distribution at $\dphi \sim \pi$ in central Au+Au collisions.

A measurement for inclusive 2- and 3-jet events of the azimuthal correlation between the two jets with the largest transverse momenta, $\Delta\phi_{12}$ , is presented. The measurement considers events where the two leading jets are nearly collinear ("back-to-back") in the transverse plane and is performed for several ranges of the leading jet transverse momentum. Proton-proton collision data collected with the CMS experiment at a center-of-mass energy of 13 TeV and corresponding to an integrated luminosity of 35.9 fb$^{-1}$ are used. Predictions based on calculations using matrix elements at leading-order and next-to-leading-order accuracy in perturbative quantum chromodynamics supplemented with leading-log parton showers and hadronization are generally in agreement with the measurements. Discrepancies between the measurement and theoretical predictions are as large as 15%, mainly in the region 177$^\circ$ $<$ $\Delta\phi_{12}$ $<$ 180$^\circ$. The 2- and 3-jet measurements are not simultaneously described by any of models.

We report the first measurements of transverse single-spin asymmetries for inclusive jet and jet + $\pi^{\pm}$ production at midrapidity from transversely polarized proton-proton collisions at $\sqrt{s} = 500$ GeV. The data were collected in 2011 with the STAR detector sampled from 23 pb$^{-1}$ integrated luminosity with an average beam polarization of 53%. Asymmetries are reported for jets with transverse momenta $6 < p_{T, jet} < 55$ GeV/c and pseudorapidity $|\eta| < 1$. Presented are measurements of the inclusive-jet azimuthal transverse single-spin asymmetry, sensitive to twist-3 initial-state quark-gluon correlators; the Collins asymmetry, sensitive to quark transversity coupled to the polarized Collins fragmentation function; and the first measurement of the "Collins-like" asymmetry, sensitive to linearly polarized gluons. Within the present statistical precision, inclusive-jet and Collins-like asymmetries are small, with the latter allowing the first experimental constraints on gluon linear polarization in a polarized proton. At higher values of jet transverse momenta, we observe the first non-zero Collins asymmetries in polarized-proton collisions, with a statistical significance of greater than $5\sigma$. The results span a range of x similar to results from SIDIS but at much higher $Q^{2}$. The Collins results enable tests of universality and factorization-breaking in the transverse momentum-dependent formulation of perturbative quantum chromodynamics.

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 results of a systematic study of the shape of the pion distribution in coordinate space at freeze-out in Au+Au collisions at RHIC using two-pion Hanbury Brown-Twiss (HBT) interferometry. Oscillations of the extracted HBT radii vs. emission angle indicate sources elongated perpendicular to the reaction plane. The results indicate that the pressure and expansion time of the collision system are not sufficient to completely quench its initial shape.

Balance functions have been measured for charged particle pairs, identified charged pion pairs, and identified charged kaon pairs in Au+Au, d+Au, and p+p collisions at $\sqrt{s_{NN}}$ = 200 GeV at the Relativistic Heavy Ion Collider using the STAR detector. These balance functions are presented in terms of relative pseudorapidity, $\Delta \eta$, relative rapidity, $\Delta y$, relative azimuthal angle, $\Delta \phi$, and invariant relative momentum, $q_{\rm inv}$. In addition, balance functions are shown in terms of the three components of $q_{\rm inv}$: $q_{\rm long}$, $q_{\rm out}$, and $q_{\rm side}$. For charged particle pairs, the width of the balance function in terms of $\Delta \eta$ scales smoothly with the number of participating nucleons, while HIJING and UrQMD model calculations show no dependence on centrality or system size. For charged particle and charged pion pairs, the balance functions widths in terms of $\Delta \eta$ and $\Delta y$ are narrower in central Au+Au collisions than in peripheral collisions. The width for central collisions is consistent with thermal blast-wave models where the balancing charges are highly correlated in coordinate space at breakup. This strong correlation might be explained either by delayed hadronization or by limited diffusion during the reaction. Furthermore, the narrowing trend is consistent with the lower kinetic temperatures inherent to more central collisions. In contrast, the width of the balance function for charged kaon pairs in terms of $\Delta y$ shows little centrality dependence, which may signal a different production mechanism for kaons. The widths of the balance functions for charged pions and kaons in terms of $q_{\rm inv}$ narrow in central collisions compared to peripheral collisions, which may be driven by the change in the kinetic temperature.

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 measurements of the nuclear modification factor, $R_{ \mathrm{CP}}$, for charged hadrons as well as identified $\pi^{+(-)}$, $K^{+(-)}$, and $p(\overline{p})$ for Au+Au collision energies of $\sqrt{s_{_{ \mathrm{NN}}}}$ = 7.7, 11.5, 14.5, 19.6, 27, 39, and 62.4 GeV. We observe a clear high-$p_{\mathrm{T}}$ net suppression in central collisions at 62.4 GeV for charged hadrons which evolves smoothly to a large net enhancement at lower energies. This trend is driven by the evolution of the pion spectra, but is also very similar for the kaon spectra. While the magnitude of the proton $R_{ \mathrm{CP}}$ at high $p_{\mathrm{T}}$ does depend on collision energy, neither the proton nor the anti-proton $R_{ \mathrm{CP}}$ at high $p_{\mathrm{T}}$ exhibit net suppression at any energy. A study of how the binary collision scaled high-$p_{\mathrm{T}}$ yield evolves with centrality reveals a non-monotonic shape that is consistent with the idea that jet-quenching is increasing faster than the combined phenomena that lead to enhancement.

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

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 report the first measurements of the moments -- mean ($M$), variance ($\sigma^{2}$), skewness ($S$) and kurtosis ($\kappa$) -- of the net-charge multiplicity distributions at mid-rapidity in Au+Au collisions at seven energies, ranging from $\sqrt {{s_{\rm NN}}}$= 7.7 to 200 GeV, as a part of the Beam Energy Scan program at RHIC. The moments are related to the thermodynamic susceptibilities of net-charge, and are sensitive to the proximity of the QCD critical point. We compare the products of the moments, $\sigma^{2}/M$, $S\sigma$ and $\kappa\sigma^{2}$ with the expectations from Poisson and negative binomial distributions (NBD). The $S\sigma$ values deviate from Poisson and are close to NBD baseline, while the $\kappa\sigma^{2}$ values tend to lie between the two. Within the present uncertainties, our data do not show non-monotonic behavior as a function of collision energy. These measurements provide a distinct way of determining the freeze-out parameters in heavy-ion collisions by comparing with theoretical models.

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.

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

We present results of analyses of two-pion interferometry in Au+Au collisions at $\sqrt{s_{NN}}$ = 7.7, 11.5, 19.6, 27, 39, 62.4 and 200 GeV measured in the STAR detector as part of the RHIC Beam Energy Scan program. The extracted correlation lengths (HBT radii) are studied as a function of beam energy, azimuthal angle relative to the reaction plane, centrality, and transverse mass ($m_{T}$) of the particles. The azimuthal analysis allows extraction of the eccentricity of the entire fireball at kinetic freeze-out. The energy dependence of this observable is expected to be sensitive to changes in the equation of state. A new global fit method is studied as an alternate method to directly measure the parameters in the azimuthal analysis. The eccentricity shows a monotonic decrease with beam energy that is qualitatively consistent with the trend from all model predictions and quantitatively consistent with a hadronic transport model.

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.

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.

Rapidity-odd directed flow($v_1$) measurements for charged pions, protons and antiprotons near mid-rapidity ($y=0$) are reported in $\sqrt{s_{NN}} =$ 7.7, 11.5, 19.6, 27, 39, 62.4 and 200 GeV Au + Au collisions as recorded by the STAR detector at the Relativistic Heavy Ion Collider (RHIC). At intermediate impact parameters, the proton and net-proton slope parameter $dv_1/dy|_{y=0}$ shows a minimum between 11.5 and 19.6 GeV. In addition, the net-proton $dv_1/dy|_{y=0}$ changes sign twice between 7.7 and 39 GeV. The proton and net-proton results qualitatively resemble predictions of a hydrodynamic model with a first-order phase transition from hadronic matter to deconfined matter, and differ from hadronic transport calculations.

We present a measurement of the first-order azimuthal anisotropy, $v_1(y)$, of deuterons from Au+Au collisions at $\sqrt{s_{NN}}$ = 7.7, 11.5, 14.5, 19.6, 27, and 39 GeV recorded with the STAR experiment at the Relativistic Heavy Ion Collider (RHIC). The energy dependence of the $v_1(y)$ slope, $dv_{1}/dy|_{y=0}$, for deuterons, where $y$ is the rapidity, is extracted for semi-central collisions (10-40\% centrality) and compared to that of protons. While the $v_1(y)$ slopes of protons are generally negative for $\sqrt{s_{NN}} >$ 10 GeV, those for deuterons are consistent with zero, a strong enhancement of the $v_1(y)$ slope of deuterons is seen at the lowest collision energy (the largest baryon density) at $\sqrt{s_{NN}} =$ 7.7 GeV. In addition, we report the transverse momentum dependence of $v_1$ for protons and deuterons. The experimental results are compared with transport and coalescence models.