Within the post Newtonian framework the fully reduced Hamiltonian (i.e., with eliminated spin supplementary condition) for the next-to-leading order spin-squared dynamics of general compact binaries is presented. The Hamiltonian is applicable to the spin dynamics of all kinds of binaries with self-gravitating components like black holes and/or neutron stars taking into account spin-induced quadrupolar deformation effects in second post-Newtonian order perturbation theory of Einstein's field equations. The corresponding equations of motion for spin, position and momentum variables are given in terms of canonical Poisson brackets. Comparison with a nonreduced potential calculated within the Effective Field Theory approach is made.

The static, i.e., linear momentum independent, part of the next-to-leading order (NLO) gravitational spin(1)-spin(1) interaction Hamiltonian within the post-Newtonian (PN) approximation is calculated from a 3-dim. covariant ansatz for the Hamilton constraint. All coefficients in this ansatz can be uniquely fixed for black holes. The resulting Hamiltonian fits into the canonical formalism of Arnowitt, Deser, and Misner (ADM) and is given in their transverse-traceless (ADMTT) gauge. This completes the recent result for the momentum dependent part of the NLO spin(1)-spin(1) ADM Hamiltonian for binary black holes (BBH). Thus, all PN NLO effects up to quadratic order in spin for BBH are now given in Hamiltonian form in the ADMTT gauge. The equations of motion resulting from this Hamiltonian are an important step toward more accurate calculations of templates for gravitational waves.

Based on recent developments by the authors a next-to-leading order spin(1)-spin(2) Hamiltonian is derived for the first time. The result is obtained within the canonical formalism of Arnowitt, Deser, and Misner (ADM) utilizing their generalized isotropic coordinates. A comparison with other methods is given.