Observation of gravitational waves (GWs) in two different frequency bands is referred to as {\it multiband GW astronomy}. With the planned Laser Interferometric Space Antenna (LISA) operating in the $10^{-4}-0.1$ Hz range, and third generation (3G) ground-based detectors such as the Cosmic Explorer (CE) and Einstein Telescope (ET), operating in the $1$--$10^4$ Hz range, multiband GW astronomy could be a reality in about a decade. In this paper we present the potential of multiband observations of intermediate mass binary black holes (IMBBHs) of component masses ${\sim}10^2$--$10^3M_{\odot}$ to test general relativity (GR). We show that mutiband observations of IMBBHs would permit multiparameter tests of GR---tests where more than one post-Newtonian (PN) coefficient is simultaneously measured yielding more rigorous constraints on possible modifications to GR. We also find that the improvement due to multibanding can often be much larger than the best of the bounds from either of the two observatories. The origin of this result, as we shall demonstrate, can be traced to the lifting of degeneracies among the various parameters when the information from LISA and 3G are taken together. We obtain the best multiband bounds for an IMBBH with a total redshifted mass of $200M_{\odot}$ and a mass ratio of 2. For single-parameter tests, this system at 1 Gpc would allow us to constrain the deviations on all the PN coefficients to below 10\% and derive simultaneous bounds on the first seven PN coefficients to below 50\% (with low spins).

Proposed third generation gravitational wave (GW) interferometers such as Cosmic Explorer will have the sensitivity to observe double neutron star (DNS) mergers up to a redshift of $\sim 5$ with good signal to noise ratios. We argue that the comoving spatial distribution of DNS mergers leaves a unique imprint on the statistical distribution of signal to noise ratios (SNRs) of the detected DNS mergers. Hence the SNR distribution of DNS mergers will facilitate a novel probe of their redshift evolution independent of the luminosity distance measurements. We consider detections of DNS mergers by the third generation detector Cosmic Explorer and study the SNR distribution for different possible redshift evolution models of DNSs and employ Anderson Darling p-value statistic to demonstrate the distinguishability between these different models. We find that a few hundreds of DNS mergers in the Cosmic Explorer era will allow us to distinguish between different models of redshift evolution.

In an earlier work [S. Kastha et al., PRD {\bf 98}, 124033 (2018)], we developed the {\it parametrized multipolar gravitational wave phasing formula} to test general relativity, for the non-spinning compact binaries in quasi-circular orbit. In this paper, we extend the method and include the important effect of spins in the inspiral dynamics. Furthermore, we consider parametric scaling of PN coefficients of the conserved energy for the compact binary, resulting in the parametrized phasing formula for non-precessing spinning compact binaries in quasi-circular orbit. We also compute the projected accuracies with which the second and third generation ground-based gravitational wave detector networks as well as the planned space-based detector LISA will be able to measure the multipole deformation parameters and the binding energy parameters. Based on different source configurations, we find that a network of third-generation detectors would have comparable ability to that of LISA in constraining the conservative and dissipative dynamics of the compact binary systems. This parametrized multipolar waveform would be extremely useful not only in deriving the first upper limits on any deviations of the multipole and the binding energy coefficients from general relativity using the gravitational wave detections, but also for science case studies of next generation gravitational wave detectors.

We propose a novel method to test the consistency of the multipole moments of compact binary systems with the predictions of General Relativity (GR). The multipole moments of a compact binary system, known in terms of symmetric and trace-free tensors, are used to calculate the gravitational waveforms from compact binaries within the post-Newtonian (PN) formalism. For nonspinning compact binaries, we derive the gravitational wave phasing formula, in the frequency domain, parametrizing each PN order term in terms of the multipole moments which contribute to that order. Using GW observations, this {\it{parametrized multipolar phasing}} would allow us to derive the bounds on possible departures from the multipole structure of GR and hence constrain the parameter space of alternative theories of gravity. We compute the projected accuracies with which the second generation ground-based detectors, such as Advanced Laser Interferometer Gravitational-wave Observatory (LIGO), the third generation detectors such as Einstein Telescope and Cosmic Explorer, as well as space-based detector Laser Interferometer Space Antenna (LISA) will be able to measure these multipole parameters. We find that while Advanced LIGO can measure the first two or three multipole coefficients with good accuracy, Cosmic Explorer and Einstein Telescope may be able to measure the first four multipole coefficients which enter the phasing formula. Intermediate mass ratio inspirals, with mass ratio of several tens, in the frequency band of planned space-based LISA mission should be able to measure all the seven multipole coefficients which appear in the 3.5PN phasing formula. Our finding highlights the importance of this class of sources for probing the strong-field gravity regime. The proposed test will facilitate the first probe of the multipolar structure of Einstein's general relativity.

Recent exoplanet surveys have predicted a very large population of planetary systems in our galaxy, more than one planet per star on the average, perhaps totalling about two hundred billion. These surveys, based on electro-magnetic observations, are limited to a very small neighbourhood of the solar system and the estimations rely on the observations of only a few thousand planets. On the other hand, orbital motions of planets around stars are expected to emit gravitational waves (GW), which could provide information about the planets not accessible to electro-magnetic astronomy. The cumulative effect of the planets, with periods ranging from few hours to several years, is expected to create a stochastic GW background (SGWB). We compute the characteristic GW strain of this background based on the observed distribution of planet parameters. We also show that the integrated extragalactic background is comparable or less than the galactic background at different frequencies. Our estimate shows that the net background is significantly below the sensitivities of the proposed GW experiments in different frequency bands. However, we notice that the peak of the spectrum, at around $10^{-5}$Hz, is not too far below the proposed space based GW missions. A future space based mission may be able to observe or tightly constrain this signal, which will possibly be the only way to probe the galactic population of exoplanets as a whole.