Fluctuations of conserved quantities such as baryon number, charge, and strangeness are sensitive to the correlation length of the hot and dense matter created in relativistic heavy-ion collisions and can be used to search for the QCD critical point. We report the first measurements of the moments of net-kaon multiplicity distributions in Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 7.7, 11.5, 14.5, 19.6, 27, 39, 62.4, and 200 GeV. The collision centrality and energy dependence of the mean ($M$), variance ($\sigma^2$), skewness ($S$), and kurtosis ($\kappa$) for net-kaon multiplicity distributions as well as the ratio $\sigma^2/M$ and the products $S\sigma$ and $\kappa\sigma^2$ are presented. Comparisons are made with Poisson and negative binomial baseline calculations as well as with UrQMD, a transport model (UrQMD) that does not include effects from the QCD critical point. Within current uncertainties, the net-kaon cumulant ratios appear to be monotonic as a function of collision energy.

Bursts of gamma ray showers have been observed in coincidence with downward propagating negative leaders in lightning flashes by the Telescope Array Surface Detector (TASD). The TASD is a 700~square kilometer cosmic ray observatory located in southwestern Utah, U.S.A. In data collected between 2014 and 2016, correlated observations showing the structure and temporal development of three shower-producing flashes were obtained with a 3D lightning mapping array, and electric field change measurements were obtained for an additional seven flashes, in both cases co-located with the TASD. National Lightning Detection Network (NLDN) information was also used throughout. The showers arrived in a sequence of 2--5 short-duration ($\le$10~$\mu$s) bursts over time intervals of several hundred microseconds, and originated at an altitude of $\simeq$3--5 kilometers above ground level during the first 1--2 ms of downward negative leader breakdown at the beginning of cloud-to-ground lightning flashes. The shower footprints, associated waveforms and the effect of atmospheric propagation indicate that the showers consist primarily of downward-beamed gamma radiation. This has been supported by GEANT simulation studies, which indicate primary source fluxes of $\simeq$$10^{12}$--$10^{14}$ photons for $16^{\circ}$ half-angle beams. We conclude that the showers are terrestrial gamma-ray flashes (TGFs), similar to those observed by satellites, but that the ground-based observations are more representative of the temporal source activity and are also more sensitive than satellite observations, which detect only the most powerful TGFs.

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

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<p_T<14$ GeV/$c$ and the total cross section are reported and compared to calculations from the color evaporation model and the non-relativistic Quantum Chromodynamics model. The dependence of $J/\psi$ relative yields in three $p_T$ intervals on charged-particle multiplicity at mid-rapidity is measured for the first time in $p+p$ collisions at $\sqrt{s}$ = 200 GeV and compared with that measured at $\sqrt{s}$ = 7 TeV, PYTHIA8 and EPOS3 Monte Carlo generators, and the Percolation model prediction.

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.

We present measurements of three-particle correlations for various harmonics in Au+Au collisions at energies ranging from $\sqrt{s_{{\rm NN}}}=7.7$ to 200 GeV using the STAR detector. The quantity $\langle\cos(m\phi_1+n\phi_2-(m+n)\phi_3)\rangle$ is evaluated as a function of $\sqrt{s_{{\rm NN}}}$, collision centrality, transverse momentum, $p_T$, pseudo-rapidity difference, $\Delta\eta$, and harmonics ($m$ and $n$). These data provide detailed information on global event properties like the three-dimensional structure of the initial overlap region, the expansion dynamics of the matter produced in the collisions, and the transport properties of the medium. A strong dependence on $\Delta\eta$ is observed for most harmonic combinations consistent with breaking of longitudinal boost invariance. Data reveal changes with energy in the two-particle correlation functions relative to the second-harmonic event-plane and provide ways to constrain models of heavy-ion collisions over a wide range of collision energies.

We present three-particle mixed-harmonic correlations $\la \cos (m\phi_a + n\phi_b - (m+n) \phi_c)\ra$ for harmonics $m,n=1-3$ for charged particles in $\sqrt{s_{NN}}=$200 GeV Au+Au collisions at RHIC. These measurements provide information on the three-dimensional structure of the initial collision zone and are important for constraining models of a subsequent low-viscosity quark-gluon plasma expansion phase. We investigate correlations between the first, second and third harmonics predicted as a consequence of fluctuations in the initial state. The dependence of the correlations on the pseudorapidity separation between particles show hints of a breaking of longitudinal invariance. We compare our results to a number of state-of-the art hydrodynamic calculations with different initial states and temperature dependent viscosities. These measurements provide important steps towards constraining the temperature dependent transport and the longitudinal structure of the initial state at RHIC.

