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.

Global hyperon polarization, $\overline{P}_\mathrm{H}$, in Au+Au collisions over a large range of collision energy, $\sqrt{s_\mathrm{NN}}$, was recently measured and successfully reproduced by hydrodynamic and transport models with intense fluid vorticity of the quark-gluon plasma. While naïve extrapolation of data trends suggests a large $\overline{P}_\mathrm{H}$ as the collision energy is reduced, the behavior of $\overline{P}_\mathrm{H}$ at small $\sqrt{s_\mathrm{NN}}<7.7$ GeV is unknown. Operating the STAR experiment in fixed-target mode, we measured the polarization of $\Lambda$ hyperons along the direction of global angular momentum in Au+Au collisions at $\sqrt{s_\mathrm{NN}}=3$ GeV. The observation of substantial polarization of $4.91\pm0.81(\rm stat.)\pm0.15(\rm syst.)$% in these collisions may require a reexamination of the viscosity of any fluid created in the collision, of the thermalization timescale of rotational modes, and of hadronic mechanisms to produce global polarization.

The Beam Energy Scan (BES) program at the Relativistic Heavy Ion Collider (RHIC) was extended to energies below $\sqrt{\textit{s}_{NN}}$ = 7.7 GeV in 2015 by successful implementation of the fixed-target mode of operation in the STAR (Solenoidal Track At RHIC) experiment. In the fixed-target mode, ions circulate in one ring of the collider and interact with a stationary target at the entrance of the STAR Time Projection Chamber. The first results for Au+Au collisions at $\sqrt{\textit{s}_{NN}}$ = 4.5 GeV are presented, including directed and elliptic flow of identified hadrons, and radii from pion femtoscopy. The proton flow and pion femtoscopy results agree quantitatively with earlier measurements by Alternating Gradient Synchrotron experiments at similar energies. This validates running the STAR experiment in the fixed-target configuration. Pion directed and elliptic flow are presented for the first time at this beam energy. Pion and proton elliptic flow show behavior which hints at constituent quark scaling, but large error bars preclude reliable conclusions. The ongoing second phase of BES (BES-II) will provide fixed-target data sets with 100 times more events at each of several energies down to $\sqrt{\textit{s}_{NN}}$ = 3.0 GeV.

In high-energy heavy-ion collisions, partonic collectivity is evidenced by the constituent quark number scaling of elliptic flow anisotropy for identified hadrons. A breaking of this scaling and dominance of baryonic interactions is found for identified hadron collective flow measurements in $\sqrt{s_{\rm NN}}$ = 3 GeV Au+Au collisions. In this paper, we report measurements of the first- and second-order azimuthal anisotropic parameters, $v_1$ and $v_2$, of light nuclei ($d$, $t$, $^{3}$He, $^{4}$He) produced in $\sqrt{s_{\rm NN}}$ = 3 GeV Au+Au collisions at the STAR experiment. An atomic mass number scaling is found in the measured $v_1$ slopes of light nuclei at mid-rapidity. For the measured $v_2$ magnitude, a strong rapidity dependence is observed. Unlike $v_2$ at higher collision energies, the $v_2$ values at mid-rapidity for all light nuclei are negative and no scaling is observed with the atomic mass number. Calculations by the Jet AA Microscopic Transport Model (JAM), with baryonic mean-field plus nucleon coalescence, are in good agreement with our observations, implying baryonic interactions dominate the collective dynamics in 3 GeV Au+Au collisions at RHIC.

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

We report on the measurements of directed flow $v_1$ and elliptic flow $v_2$ for hadrons ($\pi^{\pm}$, $K^{\pm}$, $K_{S}^0$, $p$, $\phi$, $\Lambda$ and $\Xi^{-}$) from Au+Au collisions at $\sqrt{s_{NN}}$ = 3 GeV and $v_{2}$ for ($\pi^{\pm}$, $K^{\pm}$, $p$ and $\overline{p}$) at 27 and 54.4 GeV with the STAR experiment. While at the two higher energy midcentral collisions the number-of-constituent-quark (NCQ) scaling holds, at 3 GeV the $v_{2}$ at midrapidity is negative for all hadrons and the NCQ scaling is absent. In addition, the $v_1$ slopes at midrapidity for almost all observed hadrons are found to be positive, implying dominant repulsive baryonic interactions. The features of negative $v_2$ and positive $v_1$ slope at 3 GeV can be reproduced with a baryonic mean-field in transport model calculations. These results imply that the medium in such collisions is likely characterized by baryonic interactions.

