The planned sensitivity upgrades to the LIGO and Virgo facilities could uniquely identify host galaxies of dark sirens-compact binary coalescences without any electromagnetic counterparts-within a redshift of z = 0.1. This is aided by the higher order spherical harmonic modes present in the gravitational-wave signal, which also improve distance estimation. In conjunction, sensitivity upgrades and higher modes will facilitate an accurate, independent measurement of the host galaxy's redshift in addition to the luminosity distance from the gravitational wave observation to infer the Hubble-Lema\^itre constant H0 to better than a few percent in five years. A possible Voyager upgrade or third generation facilities would further solidify the role of dark sirens for precision cosmology in the future.

The observation of a pulsar closely orbiting the galactic center supermassive black hole would open the window for an accurate determination of the black hole parameters and for new tests of General Relativity. An important relativistic effect which has to be taken into account in the timing model is the propagation delay of the pulses in the gravitational field of the black hole. Due to the extreme mass ratio of the pulsar and the supermassive back hole we use the test particle limit to derive an exact analytical formula for the propagation delay in a Schwarzschild spacetime. We then compare this result to the propagation delays derived in the usually employed post-Newtonian approximation, in particular to the Shapiro delay up to the second post-Newtonian order. For edge-on orbits we also consider modifications of the Shapiro delay which take the lensing effects into account. Our results are then used to assess the accuracy of the different orders of the post-Newtonian approximation of the propagation delay. This comparison indicates that for (nearly) edge-on orbits the new exact delay formula should be used.

We have used optical V and R band observations from the Massive Compact Halo Object (MACHO) project on a sample of 59 quasars behind the Magellanic clouds to study their long term optical flux and colour variations. These quasars lying in the redshift range of 0.2 < z < 2.8 and having apparent V band magnitudes between 16.6 and 20.1 mag have observations ranging from 49 to 1353 epochs spanning over 7.5 years with frequency of sampling between 2 to 10 days. All the quasars show variability during the observing period. The normalized excess variance (Fvar) in V and R bands are in the range 0.2% < Fvar < 1.6% and 0.1% < Fvar < 1.5%. In a large fraction of the sources, Fvar is larger in the V-band compared to the R-band. From the z-transformed discrete cross correlation function analysis, we find that there is no lag between the V and R-band variations. Adopting the Markov Chain Monte Carlo (MCMC) approach, and properly taking into account the correlation between the errors in colours and magnitudes, it is found that majority of the sources show a bluer when brighter trend, while a minor fraction of quasars show the opposite behaviour. This is similar to the results obtained from other two independent algorithms namely the weighted linear least squares fit (FITEXY) and the bivariate correlated errors and intrinsic scatter regression (BCES). However, the ordinary least squares (OLS) fit normally used in the colour variability studies of quasars, indicates that all the quasars studied here show a bluer when brighter trend. It is therefore very clear that OLS algorithm cannot be used for the study of colour variability in quasars.