We report on a measurement of the cosmic ray energy spectrum by the Telescope Array Low-Energy Extension (TALE) air fluorescence detector. The TALE air fluorescence detector is also sensitive to the Cherenkov light produced by shower particles. Low energy cosmic rays, in the PeV energy range, are detectable by TALE as "Cherenkov Events". Using these events, we measure the energy spectrum from a low energy of $\sim 2$ PeV to an energy greater than 100 PeV. Above 100 PeV TALE can detect cosmic rays using air fluorescence. This allows for the extension of the measurement to energies greater than a few EeV. In this paper, we will describe the detector, explain the technique, and present results from a measurement of the spectrum using $\sim 1000$ hours of observation. The observed spectrum shows a clear steepening near $10^{17.1}$ eV, along with an ankle-like structure at $10^{16.2}$ eV. These features present important constraints on galactic cosmic rays origin and propagation models. The feature at $10^{17.1}$ eV may also mark the end of the galactic cosmic rays flux and the start of the transition to extra-galactic sources.

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.

An intermediate-scale energy spectrum anisotropy has been found in the arrival directions of ultra-high energy cosmic rays of energies above $10^{19.2}$ eV in the northern hemisphere, using 7 years of data from the Telescope Array surface detector. A relative energy distribution test is done comparing events inside oversampled spherical caps of equal exposure, to those outside, using the Poisson likelihood ratio. The center of maximum significance is at $9^h$$16^m$, $45^{\circ}$. and has a deficit of events with energies $10^{19.2}$$\leq$$E$$<$$10^{19.75}$ eV and an excess for $E$$\geq$$10^{19.75}$ eV. The post-trial probability of this energy anisotropy, appearing by chance anywhere on an isotropic sky, is found by Monte Carlo simulation to be $9$$\times$$10^{-5}$ ($3.74$$\sigma_{global}$).

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.

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.

We present measurements of bulk properties of the matter produced in Au+Au collisions at $\sqrt{s_{NN}}=$ 7.7, 11.5, 19.6, 27, and 39 GeV using identified hadrons ($\pi^\pm$, $K^\pm$, $p$ and $\bar{p}$) from the STAR experiment in the Beam Energy Scan (BES) Program at the Relativistic Heavy Ion Collider (RHIC). Midrapidity ($|y|<$0.1) results for multiplicity densities $dN/dy$, average transverse momenta $\langle p_T \rangle$ and particle ratios are presented. The chemical and kinetic freeze-out dynamics at these energies are discussed and presented as a function of collision centrality and energy. These results constitute the systematic measurements of bulk properties of matter formed in heavy-ion collisions over a broad range of energy (or baryon chemical potential) at RHIC.

The Telescope Array (TA) experiment is located in the western desert of Utah, USA, and observes ultra high energy cosmic rays (UHECRs) in the Northern hemisphere. At the highest energies, $E>10$~EeV, the shape of cosmic ray energy spectrum may carry an imprint of the source density distribution along the line of sight different in different directions of the sky. In this study, we search for such directional variations in the shape of the energy spectrum using events observed with the Telescope Array's surface detector. We divide the TA field of view into two nearly equal-exposure regions: the "on-source" region which we define as $\pm 30^\circ$ of the supergalactic plane containing mostly nearby structures, and the complementary "off-source" region where the sources are further away on average. We compare the UHECR spectra in these regions by fitting them to the broken power law and comparing the resulting parameters. We find that the off-source spectrum has an earlier break at highest energies. The chance probability to obtain such or larger difference in statistically equivalent distributions is estimated as $6.2\pm1.1\times10^{-4}$ ($3.2\sigma$) by a Monte-Carlo simulation. The observed difference in spectra is in a reasonable quantitative agreement with a simplified model that assumes that the UHECR sources trace the galaxy distribution from the 2MRS catalogue, primary particles are protons and the magnetic deflections can be neglected.

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 first measurement of the elliptic anisotropy ($v_2$) of the charm meson $D^0$ at mid-rapidity ($|y|$\,$<$\,1) in Au+Au collisions at \sNN = 200\,GeV. The measurement was conducted by the STAR experiment at RHIC utilizing a new high-resolution silicon tracker. The measured $D^0$ $v_2$ in 0--80\% centrality Au+Au collisions can be described by a viscous hydrodynamic calculation for transverse momentum ($p_{\rm T}$) less than 4\,GeV/$c$. The $D^0$ $v_2$ as a function of transverse kinetic energy ($m_{\rm T} - m_0$, where $m_{\rm T} = \sqrt{p_{\rm T}^2 + m_0^2}$) is consistent with that of light mesons in 10--40\% centrality Au+Au collisions. These results suggest that charm quarks have achieved local thermal equilibrium with the medium created in such collisions. Several theoretical models, with the temperature--dependent, dimensionless charm spatial diffusion coefficient ($2{\pi}TD_s$) in the range of $\sim$2--12, are able to simultaneously reproduce our $D^0$ $v_2$ result and our previously published results for the $D^0$ nuclear modification factor.

