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

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.

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.

Density fluctuations near the QCD critical point can be probed via an intermittency analysis in relativistic heavy-ion collisions. We report the first measurement of intermittency in Au$+$Au collisions at $\sqrt{s_\mathrm{_{NN}}}$ = 7.7-200 GeV measured by the STAR experiment at the Relativistic Heavy Ion Collider (RHIC). The scaled factorial moments of identified charged hadrons are analyzed at mid-rapidity and within the transverse momentum phase space. We observe a power-law behavior of scaled factorial moments in Au$+$Au collisions and a decrease in the extracted scaling exponent ($\nu$) from peripheral to central collisions. The $\nu$ is consistent with a constant for different collisions energies in the mid-central (10-40%) collisions. Moreover, the $\nu$ in the 0-5% most central Au$+$Au collisions exhibits a non-monotonic energy dependence that reaches a possible minimum around $\sqrt{s_\mathrm{_{NN}}}$ = 27 GeV. The physics implications on the QCD phase structure are discussed.

We present the first measurement of charge-dependent directed flow in Cu+Au collisions at $\sqrt{s_{_{NN}}}$ = 200 GeV. The results are presented as a function of the particle transverse momentum and pseudorapidity for different centralities. A finite difference between the directed flow of positive and negative charged particles is observed that qualitatively agrees with the expectations from the effects of the initial strong electric field between two colliding ions with different nuclear charges. The measured difference in directed flow is much smaller than that obtained from the parton-hadron-string-dynamics (PHSD) model, which suggests that most of the electric charges, i.e. quarks and antiquarks, have not yet been created during the lifetime of the strong electric field, which is of the order of, or less than, 1fm/$c$.

We report the STAR measurements of dielectron ($e^+e^-$) production at midrapidity ($|y_{ee}|<$1) in Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 200\,GeV. The measurements are evaluated in different invariant mass regions with a focus on 0.30-0.76 ($\rho$-like), 0.76-0.80 ($\omega$-like), and 0.98-1.05 ($\phi$-like) GeV/$c^{2}$. The spectrum in the $\omega$-like and $\phi$-like regions can be well described by the hadronic cocktail simulation. In the $\rho$-like region, however, the vacuum $\rho$ spectral function cannot describe the shape of the dielectron excess. In this range, an enhancement of 1.77$\pm$0.11(stat.)$\pm$0.24(sys.)$\pm$0.33(cocktail) is determined with respect to the hadronic cocktail simulation that excludes the $\rho$ meson. The excess yield in the $\rho$-like region increases with the number of collision participants faster than the $\omega$ and $\phi$ yields. Theoretical models with broadened $\rho$ contributions through interactions with constituents in the hot QCD medium provide a consistent description of the dilepton mass spectra for the measurement presented here and the earlier data at the Super Proton Synchrotron energies.

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.

We report the direct virtual photon invariant yields in the transverse momentum ranges $1\!<\!p_{T}\!<\!3$ GeV/$c$ and $5\!<\!p_T\!<\!10$ GeV/$c$ at mid-rapidity derived from the dielectron invariant mass continuum region $0.10<M_{ee}<0.28$ GeV/$c^{2}$ for 0-80\% minimum-bias Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV. A clear excess in the invariant yield compared to the number-of-binary-collisions ($N_{bin}$) scaled $p+p$ reference is observed in the $p_T$ range $1\!<\!p_{T}\!<\!3$ GeV/$c$. For $p_T\!>6$ GeV/$c$ the production follows $N_{bin}$ scaling. Model calculations with contributions from thermal radiation and initial hard parton scattering are consistent within uncertainties with the direct virtual photon invariant yield.

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 present $\Lambda\Lambda$ correlation measurements in heavy-ion collisions for Au+Au collisions at $\sqrt{s_{NN}}= 200$ GeV using the STAR experiment at the Relativistic Heavy-Ion Collider (RHIC). The Lednick\'{y}-Lyuboshitz analytical model has been used to fit the data to obtain a source size, a scattering length and an effective range. Implications of the measurement of the $\Lambda\Lambda$ correlation function and interaction parameters for di-hyperon searches are discussed.