Using lattice simulations, we show that there is a phase of thermal QCD, where the spectral density ρ(λ) of Dirac operator changes as 1/λ for infrared eigenvalues λ < T .This behavior persists over the entire low energy band we can resolve accurately, over three orders of magnitude on our largest volumes. We propose that in this "IR phase", the well-known non-interacting scale invariance at very short distances (UV, λ → ∞, asymptotic freedom), coexists with very different interacting type of scale invariance at long distances (IR, λ < T ). Such dynamics may be responsible for the unusual fluidity properties of the medium observed at RHIC and LHC. We point out its connection to the physics of Banks-Zaks fixed point, leading to the possibility of massless glueballs in the fluid. Our results lead to the classification of thermal QCD phases in terms of IR scale invariance. The ensuing picture naturally subsumes the standard chiral crossover feature at "T c " ≈ 155 MeV. Its crucial new aspect is the existence of temperature T IR (200 MeV < T IR < 250 MeV) marking the onset of IR phase and possibly a true phase transition. arXiv:1906.08047v1 [hep-lat] 17 Jun 2019 1. Introduction. The study of strongly interacting matter as a function of temperature and baryon density is an active area of theoretical and experimental research (see [1] for recent review). At high energies of colliding heavy nuclei, such as those studied at LHC and the high end of RHIC, baryon densities are small enough so that the results are generally expected, among other things, to shed light on the nature and properties of thermal QCD transition in the early universe. In this regime, it has become widely accepted, largely due to the matured power of lattice QCD [2], that increasing temperature leads to a smooth crossover in properties of thermal strongly interacting matter. On the experimental side, results from RHIC [3-6] and LHC [7] based on modeling the time evolution of collisions in terms of relativistic hydrodynamics, produced a picture of strongly coupled liquid-like medium with extremely low η/s (shear viscosity/entropy density) at high temperatures. In parallel and initially independent developments, similar values of η/s were obtained in highly symmetric and strongly coupled gauge theories with large number of colors, studied by means of their holographic dual [8]. This sparked a flurry of attempts to model the medium seen in the experiments via more refined descriptions of this type.However, the physics of thermal QCD transition(s) and the nature of the discovered liquid-like state of matter are far from settled, even in the limit of vanishing net baryon density (µ = 0), the setting of our interest. Among other things, the currently favored scenario involving a single feature (crossover at"T c ") offers limited room for accommodating a dramatic change from medium described as weakly interacting hadron resonance gas to strongly interacting near-perfect fluid. In this work, we propose a hierarchy of thermal effects in QCD, based on scale inv...