  Dynamics of confinement, chiral symmetry breaking and thermal phase
transition is considered at large $N_c$. It is argued that these phenomena are
quantitatively well described within the Gaussian stochastic model of the QCD
vacuum. Selfconsistent equations are written for the field correlators of the
model, yielding important connection between gluonic correlation length of the
vacuum and the string tension. Comparison to other approaches and experimental
and lattice data is given.
