Harmonic inverse mean curvature surfaces and Painlev equations