  Theoretical understanding of experimental results from relativistic heavy-ion
collisions requires a microscopic approach to the behavior of QCD n-point
functions at finite temperatures, as given by the hierarchy of Dyson-Schwinger
equations, properly generalized within the Matsubara formalism. The convergence
of sums over Matsubara modes is studied. The technical complexity of
finite-temperature calculations mandates modeling. We present a model where the
QCD interaction in the infrared, nonperturbative domain, is modeled by a
separable form. Results for the mass spectrum of light quark flavors (u, d, s)
and for the pseudoscalar bound-state amplitudes at finite temperature are
presented.
