Topological charge density around static colour sources in lattice QCD   with dynamical quarks
  We update our numerical investigation of topological structures around static quarks in pure gauge QCD by results of the first runs including dynamical quarks. Simulations were performed on an $8^3\times4$ lattice, with $SU(3)$ Wilson action, with 3 flavours of quarks of equal mass, both in the confinement and deconfinement phase. In the confinement phase we observe indications for the existence of a flux tube between a static quark and antiquark, flux-tube breaking for large separations, and local correlation between the topological charge density and chiral condensate. In the deconfinement phase almost all configurations turn out to be topologically trivial. 
M. Faber, H. Markum, S. Olejnik and W. Sakuler