We present results for pseudo-critical temperatures of QCD chiral crossovers at zero and non-zero values of baryon (B), strangeness (S), electric charge (Q), and isospin (I) chemical potentials µ X=B,Q,S,I . The results were obtained using lattice QCD calculations carried out with two degenerate up and down dynamical quarks and a dynamical strange quark, with quark masses corresponding to physical values of pion and kaon masses in the continuum limit. By parameterizing pseudo-critical temperatures as (0)) 4 , we determined κ X 2 and κ X 4 from Taylor expansions of chiral observables in µ X . We obtained a precise result for T c (0) = (156.5 ± 1.5) MeV. For analogous thermal conditions at the chemical freeze-out of relativistic heavy-ion collisions, i.e., µ S (T, µ B ) and µ Q (T, µ B ) fixed from strangeness-neutrality and isospin-imbalance, we found κ B 2 =0.012(4) and κ B 4 =0.000(4). For µ B 300 MeV, the chemical freeze-out takes place in the vicinity of the QCD phase boundary, which coincides with the lines of constant energy density of 0.42(6) GeV/fm 3 and constant entropy density of 3.7(5) fm −3 .

We present new results on up to 6 th order cumulants of net baryon-number fluctuations at small values of the baryon chemical potential, µB, obtained in lattice QCD calculations with physical values of light and strange quark masses. Representation of the Taylor expansions of higher order cumulants in terms of the ratio of the two lowest order cumulants, MB/σ 2, allows for a parameter free comparison with data on net proton-number cumulants obtained by the STAR Collaboration in the Beam Energy Scan at RHIC. We show that recent high statistics data on skewness and kurtosis ratios of net proton-number distributions, obtained at beam energy √ s N N = 54.4 GeV, agree well with lattice QCD results on cumulants of net baryon-number fluctuations close to the pseudo-critical temperature, Tpc(µB), for the chiral transition in QCD. We also present first results from a next-to-leading order expansion of 5 th and 6 th order cumulants on the line of pseudo-critical temperatures.

Instanton-dyons, also known as instanton-monopoles or instanton-quarks, are topological constituents of the instantons at nonzero temperature and holonomy. We perform numerical simulations of the ensemble of interacting dyons for the SU(2) pure gauge theory, using standard Metropolis Monte Carlo and integration over parameter methods. We calculate the free energy as a function of the holonomy (logarithm of the Polyakov line) , the dyon densities, and the Debye mass, and find its minima as a function of those parameters. Unlike previous numerical study of such ensemble, we show that the back reaction on the holonomy potential does generate confinement, provided the density is sufficiently high (or the temperature sufficiently low). We then report various properties of the self-consistent ensembles as a function of temperature.

We present lattice non-relativistic QCD calculations of bottomonia correlation functions at temperatures T 150 − 350 MeV. The correlation functions were computed using extended bottomonia operators, and on background gauge-field configurations for 2+1-flavor QCD having physical kaon and nearly-physical pion masses. We analyzed these correlation functions based on simple theoretically-motivated parameterizations of the corresponding spectral functions. The results of our analyses are compatible with significant in-medium thermal broadening of the ground state Sand P-wave bottomonia. arXiv:1908.08437v1 [hep-lat]

We present the first lattice QCD study of up to 3S and 2P bottomonia at non-zero temperatures. Correlation functions of bottomonia were computed using novel bottomonium operators and a variational technique, within the lattice non-relativistic QCD framework. We analyzed the bottomonium correlation functions based on simple physicallymotivated spectral functions. We found evidence of sequential in-medium modifications, in accordance with the sizes of the bottomonium states.