Systematic search for islets of stability in the standard map for large parameter values

In the seminal paper (Phys. Rep. 52, 263, 1979), Boris Chirikov showed that the standard map does not exhibit a boundary to chaos, but rather that there are small islands (islets) of stability for arbitrarily large values of the nonlinear perturbation. In this context, he established that the area of the islets in the phase space and the range of parameter values where they exist should decay following power laws with exponents -2 and -1, respectively. In this paper, we carry out a systematic numerical search for islets of stability and we show that the power laws predicted by Chirikov hold. Furthermore, we use high-resolution 3D islets to reveal that the islets volume decays following a similar power law with exponent -3.