We revisit spatially flat FLRW cosmology in light of recent advances in standard model relativistic fluid dynamics. Modern fluid dynamics requires the presence of curvature-matter terms in the energy-momentum tensor for consistency. These terms are linear in the Ricci scalar and tensor, such that the corresponding cosmological model is referred to as ``Ricci cosmology''. No cosmological constant is included, there are no inflaton fields, bulk viscosity is assumed to be zero and we only employ standard Einstein gravity. Analytic solutions to Ricci cosmology are discussed, and we find that it is possible to support an early-time inflationary universe using only well-known ingredients from the Standard Model of physics and geometric properties of space-time.

Ten years ago, relativistic viscous fluid dynamics was formulated from first principles in an effective field theory framework, based entirely on the knowledge of symmetries and long-lived degrees of freedom. In the same year, numerical simulations for the matter created in relativistic heavy-ion collision experiments became first available, providing constraints on the shear viscosity in QCD. The field has come a long way since then. We present the current status of the theory of non-equilibrium fluid dynamics in 2017, including the divergence of the fluid dynamic gradient expansion, resurgence, non-equilibrium attractor solutions, the inclusion of thermal fluctuations as well as their relation to microscopic theories. Furthermore, we review the theory basis for numerical fluid dynamics simulations of relativistic nuclear collisions, and comparison of modern simulations to experimental data for nucleus-nucleus, nucleus-proton and proton-proton collisions.

At high baryon chemical potential mu_B, the equation of state of QCD allows a weak-coupling expansion in the QCD coupling alpha_s. The result is currently known up to and including the full next-to-next-to-leading order (NNLO) alpha_s^2. Starting at this order, the computations are complicated by the modification of particle propagation in a dense medium, which necessitates non-perturbative treatment of the scale Sqrt(alpha_s)*mu_B. In this work, we apply a Hard-Thermal-Loop scheme for capturing the contributions of this scale to the weak-coupling expansion, and use it to determine the leading-logarithm contribution to NNNLO: alpha_s^3 (ln alpha_s)^2. This result is the first improvement to the equation of state of massless cold quark matter in 40 years. The new term is negligibly small, and thus significantly increases our confidence in the applicability of the weak-coupling expansion.

In the limit of short mean free path, relativistic kinetic theory gives rise to hydrodynamics through a systematically improvable gradient expansion. In the present work, a systematically improvable expansion in the opposite limit of large mean free path is considered, describing the dynamics of particles which are almost, but not quite, non-interacting. This non-hydrodynamic "eremitic" expansion does not break down for large gradients, and may be useful in situations where a hydrodynamic treatment is not applicable. As applications, azimuthal anisotropies at high transverse momentum in Pb+Pb and p+Pb collisions at $\sqrt{s}=5.02$ TeV are calculated from the first order eremitic expansion of kinetic theory in the relaxation time approximation.

While QCD appears not to be accurately solvable in the regime of interest for neutron star physics, microscopic calculations are feasible at both low and very high densities. In this work, we propose using the most realistic calculations in these two regimes of nuclear physics and perturbative QCD, and construct equations of state by matching the results requiring thermodynamic consistency. We find that the resulting equations of state --- in contrast to several hadronic ones --- are able to reproduce current observational data on neutron stars without any fine tuning, and allow stable hybrid stars with masses up to 2.1M_{sun}. Using recent observations of star radii, we perform a maximum likelihood analysis to further constrain the equation of state, and in addition show that the effects of rotation on radii and masses should not be neglected in future precision studies.

We perform an O(alpha_s^2) perturbative calculation of the equation of state of cold but dense QCD matter with two massless and one massive quark flavor, finding that perturbation theory converges reasonably well for quark chemical potentials above 1 GeV. Using a running coupling constant and strange quark mass, and allowing for further non-perturbative effects, our results point to a narrow range where absolutely stable strange quark matter may exist. Absent stable strange quark matter, our findings suggest that quark matter in compact star cores becomes confined to hadrons only slightly above the density of atomic nuclei. Finally, we show that equations of state including quark matter lead to hybrid star masses up to M~2M_solar, in agreement with current observations. For strange stars, we find maximal masses of M~2.75M_solar and conclude that confirmed observations of compact stars with M>2M_solar would strongly favor the existence of stable strange quark matter.

We consider the equation of state of QCD at high density and zero temperature in perturbation theory to first order in the coupling constant $\alpha_s$. We compute the thermodynamic potential including the effect of a non-vanishing mass for the strange quark and show that corrections are sizable. Renormalization group running of the coupling and the strange quark mass plays a crucial role. The structure of quark stars is dramatically modified.