We present measurements of elliptic flow ($v_2$) of electrons from the decays of heavy-flavor hadrons ($e_{HF}$) by the STAR experiment. For Au+Au collisions at $\sqrt{s_{\rm NN}} = $ 200 GeV we report $v_2$, for transverse momentum ($p_T$) between 0.2 and 7 GeV/c using three methods: the event plane method ($v_{2}${EP}), two-particle correlations ($v_2${2}), and four-particle correlations ($v_2${4}). For Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 62.4 and 39 GeV we report $v_2${2} for $p_T< 2$ GeV/c. $v_2${2} and $v_2${4} are non-zero at low and intermediate $p_T$ at 200 GeV, and $v_2${2} is consistent with zero at low $p_T$ at other energies. The $v_2${2} at the two lower beam energies is systematically lower than at $\sqrt{s_{\rm NN}} = $ 200 GeV for $p_T < 1$ GeV/c. This difference may suggest that charm quarks interact less strongly with the surrounding nuclear matter at those two lower energies compared to $\sqrt{s_{\rm NN}} = 200$ GeV.

We report the first measurement of the longitudinal double-spin asymmetry $A_{LL}$ for mid-rapidity di-jet production in polarized $pp$ collisions at a center-of-mass energy of $\sqrt{s} = 200$ GeV. The di-jet cross section was measured and is shown to be consistent with next-to-leading order (NLO) perturbative QCD predictions. $A_{LL}$ results are presented for two distinct topologies, defined by the jet pseudorapidities, and are compared to predictions from several recent NLO global analyses. The measured asymmetries, the first such correlation measurements, support those analyses that find positive gluon polarization at the level of roughly 0.2 over the region of Bjorken-$x > 0.05$.

The STAR Collaboration reports on the photoproduction of $\pi^+\pi^-$ pairs in gold-gold collisions at a center-of-mass energy of 200 GeV/nucleon-pair. These pion pairs are produced when a nearly-real photon emitted by one ion scatters from the other ion. We fit the $\pi^+\pi^-$ invariant mass spectrum with a combination of $\rho$ and $\omega$ resonances and a direct $\pi^+\pi^-$ continuum. This is the first observation of the $\omega$ in ultra-peripheral collisions, and the first measurement of $\rho-\omega$ interference at energies where photoproduction is dominated by Pomeron exchange. The $\omega$ amplitude is consistent with the measured $\gamma p\rightarrow \omega p$ cross section, a classical Glauber calculation and the $\omega\rightarrow\pi^+\pi^-$ branching ratio. The $\omega$ phase angle is similar to that observed at much lower energies, showing that the $\rho-\omega$ phase difference does not depend significantly on photon energy. The $\rho^0$ differential cross section $d\sigma/dt$ exhibits a clear diffraction pattern, compatible with scattering from a gold nucleus, with 2 minima visible. The positions of the diffractive minima agree better with the predictions of a quantum Glauber calculation that does not include nuclear shadowing than with a calculation that does include shadowing.

The STAR Collaboration reports the measurement of semi-inclusive distributions of charged-particle jets recoiling from a high transverse momentum hadron trigger, in central and peripheral Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. Charged jets are reconstructed with the anti-kT algorithm for jet radii R between 0.2 and 0.5 and with low infrared cutoff of track constituents ($p_T>0.2$ GeV/c). A novel mixed-event technique is used to correct the large uncorrelated background present in heavy ion collisions. Corrected recoil jet distributions are reported at mid-rapidity, for charged-jet transverse momentum $p_T^\mathrm{jet,ch}<30$ GeV/c. Comparison is made to similar measurements for Pb+Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV, to calculations for p+p collisions at $\sqrt{s}$ = 200 GeV based on the PYTHIA Monte Carlo generator and on a Next-to-Leading Order perturbative QCD approach, and to theoretical calculations incorporating jet quenching. The recoil jet yield is suppressed in central relative to peripheral collisions, with the magnitude of the suppression corresponding to medium-induced charged energy transport out of the jet cone of $2.8\pm0.2\mathrm{(stat)}\pm1.5\mathrm{(sys)}$ GeV/c, for $10<p_T^\mathrm{jet,ch}<20$ GeV/c and R = 0.5. No medium-induced change in jet shape is observed for $R<0.5$. The azimuthal distribution of low-$p_T^\mathrm{jet,ch}$ recoil jets may be enhanced at large azimuthal angles to the trigger axis, due to scattering off quasi-particles in the hot QCD medium. Measurement of this distribution gives a 90% statistical confidence upper limit to the yield enhancement at large deflection angles in central Au+Au collisions of $50\pm30\mathrm{(sys)\%}$ of the large-angle yield in p+pcollisions predicted by PYTHIA.

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 present a measurement of the inclusive production of Upsilon mesons in U+U collisions at 193 GeV at mid-rapidity (|y| < 1). Previous studies in central Au+Au collisions at 200 GeV show a suppression of Upsilon(1S+2S+3S) production relative to expectations from the Upsilon yield in p+p collisions scaled by the number of binary nucleon-nucleon collisions (Ncoll), with an indication that the Upsilon(1S) state is also suppressed. The present measurement extends the number of participant nucleons in the collision (Npart) by 20% compared to Au+Au collisions, and allows us to study a system with higher energy density. We observe a suppression in both the Upsilon(1S+2S+3S) and Upsilon(1S) yields in central U+U data, which consolidates and extends the previously observed suppression trend in Au+Au collisions.