Natural boundaries for area-preserving twist maps

We consider KAM invariant curves for generalizations of the standard map of the form (x, y)=(x+y, y+f(x)), wheref(x) is an odd trigonometric polynomial. We study numerically their analytic properties by a Padé approximant method applied to the function which conjugates the dynamics to a rotation +. In the complex plane, natural boundaries of different shapes are found. In the complex plane the analyticity region appears to be a strip bounded by a natural boundary, whose width tends linearly to 0 as tends to the critical value.