We report on the first measurements of J/$\psi$ production at very low transverse momentum ($p_{T} <$ 0.2 GeV/c) in hadronic Au+Au collisions at $\sqrt{s_{\rm{NN}}} =$ 200 GeV and U+U collisions at $\sqrt{s_{\rm{NN}}} =$ 193 GeV. Remarkably, the inferred nuclear modification factor of J/$\psi$ at mid-rapidity in Au+Au (U+U) collisions reaches about 24 (52) for $p_{T} <$ 0.05 GeV/c in the 60-80$\%$ collision centrality class. This noteworthy enhancement cannot be explained by hadronic production accompanied by cold and hot medium effects. In addition, the $dN/dt$ distribution of J/$\psi$ for the very low $p_{T}$ range is presented for the first time. The distribution is consistent with that expected from the Au nucleus and shows a hint of interference. Comparison of the measurements to theoretical calculations of coherent production shows that the excess yield can be described reasonably well and reveals a partial disruption of coherent production in semi-central collisions, perhaps due to the violent hadronic interactions. Incorporating theoretical calculations, the results strongly suggest that the dramatic enhancement of J/$\psi$ yield observed at extremely low $p_{T}$ originates from coherent photon-nucleus interactions. In particular, coherently produced J/$\psi$'s in violent hadronic collisions may provide a novel probe of the quark-gluon-plasma.

Partons traversing the strongly interacting medium produced in heavy-ion collisions are expected to lose energy depending on their color charge and mass. We measure the nuclear modification factors for charm- and bottom-decay electrons, defined as the ratio of yields, scaled by the number of binary nucleon-nucleon collisions, in $\sqrt{s_{\rm NN}}$ = 200 GeV Au+Au collisions to $p$+$p$ collisions ($R_{\rm AA}$), or in central to peripheral Au+Au collisions ($R_{\rm CP}$). We find the bottom-decay electron $R_{\rm AA}$ and $R_{\rm CP}$ to be significantly higher than that of charm-decay electrons. Model calculations including mass-dependent parton energy loss in a strongly coupled medium are consistent with the measured data. These observations provide clear evidence of mass ordering of charm and bottom quark energy loss when traversing through the strongly coupled medium created in heavy-ion collisions.

Flow harmonics ($\textit{v}_{n}$) of the Fourier expansion for the azimuthal distributions of hadrons are commonly employed to quantify the azimuthal anisotropy of particle production relative to the collision symmetry planes. While lower order Fourier coefficients ($\textit{v}_{2}$ and $\textit{v}_{3}$) are more directly related to the corresponding eccentricities of the initial state, the higher-order flow harmonics ($\textit{v}_{n>3}$) can be induced by a mode-coupled response to the lower-order anisotropies, in addition to a linear response to the same-order anisotropies. These higher-order flow harmonics and their linear and mode-coupled contributions can be used to more precisely constrain the initial conditions and the transport properties of the medium in theoretical models. The multiparticle azimuthal cumulant method is used to measure the linear and mode-coupled contributions in the higher-order anisotropic flow, the mode-coupled response coefficients, and the correlations of the event plane angles for charged particles as functions of centrality and transverse momentum in Au+Au collisions at nucleon-nucleon center-of-mass energy \roots = 200 GeV. The results are compared to similar LHC measurements as well as to several viscous hydrodynamic calculations with varying initial conditions.

