We characterize and analyze rotational torsional oscillations developing in a large-eddy magnetohydrodynamical simulation of solar convection (Ghizaru, Charbonneau, and Smolarkiewicz, Astrophys. J. Lett., 715, L133 (2010); Racine et al., Astrophys. J., 735, 46 (2011)) producing an axisymmetric large-scale magnetic field undergoing periodic polarity reversals. Motivated by the many solar-like features exhibited by these oscillations, we carry out an analysis of the large-scale zonal dynamics. We demonstrate that simulated torsional oscillations are not driven primarily by the periodically-varying large-scale magnetic torque, as one might have expected, but rather via the magnetic modulation of angular-momentum transport by the large-scale meridional flow. This result is confirmed by a straightforward energy analysis. We also detect a fairly sharp transition in rotational dynamics taking place as one moves from the base of the convecting layers to the base of the thin tachocline-like shear layer formed in the stably stratified fluid layers immediately below. We conclude by discussing the implications of our analyses with regards to the mechanism of amplitude saturation in the global dynamo operating in the simulation, and speculate on the possible precursor value of torsional oscillations for the forecast of solar cycle characteristics.

Using a Legendre transformation, we compute the unconstrained Hamiltonian of a spinning test-particle in a curved spacetime at linear order in the particle spin. The equations of motion of this unconstrained Hamiltonian coincide with the Mathisson-Papapetrou-Pirani equations. We then use the formalism of Dirac brackets to derive the constrained Hamiltonian and the corresponding phase-space algebra in the Newton-Wigner spin supplementary condition (SSC), suitably generalized to curved spacetime, and find that the phase-space algebra (q,p,S) is canonical at linear order in the particle spin. We provide explicit expressions for this Hamiltonian in a spherically symmetric spacetime, both in isotropic and spherical coordinates, and in the Kerr spacetime in Boyer-Lindquist coordinates. Furthermore, we find that our Hamiltonian, when expanded in Post-Newtonian (PN) orders, agrees with the Arnowitt-Deser-Misner (ADM) canonical Hamiltonian computed in PN theory in the test-particle limit. Notably, we recover the known spin-orbit couplings through 2.5PN order and the spin-spin couplings of type S_Kerr S (and S_Kerr^2) through 3PN order, S_Kerr being the spin of the Kerr spacetime. Our method allows one to compute the PN Hamiltonian at any order, in the test-particle limit and at linear order in the particle spin. As an application we compute it at 3.5PN order.

We compute the flux of linear momentum carried by gravitational waves emitted from spinning binary black holes at 2PN order for generic orbits. In particular we provide explicit expressions of three new types of terms, namely next-to-leading order spin-orbit terms at 1.5 PN order, spin-orbit tail terms at 2PN order, and spin-spin terms at 2PN order. Restricting ourselves to quasi-circular orbits, we integrate the linear momentum flux over time to obtain the recoil velocity as function of orbital frequency. We find that in the so-called superkick configuration the higher-order spin corrections can increase the recoil velocity up to about a factor 3 with respect to the leading-order PN prediction. Furthermore, we provide expressions valid for generic orbits, and accurate at 2PN order, for the energy and angular momentum carried by gravitational waves emitted from spinning binary black holes. Specializing to quasi-circular orbits we compute the spin-spin terms at 2PN order in the expression for the evolution of the orbital frequency and found agreement with Mik\'oczi, Vas\'uth and Gergely. We also verified that in the limit of extreme mass ratio our expressions for the energy and angular momentum fluxes match the ones of Tagoshi, Shibata, Tanaka and Sasaki obtained in the context of black hole perturbation theory.

We analyze in detail the spin precession equations in binary black hole systems, when the tidal torque on a Kerr black hole due to quadrupole-monopole coupling is taken into account. We show that completing the precession equations with this term reveals the existence of a conserved quantity at 2PN order when averaging over orbital motion. This quantity allows one to solve the (orbit-averaged) precession equations exactly in the case of equal masses and arbitrary spins, neglecting radiation reaction. For unequal masses, an exact solution does not exist in closed form, but we are still able to derive accurate approximate analytic solutions. We also show how to incorporate radiation reaction effects into our analytic solutions adiabatically, and compare the results to solutions obtained numerically. For various configurations of the binary, the relative difference in the accumulated orbital phase computed using our analytic solutions versus a full numerical solution vary from about 0.3% to 1.8% over the 80 - 140 orbital cycles accumulated while sweeping over the orbital frequency range 20 - 300 Hz. This typically corresponds to a discrepancy of order 5-6 radians. While this may not be accurate enough for implementation in LIGO template banks, we still believe that our new solutions are potentially quite useful for comparing numerical relativity simulations of spinning binary black hole systems with post-Newtonian theory. They can also be used to gain more understanding of precession effects, with potential application to the gravitational recoil problem, and to provide semi-analytical templates for spinning, precessing binaries.

We study a non-dissipative hydrodynamical mechanism that can stabilize the spin of the accretor in an ultra-compact double white dwarf binary. This novel synchronization mechanism relies on a nonlinear wave interaction spinning down the background star. The essential physics of the synchronization mechanism is summarized as follows. As the compact binary coalesces due to gravitational wave emission, the largest star eventually fills its Roche lobe and accretion starts. The accretor then spins up due to infalling material and eventually reaches a spin frequency where a normal mode of the star is resonantly driven by the gravitational tidal field of the companion. If the resonating mode satisfies a set of specific criteria, which we elucidate in this paper, it exchanges angular momentum with the background star at a rate such that the spin of the accretor locks at this resonant frequency, even though accretion is ongoing. Some of the accreted angular momentum that would otherwise spin up the accretor is fed back into the orbit through this resonant tidal interaction. Two modes capable of stabilizing the accretor's spin are the l=4,m=2 and l=5,m=3 CFS unstable hybrid r-modes, which stabilize the spin of the accretor at frequency 2.6 and 1.5 times the binary's orbital frequency respectively. Since the stabilization mechanism relies on continuously driving a mode at resonance, its lifetime is limited since eventually the mode amplitude saturates due to non-linear mode-mode coupling. Rough estimates of the lifetime of the effect lie from a few orbits to millions of years.