The chiral phase transition at non-zero imaginary baryon chemical potential for different numbers of quark flavours

The so-called Columbia plot summarises the order of the QCD thermal transition as a function of the number of quark flavours and their masses. Recently, it was demonstrated that the first-order chiral transition region, as seen for $N_f\in [3,6]$ on coarse lattices, exhibits tricritical scaling while extrapolating to zero on sufficiently fine lattices. Here we extend these studies to imaginary baryon chemical potential. A similar shrinking of the first-order region is observed with decreasing lattice spacing, which again appears compatible with a tricritical extrapolation to zero.