Inspiraling supermassive black hole binary systems with high orbital eccentricity are important sources for space-based gravitational wave (GW) observatories like the Laser Interferometer Space Antenna (LISA). Eccentricity adds orbital harmonics to the Fourier transform of the GW signal and relativistic pericenter precession leads to a three-way splitting of each harmonic peak. We study the parameter estimation accuracy for such waveforms with different initial eccentricity using the Fisher matrix method and a Monte Carlo sampling of the initial binary orientation. The eccentricity improves the parameter estimation by breaking degeneracies between different parameters. In particular, we find that the source localization precision improves significantly for higher-mass binaries due to eccentricity. The typical sky position errors are $\sim1 $deg for a nonspinning, $10^7\,M_{\odot}$ equal-mass binary at redshift $z=1$, if the initial eccentricity 1 yr before merger is $e_0\sim 0.6$. Pericenter precession does not affect the source localization accuracy significantly, but it does further improve the mass and eccentricity estimation accuracy systematically by a factor of 3--10 for masses between $10^6$ and $10^7\,M_{\odot}$ for $e_0 \sim 0.3$.

We present a closed system of coupled first order differential equations governing the secular linear momentum loss of a compact binary due to emitted gravitational waves, with the leading order relativistic and spin-orbit perturbations included. In order to close the system, the secular evolution equations of the linear momentum derived from the dissipative dynamics are supplemented with the secular evolutions of the coupled angular variables, as derived from the conservative dynamics.

We analyze galactic black hole mergers and their emitted gravitational waves. Such mergers have typically unequal masses with mass ratio of the order 1/10. The emitted gravitational waves carry the inprint of spins and mass quadrupoles of the binary components. Among these contributions, we consider here the quasi-precessional evolution of the spins. A method of taking into account these third post-Newtonian (3PN) effects by renormalizing (redefining) the 1.5 PN and 2PN accurate spin contributions to the accumulated orbital phase is developed.

We give here a new third post-Newtonisn (3PN) spin-spin contribution (in the PN parameter $\epsilon $) to the accumulated orbital phase of a compact binary, arising from the spin-orbit precessional motion of the spins. In the equal mass case this contribution vanishes, but LISA sources of merging supermassive binary black holes have typically a mass ratio of 1:10. For such non-equal masses this 3PN correction is periodic in time, with period approximately $\epsilon ^{-1}$ times larger than the period of gravitational waves. We derive a renormalized and simpler expression of the spin-spin coefficient at 2PN, as an average over the time-scale of this period of the combined 2PN and 3PN contribution. We also find that for LISA sources the quadrupole-monopole contribution to the phase dominates over the spin-spin contribution, while the self-spin contribution is negligible even for the dominant spin. Finally we define a renormalized total spin coefficient $\bar{\sigma}$ to be employed in the search for gravitational waves emitted by LISA sources.

The radial component of the motion of compact binary systems composed of neutron stars and/or black holes on eccentric orbit is integrated. We consider all type of perturbations that emerge up to second post-Newtonian order. These perturbations are either of relativistic origin or are related to the spin, mass quadrupole and magnetic dipole moments of the binary components. We derive a generalized Kepler equation and investigate its domain of validity, in which it properly describes the radial motion.

In this work we consider compact binaries on eccentric orbit under the spin-spin interaction. Using the post-Newtonian formalism, the binaries undergo a perturbed Keplerian motion. Here we investigate only the radial motion and derive the contribution of the second post-Newtonian order spin-spin effect to the solution.

The angular and the radial parts of the dynamics of the perturbed Kepler motion are separable in many important cases. In this paper we study the radial motion and its parametrizations. We develop in detail a generalized eccentric anomaly parametrization and a procedure of computing a generic class of integrals based on the residue theorem. We apply the formalism to determine various contributions to the luminosity of a compact binary.

We present here a complete list of contributions to the gravitational wave frequency evolution from compact binaries on circular orbits up to the second post-Newtonian order by including the interaction of magnetic dipole moments, quadrupole-monopole interactions together with the spin-orbit, the spin-spin and the self-interaction spin contributions. We apply our results to the Manchester et al. model of the J0737-3039A-B double pulsar.

Binary systems subject to generic perturbations evolve on quasiperiodic orbits. We derive the most generic class of perturbations, which allow to evaluate secular effects via generalized complex true and eccentric anomaly parameters, by use of the residue theorem. Such perturbations include both the generic Brumberg force and various linear contributions to the conservative dynamics of compact binaries, up to the second post-Newtonian order. As an example, we compute the self-spin contribution to the luminosity due to gravitational radiation of a compact binary consisting of galactic black holes, the most important type of source for LISA. This translates to simply calculate the residue of the instantaneous energy loss in the origin of the complex parameter plane, which illustrates the power of the presented method.

Compact binaries consisting of neutron stars / black holes on eccentric orbit undergo a perturbed Keplerian motion. The perturbations are either of relativistic origin or are related to the spin, mass quadrupole and magnetic dipole moments of the binary components. The post-Newtonian motion of such systems decouples into radial and angular parts. We present here for the first time the radial motion of such a binary encoded in a generalized Kepler equation, with the inclusion of all above-mentioned contributions, up to linear order in the perturbations. Together with suitably introduced parametrizations, the radial motion is solved completely.

Gravitational radiation drives compact binaries through an inspiral phase towards a final coalescence. For binaries with\textit{spin, mass quadrupole and magnetic dipole moments}, various contributions add to this process, which is characterized by the rate of increase $df/dt$ of the gravitational wave frequency and the accumulated number $\mathcal{N}$ of gravitational wave cycles. We present here all contributions to $df/dt$ and $% \mathcal{N}$ up to the second post-Newtonian order. Among them we give for the first time the contributions due to the \textit{self-interaction of individual spins}. These are shown to be commensurable with the proper spin-spin contributions for the recently discovered J0737-3039 double pulsar, and argued to represent the first corrections to the radiation reaction in the Lense-Thirring approach.