Global polarization of $\Xi$ and $\Omega$ hyperons has been measured for the first time in Au+Au collisions at $\sqrt{s_{_{NN}}}$ = 200 GeV. The measurements of the $\Xi^-$ and $\bar{\Xi}^+$ hyperon polarization have been performed by two independent methods, via analysis of the angular distribution of the daughter particles in the parity violating weak decay $\Xi\rightarrow\Lambda+\pi$, as well as by measuring the polarization of the daughter $\Lambda$-hyperon, polarized via polarization transfer from its parent. The polarization, obtained by combining the results from the two methods and averaged over $\Xi^-$ and $\bar{\Xi}^+$, is measured to be $\langle P_\Xi \rangle = 0.47\pm0.10~({\rm stat.})\pm0.23~({\rm syst.})\,\%$ for the collision centrality 20%-80%. The $\langle P_\Xi \rangle$ is found to be slightly larger than the inclusive $\Lambda$ polarization and in reasonable agreement with a multi-phase transport model (AMPT). The $\langle P_\Xi \rangle$ is found to follow the centrality dependence of the vorticity predicted in the model, increasing toward more peripheral collisions. The global polarization of $\Omega$, $\langle P_\Omega \rangle = 1.11\pm0.87~({\rm stat.})\pm1.97~({\rm syst.})\,\%$ was obtained by measuring the polarization of daughter $\Lambda$ in the decay $\Omega \rightarrow \Lambda + K$, assuming the polarization transfer factor $C_{\Omega\Lambda}=1$.

We report on the first measurement of charm-strange meson $D_s^{\pm}$ production at midrapidity in Au+Au collisions at $\sqrt{s_{_{\rm NN}}}$ = 200 GeV from the STAR experiment. The yield ratio between strange ($D_{s}^{\pm}$) and non-strange ($D^{0}$) open-charm mesons is presented and compared to model calculations. A significant enhancement, relative to a PYTHIA simulation of $p$+$p$ collisions, is observed in the $D_{s}^{\pm}/D^0$ yield ratio in Au+Au collisions over a large range of collision centralities. Model calculations incorporating abundant strange-quark production in the quark-gluon plasma (QGP) and coalescence hadronization qualitatively reproduce the data. The transverse-momentum integrated yield ratio of $D_{s}^{\pm}/D^0$ at midrapidity is consistent with a prediction from a statistical hadronization model with the parameters constrained by the yields of light and strange hadrons measured at the same collision energy. These results suggest that the coalescence of charm quarks with strange quarks in the QGP plays an important role in $D_{s}^{\pm}$ meson production in heavy-ion collisions.

Quark interactions with topological gluon configurations can induce chirality imbalance and local parity violation in quantum chromodynamics. This can lead to electric charge separation along the strong magnetic field in relativistic heavy-ion collisions -- the chiral magnetic effect (CME). We report measurements by the STAR collaboration of a CME-sensitive observable in $p$+Au and $d$+Au collisions at 200 GeV, where the CME is not expected, using charge-dependent pair correlations relative to a third particle. We observe strong charge-dependent correlations similar to those measured in heavy-ion collisions. This bears important implications for the interpretation of the heavy-ion data.

We report the first measurement of the inclusive jet and the dijet longitudinal double-spin asymmetries, $A_{LL}$, at midrapidity in polarized $pp$ collisions at a center-of-mass energy $\sqrt{s} = 510$ GeV. The inclusive jet $A_{LL}$ measurement is sensitive to the gluon helicity distribution down to a gluon momentum fraction of $x\approx 0.015$, while the dijet measurements, separated into four jet-pair topologies, provide constraints on the $x$ dependence of the gluon polarization. Both results are consistent with previous measurements made at $\sqrt{s}= 200$ GeV in the overlapping kinematic region, $x > 0.05$, and show good agreement with predictions from recent next-to-leading order global analyses.

We present the first measurements of transverse momentum spectra of $\pi^{\pm}$, $K^{\pm}$, $p(\bar{p})$ at midrapidity ($|y| < 0.1$) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV with the STAR detector at the Relativistic Heavy Ion Collider (RHIC). The centrality dependence of particle yields, average transverse momenta, particle ratios and kinetic freeze-out parameters are discussed. The results are compared with the published results from Au+Au collisions at $\sqrt{s_{NN}} =$ 200 GeV in STAR. The results are also compared to those from A Multi Phase Transport (AMPT) model.

