A status of lattice QCD thermodynamics, as of 2013, is summarized. Only bulk thermodynamics is considered. There is a separate section on magnetic fields.
In this talk we report on our study of two-colour lattice QCD with N_f=4 staggered fermion degrees of freedom with equal electric charge q in a homogeneous magnetic field B at non-zero temperature T. We find indications for a non-monotonic behaviour
of the critical temperature as a function of the magnetic field strength and, as a consequence, for the occurence of `inverse magnetic catalysis within the transition region for magnetic fields in the range 0 < qB < 0.7 GeV^2.
We present results on the QCD equation of state, obtained with two different improved dynamical staggered fermion actions and almost physical quark masses. Lattice cut-off effects are discussed in detail as results for three different lattice spacing
s are available now, i.e. results have been obtained on lattices with temporal extent of $N_tau=4,6$ and 8. Furthermore we discuss the Taylor expansion approach to non-zero baryon chemical potential by means of an expansion of the pressure. We use the expansion coefficients to calculate various fluctuations and correlations among hadronic charges. We find that the correlations reproduce the qualitative behavior of the resonance gas model below $T_c$ and start to agree with the free gas predictions for $Tgsim 1.5T_c$.
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 $mu_{X=B,Q,S,I}$. The results were obtained usin
g 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 $ T_c(mu_X) = T_c(0) left[ 1 -kappa_2^{X}(mu_{X}/T_c(0))^2 -kappa_4^{X}(mu_{X}/T_c(0))^4 right] $, we determined $kappa_2^X$ and $kappa_4^X$ from Taylor expansions of chiral observables in $mu_X$. We obtained a precise result for $T_c(0)=(156.5pm1.5);mathrm{MeV}$. For analogous thermal conditions at the chemical freeze-out of relativistic heavy-ion collisions, i.e., $mu_{S}(T,mu_{B})$ and $mu_{Q}(T,mu_{B})$ fixed from strangeness-neutrality and isospin-imbalance, we found $kappa_2^B=0.012(4)$ and $kappa_4^B=0.000(4)$. For $mu_{B}lesssim300;mathrm{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);mathrm{GeV/fm}^3$ and constant entropy density of $3.7(5);mathrm{fm}^{-3}$.
JLQCD collaboration recently started the $N_f=3$ QCD simulations with the $O(a)$-improved Wilson fermion action employing an exact fermion algorithm developed for odd number of quark flavors. It is found that this theory has an unexpected non-trivial
phase structure in the $(beta,kappa)$ plane even at zero temperature. A detailed study is made to understand the nature of the observed phase transitions and to find the way of avoiding untolerably large lattice artifacts associated with the phase transition.
We present our results obtained from gauge cooled complex Langevin simulations in 1+1d QCD at non-zero densities in the strong coupling regime with unrooted staggered fermions. For small quark masses there are regions of the chemical potential where
this method fails to reproduce correct results. In these parameter ranges we studied the effect of different gauge cooling schemes on the distributions of the fermion determinant as well as of observables.