Ising model phase boundary

A large number of derivatives of the free energy with respect to the magnetic field are computed for the two-dimensional spin-1/2 Ising model along the phase boundary H=0, T<Tc. This series is found to be divergent, and so the phase boundary is a line of analytic singularities. The same analysis in three dimensions leads to the same conclusion, with lower precision.