We report cumulants of the proton multiplicity distribution from dedicated fixed-target Au+Au collisions at 3.0 GeV, measured by the STAR experiment in the kinematic acceptance of rapidity ($y$) and transverse momentum ($p_{\rm T}$) within $-0.5 < y<0$ and $0.4 < p_{\rm T} <2.0 $ GeV/$c$. In the most central 0--5% collisions, a proton cumulant ratio is measured to be $C_4/C_2=-0.85 \pm 0.09 ~(\rm stat.) \pm 0.82 ~(\rm syst.)$, which is less than unity, the Poisson baseline. The hadronic transport UrQMD model reproduces our $C_4/C_2$ in the measured acceptance. Compared to higher energy results and the transport model calculations, the suppression in $C_4/C_2$ is consistent with fluctuations driven by baryon number conservation and indicates an energy regime dominated by hadronic interactions. These data imply that the QCD critical region, if created in heavy-ion collisions, could only exist at energies higher than 3 GeV.

We report high-precision measurements of the longitudinal double-spin asymmetry, $A_{LL}$, for midrapidity inclusive jet and dijet production in polarized $pp$ collisions at a center-of-mass energy of $\sqrt{s}=200\,\mathrm{GeV}$. The new inclusive jet data are sensitive to the gluon helicity distribution, $\Delta g(x,Q^2)$, for gluon momentum fractions in the range from $x \simeq 0.05$ to $x \simeq 0.5$, while the new dijet data provide further constraints on the $x$ dependence of $\Delta g(x,Q^2)$. The results are in good agreement with previous measurements at $\sqrt{s}=200\,\mathrm{GeV}$ and with recent theoretical evaluations of prior world data. Our new results have better precision and thus strengthen the evidence that $\Delta g(x,Q^2)$ is positive for $x > 0.05$.

We report measurements of the longitudinal double-spin asymmetry, $A_{LL}$, for inclusive jet and dijet production in polarized proton-proton collisions at midrapidity and center-of-mass energy $\sqrt{s}$ = 510 GeV, using the high luminosity data sample collected by the STAR experiment in 2013. These measurements complement and improve the precision of previous STAR measurements at the same center-of-mass energy that probe the polarized gluon distribution function at partonic momentum fraction 0.015 $\lesssim x \lesssim$ 0.25. The dijet asymmetries are separated into four jet-pair topologies, which provide further constraints on the $x$ dependence of the polarized gluon distribution function. These measurements are in agreement with previous STAR measurements and with predictions from current next-to-leading order global analyses. They provide more precise data at low dijet invariant mass that will better constraint the shape of the polarized gluon distribution function of the proton.

The STAR collaboration presents jet substructure measurements related to both the momentum fraction and the opening angle within jets in \pp and \AuAu collisions at \sqrtsn $= 200$ GeV. The substructure observables include SoftDrop groomed momentum fraction (\zg), groomed jet radius (\rg), and subjet momentum fraction (\zsj) and opening angle (\tsj). The latter observable is introduced for the first time. Fully corrected subjet measurements are presented for \pp collisions and are compared to leading order Monte Carlo models. The subjet \tsj~distributions reflect the jets leading opening angle and are utilized as a proxy for the resolution scale of the medium in \AuAu collisions. We compare data from \AuAu collisions to those from \pp which are embedded in minimum-bias \AuAu events in order to include the effects of detector smearing and the heavy-ion collision underlying event. The subjet observables are shown to be more robust to the background than \zg~and \rg. We observe no significant modifications of the subjet observables within the two highest-energy, back-to-back jets, resulting in a distribution of opening angles and the splittings that are vacuum-like. We also report measurements of the differential di-jet momentum imbalance ($A_{\rm{J}}$) for jets of varying \tsj. We find no qualitative differences in energy loss signatures for varying angular scales in the range $0.1 < $\tsj $ < 0.3$, leading to the possible interpretation that energy loss in this population of high momentum di-jet pairs, is due to soft medium-induced gluon radiation from a single color-charge as it traverses the medium.

