Limit[Surd[n!, n] ((Sqrt[Pi] Erfi[Pi^2/6])/2 - (Sqrt[Pi] Erfi[HarmonicNumber[n, 2]])/2),n->∞]