Near-axis description of stellarator-symmetric quasi-isodynamic stellarators to second order

The near-axis description of optimised stellarators, at second order in the expansion, provides important information about the field, both of physical and practical importance for stellarator optimisation. It however remains relatively underdeveloped for an important class of such stellarators, called quasi-isodynamic (QI). In this paper we develop the theoretical and numerical framework for the construction of such solutions. We find that the case of QI stellarators calls for the careful treatment of continuity, smoothness and periodicity of the various functions involved, especially for so-called half-helicity fields, which feature prominently in existing QI designs. The numerical implementation of necessary elements is described, and several examples are constructed and quantitatively verified in detail. This work establishes a basis for further systematic exploration of the space of QI stellarators, and the development of both theoretical and practical tools to facilitate effective optimisation of QI stellarators.