The transversity distribution, which describes transversely polarized quarks in transversely polarized nucleons, is a fundamental component of the spin structure of the nucleon, and is only loosely constrained by global fits to existing semi-inclusive deep inelastic scattering (SIDIS) data. In transversely polarized $p^\uparrow+p$ collisions it can be accessed using transverse polarization dependent fragmentation functions which give rise to azimuthal correlations between the polarization of the struck parton and the final state scalar mesons. This letter reports on spin dependent di-hadron correlations measured by the STAR experiment. The new dataset corresponds to 25 pb$^{-1}$ integrated luminosity of $p^\uparrow+p$ collisions at $\sqrt{s}=500$ GeV, an increase of more than a factor of ten compared to our previous measurement at $\sqrt{s}=200$ GeV. Non-zero asymmetries sensitive to transversity are observed at a $Q^2$ of several hundred GeV and are found to be consistent with the former measurement and a model calculation. %we observe consistent with the former measurement are observed.} We expect that these data will enable an extraction of transversity with comparable precision to current SIDIS datasets but at much higher momentum transfers where subleading effects are suppressed.

In relativistic heavy-ion collisions, a global spin polarization, $P_\mathrm{H}$, of $\Lambda$ and $\bar{\Lambda}$ hyperons along the direction of the system angular momentum was discovered and measured across a broad range of collision energies and demonstrated a trend of increasing $P_\mathrm{H}$ with decreasing $\sqrt{s_{NN}}$. A splitting between $\Lambda$ and $\bar{\Lambda}$ polarization may be possible due to their different magnetic moments in a late-stage magnetic field sustained by the quark-gluon plasma which is formed in the collision. The results presented in this study find no significant splitting at the collision energies of $\sqrt{s_{NN}}=19.6$ and $27$ GeV in the RHIC Beam Energy Scan Phase II using the STAR detector, with an upper limit of $P_{\bar{\Lambda}}-P_{\Lambda}<0.24$% and $P_{\bar{\Lambda}}-P_{\Lambda}<0.35$%, respectively, at a 95% confidence level. We derive an upper limit on the na\"ive extraction of the late-stage magnetic field of $B<9.4\cdot10^{12}$ T and $B<1.4\cdot10^{13}$ T at $\sqrt{s_{NN}}=19.6$ and $27$ GeV, respectively, although more thorough derivations are needed. Differential measurements of $P_\mathrm{H}$ were performed with respect to collision centrality, transverse momentum, and rapidity. With our current acceptance of $|y|<1$ and uncertainties, we observe no dependence on transverse momentum and rapidity in this analysis. These results challenge multiple existing model calculations following a variety of different assumptions which have each predicted a strong dependence on rapidity in this collision-energy range.

We report the first multi-differential measurements of strange hadrons of $K^{-}$, $\phi$ and $\Xi^{-}$ yields as well as the ratios of $\phi/K^-$ and $\phi/\Xi^-$ in Au+Au collisions at ${\sqrt{s_{\rm NN}} = \rm{3\,GeV}}$ with the STAR experiment fixed target configuration at RHIC. The $\phi$ mesons and $\Xi^{-}$ hyperons are measured through hadronic decay channels, $\phi\rightarrow K^+K^-$ and $\Xi^-\rightarrow \Lambda\pi^-$. Collision centrality and rapidity dependence of the transverse momentum spectra for these strange hadrons are presented. The $4\pi$ yields and ratios are compared to thermal model and hadronic transport model predictions. At this collision energy, thermal model with grand canonical ensemble (GCE) under-predicts the $\phi/K^-$ and $\phi/\Xi^-$ ratios while the result of canonical ensemble (CE) calculations reproduce $\phi/K^-$, with the correlation length $r_c \sim 2.7$ fm, and $\phi/\Xi^-$, $r_c \sim 4.2$ fm, for the 0-10% central collisions. Hadronic transport models including high mass resonance decays could also describe the ratios. While thermal calculations with GCE work well for strangeness production in high energy collisions, the change to CE at $\rm{3\,GeV}$ implies a rather different medium property at high baryon density.

