Spin-squared Hamiltonian of next-to-leading order gravitational interaction

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 three-dimensional 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.