  Complex nature of finite density QCD with heavy quarks in the strong coupling
region is studied. For this purpose, we consider the effective potential as a
function of Polyakov line, and study thermodynamic singularities and associated
Stokes boundaries in the complex chemical potential plane. We also perform an
explicit analytic continuation of the first order transition and crossover
lines in the complex chemical potential plane.