The deconfined quark-gluon plasma (QGP) created in relativistic heavy-ion collisions enables the exploration of the fundamental properties of matter under extreme conditions. Non-central collisions can produce strong magnetic fields on the order of $10^{18}$ Gauss, which offers a probe into the electrical conductivity of the QGP. In particular, quarks and anti-quarks carry opposite charges and receive contrary electromagnetic forces that alter their momenta. This phenomenon can be manifested in the collective motion of final-state particles, specifically in the rapidity-odd directed flow, denoted as $v_1(\mathsf{y})$. Here we present the charge-dependent measurements of $dv_1/d\mathsf{y}$ near midrapidities for $\pi^{\pm}$, $K^{\pm}$, and $p(\bar{p})$ in Au+Au and isobar ($_{44}^{96}$Ru+$_{44}^{96}$Ru and $_{40}^{96}$Zr+$_{40}^{96}$Zr) collisions at $\sqrt{s_{\rm NN}}=$ 200 GeV, and in Au+Au collisions at 27 GeV, recorded by the STAR detector at the Relativistic Heavy Ion Collider. The combined dependence of the $v_1$ signal on collision system, particle species, and collision centrality can be qualitatively and semi-quantitatively understood as several effects on constituent quarks. While the results in central events can be explained by the $u$ and $d$ quarks transported from initial-state nuclei, those in peripheral events reveal the impacts of the electromagnetic field on the QGP. Our data put valuable constraints on the electrical conductivity of the QGP in theoretical calculations.

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.

Quark interactions with topological gluon configurations can induce local chirality imbalance and parity violation in quantum chromodynamics, which can lead to the chiral magnetic effect (CME) -- an electric charge separation along the strong magnetic field in relativistic heavy-ion collisions. The CME-sensitive azimuthal correlator observable ($\Delta\gamma$) is contaminated by background arising, in part, from resonance decays coupled with elliptic anisotropy ($v_{2}$). We report here differential measurements of the correlator as a function of the pair invariant mass ($m_{\rm inv}$) in 20-50% centrality Au+Au collisions at $\sqrt{s_{_{\rm NN}}}$= 200 GeV by the STAR experiment at RHIC. Strong resonance background contributions to $\Delta\gamma$ are observed. At large $m_{\rm inv}$ where this background is significantly reduced, the $\Delta\gamma$ value is found to be significantly smaller. An event-shape-engineering technique is deployed to determine the $v_{2}$ background shape as a function of $m_{\rm inv}$. We extract a $v_2$-independent and $m_{\rm inv}$-averaged signal $\Delta\gamma_{\rm sig}$ = (0.03 $\pm$ 0.06 $\pm$ 0.08) $\times10^{-4}$, or $(2\pm4\pm5)\%$ of the inclusive $\Delta\gamma(m_{\rm inv}>0.4$ GeV/$c^2$)$ =(1.58 \pm 0.02 \pm 0.02) \times10^{-4}$, within pion $p_{T}$ = 0.2 - 0.8~\gevc and averaged over pseudorapidity ranges of $-1 < \eta < -0.05$ and $0.05 < \eta < 1$. This represents an upper limit of $0.23\times10^{-4}$, or $15\%$ of the inclusive result, at $95\%$ confidence level for the $m_{\rm inv}$-integrated CME contribution.

The elliptic ($v_2$) and triangular ($v_3$) azimuthal anisotropy coefficients in central $^{3}$He+Au, $d$+Au, and $p$+Au collisions at $\mbox{$\sqrt{s_{\mathrm{NN}}}$}$ = 200 GeV are measured as a function of transverse momentum ($p_{\mathrm{T}}$) at mid-rapidity ($|\eta|<$0.9), via the azimuthal angular correlation between two particles both at $|\eta|<$0.9. While the $v_2(p_{\mathrm{T}})$ values depend on the colliding systems, the $v_3(p_{\mathrm{T}})$ values are system-independent within the uncertainties, suggesting an influence on eccentricity from sub-nucleonic fluctuations in these small-sized systems. These results also provide stringent constraints for the hydrodynamic modeling of these systems.