Analytic continuation of finite density QCD with heavy quarks in the strong coupling region

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.